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APRESENTAÇAO 


Estes apontamentos foram sendo escritos ao longo dos semestres em que dei Análise Com¬ 
plexa e Equações Diferenciais no iST/TagusPark. Na terceira iteração, decidi separar a 
análise complexa das equações diferenciais. Nesta (a quarta), acrescentei alguns tópicos, 
figuras e correções. Poderão encontrar informação atualizada em http : //web .técnico 
ulisboa.pt/joao .boavida/2015/ACED/, 

Quis que o resultado final ficasse relativamente curto, pois é fácil ser-se inundado com 
imensa informação e perder de vista as ideias centrais. Assim, tentei evitar detalhes exagera¬ 
dos e deixei alguns tópicos para os exercícios (isto não significa que esses tópicos são menos 
importantes; só significa que explicar os detalhes por escrito não ia acrescentar muito àquilo 
que conseguiriam fazer sozinhos, seguindo os passos sugeridos). 

A falta de melhor nome, chamei “episódio” a cada conjunto de apontamentos, corres¬ 
pondente a cerca de uma semana de aulas. Tentei escrever informalmente, como se estivés¬ 
semos a conversar (daí usar o plural da primeira pessoa). Como não é um monólogo, como 
a ideia é que participem ativamente na viagem, deixei algumas perguntas e exercícios pelo 
caminho. Quando algum exercício for inesperadamente exigente, eu assinalo-o com *. 

Ao contrário de análise complexa (em que as ideias centrais são tão surpreendentes que 
se justifica começar a explorá-las desde logo, em vez de só as descobrir mesmo no final), em 
equações diferenciais, em geral, a ordem “lógica” coincide com a ordem natural, e por isso 
eu vou normalmente segui-la. Os dois primeiros episódios são a única exceção relevante: 
a ordem “lógica” seria começar pelas demonstrações no episódio [2j enquanto que a ordem 
mais intuitiva é começar pelos exemplos no episódio [T[ 

Nada do que aqui lerem é original; estas ideias já estavam arrumadas, aproximadamente 
no formato atual, há mais de meio século. 

Uma equação diferencial é uma equação cuja incógnita é uma função y(t) e em que apa¬ 
recem derivadas de y. Por exemplo, 


y'-t-y = o 

é uma equação diferencial de primeira ordem (a “ordem” é a ordem da derivada mais alta 
que aparece). Substituindo y(t) = t 2 , podemos verificar que não é uma solução, pois 

(t 2 y-t-(t 2 )# 0. 

Por outro lado, y(t) = 0 é uma solução, pois 

(O)' - t • 0 = 0. 

E comum omitir a variável independente na equação (ou seja, escrever y em vez de y(t)). E 
claro que isso não causa ambiguidade, pois o facto de escrevermos y' sugere imediatamente 
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que y depende da outra variável, e não o contrário. Por vezes, usamos a notação mais 
explícita Tjj, para indicar claramente qual a variável dependente e a variável independente. 

Também distinguimos entre equações diferenciais ordinárias (EDOs, as que só têm uma 
variável independente) e as equações diferenciais parciais (EDPs, as com mais que uma va¬ 
riável independente, e que portanto envolvem derivadas parciais). Vamos trabalhar quase 
só com EDOs; só no último episódio é que faremos alguma coisa com EDPs. 

A primeira vista, não é nada óbvio como encontrar soluções para uma equação diferen¬ 
cial. E verdade que algumas equações diferenciais podem ser resolvidas de forma relativa¬ 
mente simples; para dar um exemplo extremo, podemos usar integração para resolver 

/ = t 2 , com y(0) = 2. 

Aqui, y(0) (ou y noutro ponto específico, se for isso que temos) é a condição inicial. 

Mas esse método só raramente funciona. Por exemplo, não funciona com a equação 

xy'-2yy' + y = 0, y( 0) = 3. 

O episódio [T] é dedicado a várias estratégias para resolver EDOs de primeira ordem. 

E razoável querer um critério simples para determinar a priori se uma EDO tem solução, 
e se, fixada a condição inicial, a solução é única. E que em geral esperamos que, se dois 
carros começam a viagem lado a lado e as suas velocidades variam da mesma forma ao 
longo do tempo, então deveriam continuar lado a lado ao longo de toda a viagem. No 
episódio [2] discutimos um critério. 

Se tivermos um sistema como 

I *í = x 2 , 

| x' 2 = —42x 1 + 13x 2 , 

podemos escrevê-lo na forma matricial 


*1 

/ 

0 

1 ' 


*1 

. X 2 . 


-42 

13 


. X 2 . 


A 


obtendo x' = Ax. Imaginemos que x = Su é uma mudança de variável (exigimos em 
particular que a matriz S seja invertível), permitindo escrever Su 7 = ASu. Se escolhermos 
J = S~ 1 AS, temos SJ =AS, ou seja, a equação fica Su 7 = SJu ou (multiplicando à esquerda 
por S _1 , de ambos os lados, e simplificando) u 7 = Ju. Acontece que todas as matrizes têm 
uma forma canónica de Jordan J, para a qual u 7 = Ju é fácil de resolver. A partir dessa 
solução u, podemos obter a solução x = Su. No episódio [3] vamos discutir essas e outras 
EDOs lineares em M n . 

Uma EDO de ordem superior pode ser transformada numa EDO em M. n . Por exemplo, 
f" — 13/ 7 + 42/ = 0 (a variável independente continua a ser t) pode ser simulada usando 
x 1 = f e x 2 = f. O sistema fica então 


/ *í =f' = x 2 , 

*i 

/ 

0 1 “ 


*i 

{ , , OU 


-42 13 


{ x' 2 = f" = -42/ + 13/ = -42*! + 13*2, 

L X 2 . 



. X 2 J 



A 
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Obtemos precisamente a mesma equação de há pouco, o que significa que os métodos para 
equações lineares em M n são suficientes para resolver equações de ordem superior. Porém, 
raramente são convenientes, e por isso no episódio [4] vamos discutir outros métodos, com 
especial ênfase no método dos aniquiladores. 

Entre outras propriedades, a transformada de Laplace converte soluções de EDOs de 
ordem superior em funções no plano complexo com um número finito de polos, e é possível 
usar o teorema dos resíduos para recuperar a função original. Isto só por si não seria muito 
interessante. O que é interessante é que a transformada converte equações diferenciais em 
equações algébricas, mais fáceis de resolver. Vamos falar de tudo isto no episódio [õ] 
Finalmente, imaginem que querem resolver esta equação diferencial parcial: 

du _ d 2 u 
dt dx 2 ’ 

com 0 < r < 1 e t > 0. Como é que alguém poderia resolver uma equação desse género? 
Na falta de outras ideias, uma opção é tentar encontrar soluções por tentativa e erro. Por 
exemplo, podemos substituir u(t,x) = T(t) -X(x). Com algum cuidado, acabamos por 
concluir que 

X(x) = Ci cos(7tfcx) + C 2 sin(7tfcx), onde k e 7L. 

De facto, para as condições de fronteira (isto é, condições para u ou em x = 0 e x = 1) 
que vamos considerar, é possível mostrar que a solução geral é da forma 

u(t,x) = E a fc (t) cos(ttícx) + ^ b fc (t) sin(7iícx), 
fc>o k>0 

onde as funções a k e b k podem ser obtidas (com as técnicas dos episódios anteriores) a 
partir da condição inicial iz(0, x). No episódio [ój vamos falar das séries de Fourier (isto 
é, séries semelhantes à de iz(t,x)) e algumas variantes. No episódio [ 7 ] discutimos vários 
exemplos de EDPs, resolvidas aplicando técnicas de todos os episódios anteriores. 


João Pedro Boavida <joao . boavidaStecnico .ulisboa.pt>, 
Departamento de Matemática, Instituto Superior Técnico, 
Universidade de Lisboa, 

Outubro de 2015 
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EPISÓDIO 1 


EQUAÇÕES DIFERENCIAIS ORDINÁRIAS 


Neste episódio, vamos discutir as estratégias mais importantes para resolver equações dife¬ 
renciais ordinárias (EDOs) de primeira ordem. Aqui, mencionamos só algumas ideias. 
Considerem a equação 

xy' - 2yy' + y = 0, y(0) = 3. 

À primeira vista, não é óbvio como resolvê-la. Por outro lado, quem é j^(xy — y 2 )? (Notem 
que x é a variável independente, por isso (y 2 )' = ^;(y 2 ) não é só 2y!) Qual é então a 
solução da equação? 

Uma equação diferencial diz-se exata se existe uma função 4>(x, y) que permite escrever 
a equação na forma ^($(x,y)) = 0. Usando a regra da função composta, isso fica + 
^ = 0. Assim, se y(x) é uma solução da equação, $(x,y(x)) é constante. 

Pode acontecer que uma equação não seja exata, mas que seja redutível a exata. Por 
exemplo, a equação 

xe x y' + e x y - 1 = 0 

não é exata (como é que podemos ter a certeza?). Mas se multiplicarmos tudo por e~ x 
(porque é que podemos?), obtemos 

xy' + y-e~ x = 0, 

que é exata. Conseguem resolvê-la? 


(1.1) Alguns exemplos para ilustrar como é possível resolver certas equações, mesmo sem 
grande estratégia. Na equação 

x + xy 7 = 2y, 

podemos tentar que xy' e 2y sejam parcelas de uma derivada do quociente (porquê?). Para 
tal, precisaríamos que fossem x 2 y / e 2xy, algo que podemos conseguir multiplicando tudo 
por x e fazendo algumas manipulações simples: 

2 Ir 2V 2 

x y (x j y = -X z . 


Daqui, basta-nos dividir tudo por (x 2 ) 2 , que nos leva a 

(x 2 )/-(x 2 yy 1 


(x 2 ) 2 


e a 


&)'=&' 


ou 


y 1 

Á = - + c > 

x z x 


para alguma constante C. Assim, desde que o intervalo de definição que queremos exclua 
x = 0 (porquê?), concluímos que y = x + Cx 2 . 
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Episódio 1: Equações diferenciais ordinárias 


(1.2) Vejamos um exemplo com uma equação de segunda ordem (i.e., envolvendo a segunda 
derivada). Se 


(2x + 1 )/" - 2f = (1 - 4x 2 ) (2x + 1), 


é tentador ver o lado esquerdo como numerador de uma regra do quociente. Se dividirmos 
tudo por (2x + l) 2 , ficamos com 

(2x + l)/ // -(2x + l)7 / 1 — 4x 2 

(2x + l) 2 ~~ 2x + 1 

Como 1 — 4x 2 = (1 — 2x)(l + 2x), podemos simplificar um pouco mais e chegar a 


( 


f 

2x + 1 


)' 


= 1 - 2x. 


Assim, para intervalos de definição excluindo 2x + 1 ^ 0 (porquê?), vemos que 


f 

2x + 1 


= x - x 2 +A, 


para alguma constante A. Escrevendo em ordem a f'e primitivando uma segunda vez, 
concluímos que 

x 2 x 3 X 4 „ 

/ = y + y-y +A(x 2 + x) + B, 

para uma segunda constante B. 


(1.3) Exercício. Considerem uma solução definida em todo o K. Acabámos de ver que para 
x < — | (porquê — |?) a solução é da forma 

r 2 r 3 Y 4 

f = y + 7- 2 +Al + + Bl 

para constantes A 1 e B 1 , e que para x > — | é da forma 

2 3 4 

-y* V ’ 

/ = y + “^-y +^2 (x 2 + x) + B 2 

para constantes A 2 e B 2 . Quais as relações entre as várias constantes? Conseguem determi¬ 
nar todas as soluções da equação? Pode ser útil voltar à equação original e tentar determinar 

n-li 


(1.4) Exercício. Resolvam estas equações diferenciais tentando reconhecer regras de deri¬ 
vação. Determinem também os intervalos em que as soluções são válidas. (Pode ser útil 
começarem por identificar a variável dependente e a independente, já que variam de equa¬ 
ção para equação.) Se tiverem informação suficiente, determinem as constantes. 

(a) y' = t 2 , y(0) = 2; 

(b) 2 yy' = x 2 , y(l) = l; 

(c) 2 yy' = x 2 , y( 1) = -1; 

(d) f + tf' = t, f (—2) = —1; 

/( 2 ) = — 2 ; 


(e) f + tf' = t, 
1 


u 

(f) - = 


u(l) = -l. 
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(1.5) As equações separáveis são aquelas em que, através de multiplicações ou divisões 
(de ambos os membros pela mesma função), conseguimos separar as variáveis em lados 
opostos da equação. Por exemplo, 


d y 2 

—-x = y • logx, 
dx 


y(i) = -2, 


pode ser escrita na forma 


1 dy logx 
y 2 dx x 


y( i) = -2, 


desde que não haja divisões por zero. (Já agora, log é o logaritmo natural, já que não há 
qualquer outro logaritmo relevante daqui em diante.) Ou, escrito de forma mais sugestiva, 

dy 1 

—y = — • log x dx. 
y x 

Integrando de ambos os lados, obtemos (é preciso fazer mudanças de variável—quais?) 

1 _ (logx) 2 ^ 

y 2 

(não se esqueçam da constante!). Substituindo a condição inicial y(l) = —2, descobrimos 
que C = 1/2. Portanto, 

2 

(logx) 2 + l 


(1.6) Exercício. Aqui está uma equação mais simples: usando a mesma estratégia, procurem 
soluções de 

y'-xy + y = 0, y(0) = 2. 

Notem que é preciso algum cuidado a lidar com os casos y = 0 ou y < 0. 


(1.7) Exercício. Resolvam a equação 


yV = l + y, y(0) = e — 1. 


(1.8) Exercício. Resolvam a equação 

ty / + y 2 + l =y' + t 2 y 2 + t 2 , y(0) = 0. 

(1.9) Exercício. Considerem agora a equação 

/ = (3 + t)y 2 . 

(a) Assumindo que y(t) ^ 0, qual a solução desta equação? Notem que não estamos 
a especificar a sua condição inicial. Isso significa que, depois de fazer as primitivas, não 
temos como determinar o valor da constante, e por isso deixamo-la ficar, sem tentarmos 
especificar o seu valor. 

(b) Qual o intervalo de definição da solução da alínea anterior? 
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Episódio 1: Equações diferenciais ordinárias 


(c) Encontrem uma solução com y(0) = 0. 

(d) Imaginem que y(t) é uma solução com y(t 0 ) > 0- Mostrem que, no seu intervalo 
de definição, temos sempre y(t) > 0. 

(e) Façam o mesmo para uma solução com y(t 0 ) < 0. 

(f) Concluam que só existe uma solução com y(t 0 ) = 0- Qual é? 

(1.10) Exercício. Façam um estudo (como o do exercício anterior) das soluções de cada 
uma destas equações: 

(a) (3 +t)y' = y; (b) (3+t)y / = y 2 ; (c)/-yt = /t 2 . 


(1.11) Equações exatas. Há pouco, considerámos a equação 


xy'-2yy' + y = 0 


(a variável independente é x) com condição inicial y(0) = 3. Observámos que 

^ {xy - y 2 ) = xy' - 2 yy’ + y, 

o que nos leva a concluir (porquê?) que xy — y 2 é constante. Substituindo a condição 
inicial, vemos que a constante é —9 e obtemos a fórmula xy — y 2 = —9, que pode ser 
resolvida para obter y. 

Dissemos que a equação (com variável independente x) 


P + Q.y = o 


é exata se existir uma função 3>(x,y) cujo gradiente seja (P, Q), ou seja, cujas derivadas 
parciais sejam = P e = Q. Nesse caso, 3>(x,y) é constante ao longo das soluções 
y(x) da equação (notem que x é a variável independente), pois a regra da função composta 
diz-nos que 


d 

dx 


($(x,y(x))) 




dy 



= P + Q-y' = 0. 


Se um tal $ existir, então teremos 


d 3<í> 
dx dy 
Q 


d 3<í> 
dy dx 
p 


dQ dP 
dx dy 


(Se isto vos lembra alguma coisa, não é coincidência.) Ou seja, a equação é exata apenas se 

ÕQ _ dP 
dx dy ' 
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(1.12) Por exemplo, considerem a equação diferencial 

- 2 x 2 y y' -2xy 2 + e x = 0 
Q p 

com condição inicial y(l) = 2. Queremos encontrar uma função $ cujo gradiente (|^, |y) 
seja (P, Q). Para isso consideramos 

—2 xy 2 + e x => 

- 2 x 2 y => 

= —x 2 y 2 +e x é compatível com ambas as condições. Como 4>(1,2) = e—4 
e o 2 ?), concluímos que a solução satisfaz a equação 

—x 2 y 2 + e x = e — 4. 

Resolvendo em ordem a y, obtemos 

e x — e + 4 
x 2 

Só o sinal positivo é compatível com a condição inicial; essa é a fórmula para a solução. 

(1.13) Exercício. Experimentem resolver estas equações exatas (nalguns casos pode não ser 
prático resolver em ordem a y; vão tão longe quanto conseguirem): 

(a) 2 x + 2yy' = 1 , y( 0 ) = l; 

(b) 2yy' + y sinx = y' cosx, y( 0 ) =- 2 ; 

(c) y 2 + 2 xyy' + y sinx = 3 /cosx — 2 x, y{n) = 0 ; 

(d) 2xyy' + y 2 = 4x 3 , y(l) = l; 

(e) x 4 y' + 4x 3 y + 1 = 2xyy / + y 2 + y', y(l) = 2; 

(f) ny sin( 7 tx) + y / = x 2 y' 4- 2xy + y' cos(trx), y(2) = 1; 

(g) 4x 3 + xy' + 2yy' + y + 9 = 2x 2 y' + 3x 2 + 4xy + 4x, y(0) = 3. 


(1.14) Equações redutíveis a exatas. Também vimos que a equação 

xe x y' + e x y - 1 = 0 

não é exata e observámos que, se a multiplicarmos por e~ x , obtemos 

xy' + y-e~ x = 0, 

que é uma equação exata com as mesmas soluções que a primeira. (Conseguiram resolvê- 
-la?) Por isso, dizemos que a equação original é redutível a exata. A função e~ x usada para 
fazer a redução é chamada/ator integrante ou fator de integração. (Como usamos integração 



4> = —x 2 y 2 + e x + (função de y), 
$ = —x 2 y 2 + (função de x). 


54 > 

dx 

5 $ 

5y 

A escolha 4>(x,y) 
(de onde vêm o 1 
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Episódio 1: Equações diferenciais ordinárias 


para resolver uma equação exata, também há quem fale em “integrar a equação” em vez de 
“resolver a equação”—e por isso o fator ajuda a integrar a equação. Daí o nome.) 

E como é que encontramos fatores integrantes? Boa pergunta. Se tivermos uma equação 


P + Q-y' = 0 


e pi(x, y) for um fator integrante (é tradicional chamar pi ao fator), então a equação 

pi • P + pi • Q ■ y 7 = 0 


é exata. Isso significa que as derivadas cruzadas são iguais, ou seja, 


d(/iQ) _ djpP) 
dx dy 

Simplificando estas derivadas (com a regra do produto), obtemos 


d ii dQ <9pi dP 

x~ ’Q + P’ x~ = -X— - P + H- — ■ 

d x dx ay dy 


Notem que tanto P e Q (conhecidas) como pi (desconhecido) são funções de x e de y. Se 
calhar isto não é lá muito prático... 


(1.15) De acordo: se não soubermos nada sobre pi, não vamos longe. Mas podemos tentar 
casos especiais. Por exemplo, se pi só depender de x, então a equação anterior fica apenas 


dpi 

dx 


dQ dP 

■Q + P- x- = P- x~- 
dx dy 


Se, após simplificação, for possível eliminar y, então esta equação é separável (!) para pi 
(se tal não for possível, então não existe pi dependente só de x). Voltemos ao exemplo de 
há pouco: 

x_e^y' + e x y -1 = 0 

A equação para pi fica 

dpi 

— • x e x +pi • (e x + x e x ) = pi 
dx v -'v^ <-- v -> 

Q SQ 

dx 



Concordo que isto parece cada vez pior. Mas podemos fazer alguns cancelamentos e simpli¬ 
ficar um pouco a equação: 


dpi 

dx 


= -pi. 


Resolvam em ordem a pi (há mais que uma solução, escolham uma qualquer que não seja 
nula, por exemplo, a que tem pi(0) = 1). Que função encontram? (A sério; parem aqui para 
calcular.) E qual é o fator integrante que tínhamos usado antes? Coincidência? 
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(1.16) Se quisermos um fator integrante p dependente apenas de y, a equação no final de 
d 1.1 4D reduz-se (porquê?) a 


Por exemplo, a equação 



d/i 

— -P + /Í- 
dy 


dP 

dy' 


x 4 y 2 y' + 4x 3 y 3 = 2y 3 y', y(0) = -4 


tem um fator integrante que só depende de y. Se escreverem a equação diferencial para p 
(experimentem!) e simplificarem, chegam a algo como 

d/i 2 /i 

dy y ' 

(Já agora, reparem que depois de simplificada a equação envolve apenas p e y. Se ainda 
restasse algum x, era impossível resolvê-la com um p dependente apenas de y. Percebem 
porquê?) Esta equação é separável; qual é a sua solução? 

Podemos ver que p(y) = y ~ 2 é um fator integrante. (Conseguem, sem quaisquer contas, 
indicar outros fatores integrantes?) Multiplicando-o na equação original, obtemos 

x 4 y' + 4x 3 y - 2yy' = 0, y(0) = -4, 

cuja integração leva a x 4 y — y 2 = —16. Resolvendo em ordem a y, obtemos 

x 4 - \j x 8 + 64 


(reparem que escolhemos o sinal negativo, pois só esse é compatível com y(0) = —4). 


(1.17) Exercício. Encontrem fatores integrantes para as seguintes equações, e resolvam-nas. 

(a) x 2 ye x + x 2 e x y / + xy / = y, y(l) = e + l; 

(b) y 2 cosx + yy / sinx = 0, y(|) = 2; 

(c) 2x 2 y 2 y , + xy 3 + y / = 3xyy' + y 2 , y(l) = 1; 

xy' 

(d) yy + y e x + 1 = -, y(0) = 2. 

y 

As equações podem ser resolvidas com fatores da forma p(x) ou p(y). Se houver fatores de 
ambos os tipos, determinem-nos, e resolvam a equação original usando ambos os fatores. 


(1.18) Mas sejamos um pouco mais ambiciosos: haverá alguma forma sistemática de verifi¬ 
car se existe um fator integrante da forma p(x + y) ou p(xy) ou p(x — y)? (Aqui estamos 
a cometer um pequeno abuso de linguagem: em vez de falar em p(x + y), deveríamos falar 
em p(x,y) = v(x + y), onde v é uma nova função. Mas tal não seria necessariamente mais 
claro.) 

Há. Se o fator integrante for da forma p(s) (onde s = x + y, ou s = xy, ou s = x — y, 
ou outra expressão que queiramos experimentar), então a equação no final de Q1.14P pode 


dQ 


d p 


ds 


escrever-se 


dp ds 


dP 
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Episódio 1: Equações diferenciais ordinárias 


Separando o /i à esquerda e abreviando p! = podemos ainda obter (verifiquem) 

/ SQ _ ÕP_ 

/I õx õy 


Se for possível simplificar o lado direito de forma a envolver apenas s, então a solução desta 
equação é um fator integrante. De contrário, a equação é impossível (pois o lado esquerdo 
poderia ser escrito apenas em termos de s e o lado direito não, o que seria uma contradição). 

(1.19) Por exemplo, considerem a equação 

x 5 / + x 4 yy' + 4x 4 y + 4x 3 y 2 = 2x 2 y y’ + 2x y 2 y’ + xy 2 + y 3 , y(l) = 2, 

que tem um fator integrante da forma /i(x + y) (portanto, s = x + y, |^ = le|^ = l). 

Como P = 4x 4 y + 4x 3 y 2 — xy 2 — y 3 e Q = x 5 + x 4 y — 2x 2 y — 2xy 2 , a fórmula para 
— fica 

fj,' x 4 - 4x 3 y - 2xy+ y 2 x 4 - 4x 3 y - 2xy+y 2 

Ijl il.Q — il.p x 5 — 3x 4 y — 2x 2 y — xy 2 — 4x 3 y 2 + y 3 

õx 3 y 

Se existir um fator integrante que dependa apenas de s = x + y, devíamos ser capazes de 
simplificar esta expressão de forma a depender apenas de s. A tarefa não parece promete¬ 
dora, mas se for possível, o x 4 de cima deve (junto com o resto do numerador) conseguir 
cancelar parte do denominador, de forma a obter uma expressão com x + y. Assim, um 
palpite razoável neste caso é: será que o denominador é apenas x + y multiplicado pelo 
numerador? Experimentemos: 

(x + y ) (* 4 - 4x 3 y - 2xy + y 2 ) = x 5 - 4x 4 y - 2x 2 y + xy 2 + x 4 y - 4x 3 y 2 - 2xy 2 + y 3 . 
É realmente o denominador. Assim, concluímos que 

\i! _ 1 1 

fj, x + y s 

Uma solução desta equação é pi(s) = s _1 , mostrando que l/(x +y) é um fator integrante. 
Dito de outra forma, devemos conseguir pôr (x + y) em evidência na equação original 

x 5 y' + x 4 yy' + 4x 4 y + 4x 3 y 2 = 2x 2 yy' + 2 xy 2 y' + xy 2 + y 3 , y(l) = 2. 


Tentemos: 

(x + y) x 4 y' + 4(x + y) x 3 y = 2(x + y) xyy' + (x + y ) y 2 . 

Realmente é verdade. Dividindo tudo por x + y obtemos então a equação exata 

x 4 y / + 4x 3 y-2xyy / -y 2 = 0, y(l) = 2, 

que, por integração, nos leva a x 4 y — xy 2 = —2. Resolvendo em ordem a y e usando a 
condição inicial para decidir o sinal, chegamos a 

x 4 + \/ x 8 + 8x 
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(1.20) Mais um exemplo. A fórmula Ql . 18D para a equação 

2xy 2 y' + 2y 3 = 4x 2 y, y( 1) = 1, 


fica (depois de simplificada) 

\il -4 y 2 + 4x 2 

Se s = x, obteríamos 

/j,' -4y 2 + 4x 2 

M 2 xy 2 

Veem alguma forma de simplificar isto de forma a eliminar o y? Eu também não. 

Com s = y, teríamos 

n' -4 y 2 + 4x 2 
/i —2y 3 + 4x 2 y 

Infelizmente, por pouco (bastava que fosse um 4 em vez de um 2), também não consegui¬ 
mos livrar-nos do x. 

Consideremos então s = y/x, com = —y/x 2 e = l/x. (Em geral, está longe de 
ser evidente qual um s que funcionará. Não se tratando de s = x nem de s = y, é mesmo 
preciso alguma pista adicional.) Nesse caso, a equação para /i fica 

// -4y 2 + 4x 2 4x 2 (y 2 -x 2 ) 4x 2 (y 2 -x 2 ) x 

|Li y • 2xy 2 2y 3 — 4x 2 y — 2xy 3 — x (2y 3 — 4x 2 y) —4xy (y 2 — x 2 ) y 

x 2 x 

Ou seja, — = — j. Assim, /i(s) = s -1 = x/y é um fator integrante (verifiquem, e tenham 
cuidado com as coincidências). Multiplicando a equação original pelo fator integrante, 
obtemos 

2x 2 yy' + 2xy 2 — 4x 3 = 0, y(l) = l. 

Integrando, obtemos x 2 y 2 — x 4 = 0, que tem duas soluções diferenciáveis (verifiquem!): 
y = x e y = —x. Destas, só y = x é compatível com a condição inicial. 

(1.21) Exercício. Resolvam a equação 

6 xy 3 y / + 3y 4 = 5x 2 y 2 + 4x 3 yy', y( 2) = 2, 


sabendo que tem um fator integrante da forma /i(xy). 


(1.22) Equações lineares de primeira ordem. Uma equação diferencial diz-se linear se 
corresponde a um problema linear (homogéneo ou não homogéneo). Em particular, no 
caso homogéneo, o espaço das soluções é um espaço vetorial. 

Assim, a equação (nesta secção vamos usar t como variável independente) 


y' + ty = 0 
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Episódio 1: Equações diferenciais ordinárias 


é linear, pois a transformação T (do espaço C 1 das funções com derivada contínua para o 
espaço C° das funções contínuas) dada por 

Tíy) = / + t-y 


é linear (sim?) e a equação original corresponde ao sistema (homogéneo) T(y) = 0. 
Mais geralmente, uma equação linear de primeira ordem pode ser escrita na forma 

/(t) = a(t)-y(t) + b(t) 

e é sempre redutível a exata. Como habitualmente, abreviamo-la assim: 

—a{_t) ■ y — b(t)+ 1 -y' = 0. 


Usando a equação Ql. 18D com s = t, |^ = l,e|^ = 0, obtemos (recordem que estamos a 
usar t, e não x, como variável independente) 


ds 


õQ 

dt 


dP_ 

dy 


fJi 


0 — (—a(t)) 


= —a(0- 


1 • Q — 0 ■ P 

Nitidamente, o lado direito depende apenas de t, pelo que existe um fator integrante da 
forma pi(t). Resolvendo a equação, obtemos 


log|/i(t)| = 


a(t)dt. 


Assim, podemos escolher 


/ i(t) = e-f aCt)dt 


como fator integrante. Se tivermos uma condição inicial em t = t 0 , será preferível escolher 
uma primitiva mais específica: 

r >. — I a(s)ds 

/j(t) = e Jt o 

(reparem que tivemos que mudar a variável de integração—veem porquê?). 

Para simplificar as contas que se seguem, observamos desde já (a partir da equação 
acima para pd/í-O que pd(t) = —a(t)-/i(t). Multiplicando pi(t) na equação original, ficamos 
com 

/i(t) • / - pi(t) • a(t) • y = [i (t) • b(t). 

Mas isto é equivalente (verifiquem!) a 

( m (0 -y) 

cuja solução é 


Kt)-y = 


/i(t) • b(t) dt. 


Se tivermos uma condição inicial y(t 0 ) = y 0 , podemos ser mais específicos e dizer que 


r t 


■ y(0 - M(to) • yo = 


pi(r) • b(r) dr 
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(de novo, precisámos de mudar a variável de integração). Resolvendo em ordem a y(t), 
chegamos a (verifiquem!) 

r £ 

y(t) = e-^o a ^' )ds • y 0 + J e^- a( - s - )ds ■ b(r)dr . 

J t 0 

Esta é a fórmula de variação das constantes. Se olharem para uma equação diferencial e 
reconhecerem que é linear, podem usar esta fórmula diretamente, sem terem de perder 
tempo à procura de um fator integrante. 

(1.23) A primeira vista, a fórmula de variação das constantes parece não ter um padrão 
simples. Por isso, vale a pena explicar um pouco como traduzi-la em intuições sobre a 
passagem do tempo. Com efeito, podemos pensar que a função 

e / r ta <» ds 


codifica o efeito da passagem do tempo entre o instante r e o instante t. Chamemos-lhe por 
isso passagem*^'. Recordemos que a equação original era 

/(t) = a(t)-y(t) + b(t), j(t 0 ) = Jo- 


Intuitivamente, isto indica que a variação de y tem uma componente que depende do pró¬ 
prio y (portanto deveria envolver uma exponencial e a condição inicial—daí passagem í<_r 
envolver a exponencial) bem como uma contribuição adicional em cada instante r anterior 
a t, controlada por b(r). Isto é, a contribuição de y 0 vê o tempo passar de t 0 a t, enquanto 
que cada contribuição b(r) vê o tempo passar de r a t. Assim, o resultado deveria ser 


r t 


y(t) = passagem^ 0 • y 0 + 


passagem 


b(r)dr. 


Se compararem com a fórmula de variação das constantes, veem que é exatamente isso o 
que ela diz. 

Embora esta mnemónica possa parecer pouco fiável, a verdade é que há um conceito 
muito geral de exponencial que codifica esta ideia de passagem do tempo. Por exemplo, a 
exponencial habitual (na reta real ou no plano complexo) corresponde a passagem t ^ r = 
e t ~ r . As propriedades da exponencial correspondem a propriedades da passagem do tempo: 
dizer e t-t = 1 é dizer passagem 1 ^ = 1, e dizer e t-s • e s ~ r = e t-r é dizer passagem^ 5 • 
passagem s<_r = passagem t<_r (ou seja, é dizer que fazer o tempo passar de r a s e depois de 
s a t é equivalente a fazê-lo passar de r a t). 


(1.24) Exercício. Aqui estão mais algumas equações de primeira ordem para resolverem: 

(a) / —sint-y = sint, y(0) = -2; 

(b) 3y 2 y' + 2xy / + 2x + 2y = 0, y(0) = 2; 

(c) 6tyy' + 4y 2 y' + 2y 2 = ye t + 2e t y', y(4) = 0; 

(d) 2t 2 y + tyy' + ty' + 3t e f = t 3 + 2y 2 y' + 2 yy' + 6 y e c , y(0) = 2. 

Na última equação, precisam de um fator integrante do tipo /i(t — 2y). 
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(1.25) Mudança de variável. Por vezes, uma mudança de variável bem escolhida pode 
ajudar a encontrar soluções de uma equação diferencial. Por exemplo, na equação 

v + xy' = 2 y 

podemos substituir y(x) = u(x)x 2 . (Como esta mudança não é invertível em x = 0, sabe¬ 
mos desde logo que podemos não encontrar todas as soluções, ou que as que encontrarmos 
podem ter problemas, próximo de x = 0.) Como y' = u' x 2 + 2ux, substituindo na equação 
original obtemos 

x + x (tí 7 x 2 + 2 ux) = 2ux 2 , 
que depois de simplificada se reduz a 

x + u / x 3 = 0, ou u =- 

x A 

Somos assim levados a concluir que u = ^ + C, e que 

y = ux 2 = x + Cx 2 . 


(1.26) Exercício. Usando a equação original em Q1.25D , o que podem concluir sobre y(0)? 
Mostrem que é possível escolher a solução em x < 0 independentemente da solução em 
x > 0 e determinem todas as soluções da equação diferencial em R. 


(1.27) As substituições úteis nem sempre são óbvias, mas por vezes algum padrão na equa¬ 
ção pode sugerir opções a tentar. Por exemplo, em 

y' = 3y + 3 y 2/3 

é tentador livrarmo-nos da raiz cúbica usando y = u 3 (que, naturalmente, pode esconder 
problemas quando y = 0, pois embora nesse ponto a mudança de variável seja invertível— 
u = y 1 / 3 —, a inversa não é diferenciável). Substituindo y' = 3u 2 u' na equação original, 
obtemos (é preciso algum cuidado com os sinais na aritmética com expoentes fracionários) 

3u 2 u'= 3u 3 + 3u 2 , ou u'= u + 1. 


(1.28) Exercício. Encontrem um fator integrante adequado para u' = u + 1 e determinem 
todas as soluções dessa equação. Quais as soluções correspondentes para a equação original 
y' = 3y + 3y 2 / 3 ? Há mais soluções? 


(1.29) Exercício. Usem substituições adequadas para resolver as seguintes equações. 

(a) y / = 2y + 4y 3/4 , y(0) = l. 

(b) y / = 4y + 2y 1 / 4 , y(l) = 1. 

(c) cosyy 7 = e l + siny, y(l) = 0. 

(d) y' = te 2t y + 2y logy, y(2) = l. 

Basta identificarem uma solução, mas tentem perceber se pode haver mais. 
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(1.30) Respostas de exercícios. 


dl.3D basta exigir — ^A 1 + B 1 = — ^ A 2 + B 2 . 

d> (a)y(0=|f 3 + 2, tel. 

(b) y(x) = s m±,x>-V2. 

(c) yM = -fi^,x>-V2. 

Cd)/(t) = teR. 

(e) /(t) = ^ + - t ,set <0; f(t ) = § + j, se 
t > 0. 

(f) u(x) = C exp(— -) (com C V 0), se x < 0; 
u(x) = — exp(l — i), se x > 0 . 


<m> 

y(*) 

= 2 exp( 

f-x). 


y(*) 

= exp (2 

— e~ x ) — 1 . 

CÜ> 

y (0 = 

= tan(| 

+ t), para 

-1 - 

\/i + 

n < t < 

-1 + Vl + 71. 


(a) y(t) = 

C- 

1 

-31- 4 ' 


(b) o maior intervalo contendo a condição 
inicial e excluindo os pontos — 3 ± V9 + 2C. 

(c) y(t) = 0 . (f) y(t) = 0 . 

( TTTTÕl ) (a) y(t) = C(t + 3) (com C € R). 

(b) y(t) =- (com C e R, e no 

C — log|t + 3| 

maior intervalo excluindo zeros do 

denominador) ou y(t) = 0 . 

C 

(c) y(t) = - (com C V 0, e no maior 

intervalo excluindo t = ±l) ouy(t) = 0 . 


1 -x 4 
(e) y(x) =- 


-/(x 4 - l ) * 1 2 + 4x (x + 3) 

— 2 x 


(f) y(x) 

(g) yQO 


4 

C0S(7IX) — 1 + X 2 

x 2 — 2x + 3. 


( |1.17| ) (alguns dos fatores integrantes possíveis, 

seguidos da solução) 

, , _2 r , (e + 1) 2 

(a) x 2 ; y(x) = 


e x + 1 


(b) sinx ou y , y(x) =-. 

sinx 

1 

(c) y; y(x) = -. 


(d) y y(x): 


3 - e* + ^(e x -3) 2 -4x 


dl. 21 | ) um fator possível é x 2 y 2 ; y(x) = x. 


dl.24D (a) y(t) = —1 — exp(l — cos t). 

(b) solução implícita de y 3 + 2 xy + x 2 = 8 . 

(c) um fator possível é y; y(x) = 0 . 

(d) um fator possível é (t — 2 y) _1 ; 

y(t) = -1 + t 3 . 


dl.26D y(0) = 0. y(x) = x + C 1 x 2 se x < 0; 
y(x) = x + C 2 x 2 se x > 0. 


( fL28l ) u(t) = Cé - 1 (com C e R). 
y(t) = (Ce c — l) 3 . y(t) = 0 e ainda colagens de 
soluções dos dois tipos (vejam o exercício d2.26| ) 
para mais pistas). 


d!.13D (a) y(x) = \j 1 + x-x 2 . 

cosx — \/ (cosx ) 2 + 24 


(b) y(x) = 


(c) y(x) = 


cos x + \/ (cosx ) 2 — 4x 3 + 4 X 7 I 2 


(d) y(x) = x 


- v 3 /2 


2 x 


dR29l ) (a) y(t) = (3e ^ 2 — 2) 4 . 

(b) y(t)=(fe 3 t - 3 -I) 4/3 . 

(c) y(t) = arcsin(te t — e f ). 

(d) y(t) = exp(^e 2t - 2 e 2t ). 

O exercício ( |2.27| ) dá algumas pistas sobre as 
outras soluções. 
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EPISÓDIO 2 


EXISTÊNCIA DE SOLUÇÕES DE EDOs 


Nada do que vimos até aqui nos dá indicações sobre a existência ou unicidade das soluções 
de uma equação diferencial. Como discutimos na introdução, em geral esperamos que, se 
dois carros começam a viagem lado a lado e as suas velocidades variam da mesma forma 
ao longo do tempo, então deveriam continuar lado a lado ao longo de toda a viagem. 

Começamos por falar um pouco dos retratos de fase. Depois vemos que as equações 
diferenciais podem ser convertidas em equações integrais, para as quais uma discussão algo 
técnica nos dá um critério simples com que assegurar existência e unicidade das soluções. A 
primeira vista, o critério parece restringir-nos a intervalos mais limitados do que queremos, 
mas vamos ver que por vezes é possível prolongar as soluções (sempre de forma única) a 
intervalos maiores. 


(2.1) Retratos de fase. Para percebermos melhor o que acontece, vamos recorrer aos re¬ 
tratos de fase. A ideia fica mais clara se começarmos com duas variáveis dependentes. Ima¬ 
ginem que queremos perceber as soluções de 

xí(t) =-x 2 (t), xí,(t) = X^t). 

Podemos escrever isto de forma mais compacta pondo x(t) = (x^t), x 2 (t)) e dizendo que 
a equação é 

x / (t) = F(x(t)), 

onde F é um campo vetorial dependente de x(t) (F poderia depender também de t, mas 
neste exemplo isso não se passa). Nestes termos, o que queremos é encontrar um caminho 
x parametrizado por t, cujo vetor velocidade seja F(x(t)). Equações com a mesma F mas 
condições iniciais diferentes, corresponderão a caminhos diferentes no retrato de fase (que é 
essencialmente uma representação do campo vetorial F). Em geral, não convém representar 
as velocidades na mesma escala que as posições (de outra forma, a figura tende a ficar 
confusa). 


(2.2) Quando temos só uma variável dependente, faz sentido incluir também a variável in¬ 
dependente no retrato de fase. Não só é mais fácil percebermos o gráfico das soluções dessa 
forma (por exemplo, se t for a variável independente e y a variável dependente, podemos 
ler diretamente os gráficos das soluções y(t)), mas também porque F pode perfeitamente 
depender da variável independente. 

Por exemplo, voltemos à equação do exercício Q1.9D : 

/ = (3 + t)y 2 . 


21/99 



João Pedro Boavida, Episódios de Equações Diferenciais, 2015 



Se pusermos v(t) = (t,y(t)) € M 2 , temos v^t) = ( t',y') = (1 ,y'). Assim, a equação pode 
ser reescrita como 

V(t) = F(v(t)), onde F(t,y) = (1, (3 + t)y 2 ). 

Isto dá-nos um novo retrato de fase, que revela as soluções correspondentes. 


(2.3) O teorema de Picard-Lindelõf dá-nos condições concretas sob as quais uma equa¬ 
ção diferencial (com condição inicial) tem solução única. Para fixar ideias, vamos trabalhar 
com vetores em M n . (Nas figuras vamos usar n = 1, apenas por uma questão de clareza.) 
Supondo que F(t, x) está definida nos pontos que usarmos, o nosso objetivo é decidir se a 
equação 

x / (t) = F(t,x(t)), x(t 0 ) = x 0 



como definição de uma nova função, à qual chamaremos (à falta de melhor nome) < f>(x). 
Nessa linguagem, x é uma solução da equação original precisamente quando <f>(x) = x for 
satisfeita. 


Ora, isto sugere a seguinte estratégia: começamos com uma função qualquer (por 
exemplo, X( 0 j(t) = 0), e definimos recursivamente X( n+1 ) = <f>(x(„)). Se tudo correr bem, o 
limite destas iteradas de Picard é uma solução da equação original. 
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(2.4) Como isto é algo confuso, é melhor vermos um exemplo. Considerem a equação (não 
vamos usar notação vetorial, porque estamos em R 1 ) 


x / (t) = x(t), x( 0 ) = l. 


A função F correspondente é dada por F(t,x) = x, enquanto que <f> é definida por (notem 
que, conforme o contexto, usamos x para denotar um real, ou para denotar uma função de 
valores reais) 


<f>(x)(t) = 1 + x(s)ds. 


Jo 

(Verifiquem que tudo isso faz sentido, e que nesta fórmula x é sempre uma função.) 
Assim, se começarmos com x^(t) = 0, temos 


*(i)(t) = l+ *(o) 00 ds = l + 


Ods = 1; 


Jo 


Jo 



(2.5) Exercício. Determinem ainda X( 4 ) e X( 5 ). Qual acham que é a fórmula para x^(t)? 
Conseguem prová-la? Qual a função limite? 

(2.6) Exercício. Suponham que tínhamos começado com x^(t) = 1. Quais seriam as fun¬ 
ções X(q), X( 2 ), V( 3 ) correspondentes? E se começássemos com *a»(0 = t? 

(2.7) Exercício. E se a condição inicial fosse x(2) = —3, quem seria <f>(x)(t)? Começando 
com X( 0 )(t) = 0 , quem seriam as iteradas x^p X( 2 ), correspondentes? 

(2.8) Exercício. Considerem a equação 

x\t) = 4 - x(t), x(l) = 3. 

Quem são F(t,x) e 4>(x)(t)? Partindo de x^(t) = 1, determinem X( 3 ). 

(2.9) Exercício. Considerem a equação 

y'e x = l+y, y( 0 ) = e-l. 

do exercício ( |1.7P . Quem são F(x,y) e <f>(y)(x)? Partindo de y(o)( x ) = determinem y^ 2 )- 

(2.10) Exercício. Qual a equação diferencial e respetiva condição inicial associadas a 


4>(u)(t) = —2 + s + iz(s)ds? 


J2 


Se ri ( 0 ) (t) = 1 e u (n+1) = 4>(u (n) ), qual o valor de lim,,^^ u (n) (t)? 


23/99 








João Pedro Boavida, Episódios de Equações Diferenciais, 2015 


(2.11) Constantes de Lipschitz e operadores contrativos. É claro que nada garante que 
as coisas corram bem. Para isso, precisamos de uma discussão algo técnica. Se a quiserem 
saltar, continuem no Q2.21P . 

Vamos então assumir que F(t,x) está definida num aberto contendo (t 0 ,x 0 ). Além disso, 
vamos assumir que F é contínua, e localmente Lipschitz em relação a x. Recordem que a 
última condição significa que, para cada ponto no domínio de F, existe uma constante C e 
uma vizinhança do ponto tais que 

l|F(t,Yi) — F(t,y 2 )|| < C Ik — y 2 || 


para quaisquer dois vetores (t,yi) e (t,y 2 ) na vizinhança. (Uma tal constante é chamada 
uma constante de Lipschitz .) Notem que esta condição é satisfeita se F for C 1 (e este é o caso 
mais frequente). 


(2.12) Exercício *. Mostrem que se F é diferenciável, então é localmente Lipschitz se e só 
se as suas derivadas são localmente limitadas (isto é, se existe uma vizinhança de cada t, x 
onde as derivadas são limitadas). 


(2.13) Assumindo por um momento que não há problemas com os domínios das funções, 
dadas duas funções yXO e y 2 (t), temos então 

r t 


< KyiXO = x 0 + 


F0byi0ü)ds, 


r f 


$(5^X0 = x 0 + 


FO,y 2 Oü)ds, 


r f 


■KyiXO - $(y2Xt) = 


FObyiO)) -F(s,y 2 (s))ds, 


donde se segue que 


11^1X0-^2X011 = 


r t 


F(s,yi(0) -F(s,y 2 (s))ds 


r t 


< 


HFOqyXs)) -F(s,y 2 (s))||ds 


r t 


< C 


IIO,yiO)) - 0,y 2 (0)|| ds<CD\t- t 0 |, 


se chamarmos D à distância máxima entre y x (s) e y 2 (s) (para os valores de s usados). Além 
disso, 


IkCyiXO-xoll = 


r f 


FObyXO)^ 


r t 


l|F(s,yi(s))|| ds <M\t- t 0 |, 


se chamarmos M ao máximo de ||F(s,y 1 (s))|| (para os valores de s usados). 

Voltemos então ao t 0 e ao x 0 que definem a condição inicial. Como F está definida e é 
Lipschitz numa vizinhança de (t 0 ,x 0 ), podemos escolher D > 0 e e > 0 tais que 

\t — t 0 | < e e ||x —x 0 ||<D => (t,x) e domínio de F. 


Sendo M o supremo de ||F(t,x)|| (para \t — t 0 | < e e ||x — x 0 || < D) e C a constante de 
Lipschitz de F nessa vizinhança, encolhemos e (se necessário) até que Ce < 1 e Me < D 
(veem como?). 
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(2.14) Consideramos então o espaço vetorial V das funções ] t 0 — e, t 0 + £■[ —»R n contínuas 
e limitadas, e nele definimos a seguinte norma: 

llyllv = su p lly(t)ll- 

t:|t-t 0 |<£ 

Podemos assim falar em convergência de sucessões (y n —» y significa ||y n — y|| v —* 0) e de 
sucessões de Cauchy (y n é uma sucessão de Cauchy se ||y m — y n ||y —^ 0 quando m, n —* oo). 

(2.15) Exercício. Mostrem que V é realmente um espaço vetorial. 

(2.16) Exercício. Mostrem ainda as três propriedades relevantes da norma: 

(a) se ||y|| v = 0 , então y(t) = 0 ; 

(b) llyi +y2Üv ^ llyillv + Ily2llv; 

(c) seAel, então || Ay|| v = |A| ||y|| v . 

(2.17) Exercício*. Mostrem que ||y n — y|| v —» 0 (aquilo que acima escrevemos como y n —* y) 
significa que y n converge uniformemente para y. 

(2.18) Exercício*. Mostrem que as sucessões de Cauchy em V convergem. Para isso, consi¬ 
deramos uma sucessão de Cauchy y n . Isto significa que 

||y m - yjy —> 0 quando m, n —» oo. 

Queremos construir uma função yeVe mostrar que y n converge uniformemente para y. 

(a) Para cada t e ]t 0 — e, t 0 + e[, mostrem que y n (t) é uma sucessão de Cauchy em M. 

(b) Construam uma função y : ] t 0 — e, t 0 + £[ — 1 * K que seja limite pontual de y n . 

(c) Mostrem que y n —*> y uniformemente. 

(d) Mostrem que y e V. 


(2.19) Consideramos então o subconjunto Y ç V das funções y cujos valores satisfazem 
||y(t) — x 0 1| < D. Tudo o que fizemos acima assegura-nos que se y está em Y então também 
<f>(y) está em Y. (Verifiquem!!) 

Dadas duas funções y 1; y 2 e Y, vimos há pouco que 


r [ 


ll*(yi)(0-*(y2X0ll<c 


|yiCs) — y 2 C5)|| ds. 


J tr) 


ou seja (verifiquem!), 

Il*(y 1 )-*(y 2 )llv - Ce Hyj -y 2 |lv, 

onde Ce < 1. (Quando esta desigualdade é válida com uma constante menor que 1, dizemos 
que <f> é um operador contrativo .) 

Mas nós tínhamos dito que x era uma solução da equação diferencial precisamente se 
4>(x) = x. Suponham então que y x e y 2 são duas soluções. Nesse caso, temos 


Il^(yi)-^(y 2 )llv = llyi -Y 2 llv < Ce llyi -y 2 || v . 

Yi Y2 
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Mas como Ce < 1, isto só é possível se H)^ — y 2 || v = 0, ou seja, se 3 ^ = y 2 . Isto mostra que, 
a existir solução da equação diferencial, essa solução está unicamente definida (pelo menos 
no intervalo [t - t 0 1 < e). 

Mas será que existe alguma solução? Fazemos o seguinte: escolhemos um qualquer 
ponto de partida X( 0 ) e construímos a sucessão definida por X( n+1 ) = <f>(x( n )). Vamos verificar 
que é uma sucessão de Cauchy. 

(2.20) Exercício. Para abreviar um pouco, pomos R = Cs. 

(a) Mostrem que ||x (n+1) - x (rt+2) || v < R||x (n) - x (n+1) || v . 

(b) Usem a desigualdade triangular para mostrar que, para k inteiro positivo, 

ll x (n) x (n+fc)IIV ^ ll x (n) _x (n+l)llv + ll x (n+l) _ x (n+2)IIV d-^ ll x (n+íc-l) _x (n+fc)llv- 

(c) Mostrem que, para k inteiro positivo, 

R n -R n+k R n 

H X ( n ) x (n+fc)llv ^ 1 _ R ■ H x ( 0 ) — x Cl)llv < -j—j^ ' H X (0) — X (l)Hv- 

(d) Mostrem que X( n ) é uma sucessão de Cauchy. 

(e) Mostrem que o limite de Xç n j também está em Y. 


(2.21) Prolongamento de soluções. Resumindo... Uma equação diferencial 

x(t) = F(t,x(t)), x(t 0 ) = x 0 , 

onde F é localmente Lipschitz, tem uma única solução definida nalguma vizinhança de t 0 . 
Vamos mostrar que existe um intervalo máximo (contendo t 0 ) no qual a solução é única. 


Lembram-se que no final de Q2. 13[ ) encolhemos o intervalo de definição de modo a ter 
Ce < 1? De facto, em subintervalos de comprimento e ou menor, o teorema de Picard- 
Lindelõf aplica-se (veem porquê?). Fazemos então o seguinte: dividimos o intervalo ori¬ 
ginal em subintervalos de comprimento e. No subintervalo que contém a condição inicial, 
encontramos a (única) solução, e usamo-la para decidir a condição inicial no subintervalo à 
esquerda e no à direita. Dessa forma, vamos prolongando a solução para a esquerda e para 
a direita, sempre com solução única (porquê?). Só somos forçados a parar quando saímos 
do intervalo original, ou do conjunto onde a constante de Lipschitz C é válida. 

Assim, vemos que há apenas duas formas de perder unicidade de solução: ou ao longo 
de t chegamos a um ponto onde F não é localmente Lipschitz, ou x “foge” de qualquer 
conjunto limitado (isto é, tem uma assintota em alguma das variáveis). 


(2.22) Exercício. Considerem a equação y' = —e~ y . 

(a) Encontrem conjuntos nos quais se aplique o teorema de Picard-Lindelõf. 

(b) Qual a solução geral da equação? Contradiz alguma das nossas conclusões sobre 
possibilidade de prolongamento de soluções? 
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(2.23) Exercício. Considerem a equação diferencial 


(a) Qual a função F(t, y) correspondente a esta equação? 

(b) Encontrem a solução desta equação cuja condição inicial é y(l) = 1. Mostrem que 
essa solução está definida para todo otei. Qual o seu valor em t = 0? 

(c) Encontrem a solução desta equação cuja condição inicial é y(2) = 1. Mostrem que 
essa solução está definida para todo otei. Qual o seu valor em t = 1? E em t = 0? 

(d) Encontrem a solução desta equação cuja condição inicial é y(3) = 1. Mostrem que 
essa solução está definida para todo oteM. Qual o seu valor em t = 3? E em t = 0? 

(e) Encontrem todas as outras soluções desta equação. Quais delas têm y(0) = 0? Isto 
contradiz o teorema de Picard-Lindelõf? 


(2.24) Exercício. Considerem a equação 

t-y' = 2y. 

(a) Qual a solução geral desta equação quando t > 0? E quando t < 0? 

(b) Se y for uma solução da equação definida em t = 0, qual é o valor de y (0)? Mostrem 
que coincide com o limite quando t —» 0 das soluções da alínea anterior. 
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(c) Usem as alíneas anteriores para construir todas as soluções da equação definidas em 
E. Em que pontos do domínio se aplica o teorema de Picard-Lindelõf? 


(2.25) Exercício. Encontrem todas as soluções de 

y' = 2 y 1//2 se y > 0 e y' = 0 caso contrário. 

(2.26) Exercício. Considerem de novo a equação 

y' = 3y + 3 y 2/3 


do exemplo fll.27[ ) e dos exercícios que se seguem. Que valores de t e y satisfazem as 
hipóteses do teorema de Picard-Lindelõf? Determinem todas as soluções da equação. 


(2.27) Exercício. Considerem de novo as soluções que encontraram para as equações do 
exercício fll.29D - 

(a) y' = 2y + 4y 3/4 , y(0) = l. 

(b) / = 4y + 2y 1 / 4 , y(l) = 1. 

(c) (cosy)y / = e t + siny, y(l) = 0. 

(d) y' = te 2t y + 2ylogy, y(2) = l. 

Qual o maior intervalo onde as soluções são únicas? Quais as soluções fora desse intervalo? 

(2.28) Exercício. Considerem a equação 

y' = l + |y| 1/2 . 

(a) Mostrem que todas as soluções são funções crescentes. 

(b) Encontrem expressões para as soluções, quando y > 0. (Pode ser útil usar y = u 2 .) 

(c) Encontrem expressões para as soluções, quando y < 0. (Pode ser útil usar y = — u 2 .) 

(d) Encontrem todas as soluções, e mostrem que, apesar desta equação não satisfazer 
as hipóteses do teorema de Picard-Lindelõf, a conclusão é verdadeira. 


(2.29) Exercício. De facto, é possível mostrar a existência e unicidade de soluções num 
aberto O no qual 


(cyi),(t,y2) e 


l|F(t,yi) — F(t,y 2 )|| < c ||yi — y 2 ll- 


Para tal, consideramos o espaço V das funções definidas num intervalo ]a, b[ contendo t 0 . 
(a) Se a < 0 < b e a < t < b, mostrem que 


r t 


Ce c| s | -c | t| ds 


0 


< 1 — e 


—C max{|a|,|ò|} 


(b) Se a < t 0 < b e a < t < b, mostrem que 

r t 

ÇgC|s— 1 0 | — C|t—t 0 l d5 


< l_ e -C|b-a|_ 
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(c) Verifiquem que 

||y||;= sup (V c|t ~ to l ||y(t)lQ- 

a < t < b v y 

é uma norma em V. Tenham cuidado para não confundir ||y||y com ||y(t)||. 

(d) Mostrem que 

H^yr) - *(y 2 )||v < (i - e- c|i, - Q| ) ||yi -y 2 ||v• 

(e) O que podem concluir sobre existência e unicidade de soluções de x'(t) = F(t,x(t)) 
em íl? Que detalhes é preciso ter em conta sobre a relação entre ]a, b[ e íl? 


(2.30) Comparação de soluções. Digamos que y(t) é solução de y'(t) = F(t,y{t)) com 
condição inicial y(t 0 ) = Jo e fi ue z(0 é solução de z\t) = G(t,z(t)) com condição inicial 
z(t 0 ) = z 0 . Se y 0 — z o e F(t,y) — G(t,z ) sempre que y < z (notem que aqui y e z são 
valores, não funções), conseguimos mostrar que y(t) < z(t) sempre que t > t 0 ? 

Se pudermos usar o teorema de Picard-Lindelõf, a resposta é sim. 

(2.31) Exercício. Começamos por ver o que se passa quando y 0 < z 0 . 

(a) Mostrem que se F(t,y(t)) > G(t,z(t)) para algum t > t 0 , então z(t) — y(t) não 
pode ser uma função crescente. 

(b) Mostrem que nesse caso existe algum t 1 > t 0 com z^tq) — y / (ti) < 0, mas z(t) — 
y(t) > 0 no intervalo ]t 0 , tq[. 

(c) Mostrem que isso contradiria as condições sobre F e G. 

(2.32) Exercício. Consideramos agora o que se passa se y 0 = z 0 , e assumimos que G satisfaz 
as condições do teorema de Picard-Lindelõf. 

(a) Mostrem que se z(q) < y(Q) para algum t 1 > t 0 , é possível escolher uma solução 
w(t) de w/(t) = G(t, w(t)) (notem que é a mesma equação satisfeita por z(t)) com w(t 0 ) > 
z (to) = z o = Jo e z(ti) < w(q) < y(t a ). 

(b) Mostrem que isso contradiria as conclusões do exercício anterior. 

(2.33) Exercício. Considerem as equações w 1 = —w, y' = y siny, e z' = z, com condição 
inicial w(0) = y(0) = z(0) > 0. Mostrem que se t > 0, então w(t) < y(t) < z(t). 

(2.34) Exercício. Mostrem que as soluções de y' = e~ y2 e as de z' = ze~ z2 estão definidas 
em toda a reta real. 
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(2.35) Respostas de exercícios. 


! + t + f + f ; 1 + t +f + ¥ + ã 


V il 
k! ' 

0 <k<n 


t —> e . 


( |231 ) 


as anteriores x™, x. 


( 2 )’ a ( 3 )> 


( 4 )- 

.4 


i + ^l + t + T^ + f + T + i 


( f2?7l ) 4>(x)(t) = —3 + x(s)ds. 

-3; 3-3t; -3 + 3t- |t 2 . 

( pT8l ) F(t,x) = 4 — x. 

4>(x)(t) = 3 + 4 — x(s)ds. 

x ( 3 ) (t)=l + |t- 2 t 2 +±t 3 . 

( |2.9D F(x,y) = e~ x + ye~ x . 

Hy)M = e-l + fj e“ s + y(s) e" s ds. 
y ( 2 ) (x) = 2 e — 2 e _ * — e 1_x + e~ 2x . 


( | 2 JÕ 1 ) u' = t + u; ií(2) = - 2 . 
u(t) = e t_2 — t — 1 . 

d232t (b)y(t) = log(C-t). 

(a)E(t,y) = y 1 / 3 . 

(e) y(t) = (|t + C ) 3/2 se t > -\C e y(t) = 0 
caso contrário; ou y(t) = 0 em todo o R. 


Q2.24D (a) y(t) = Ct 2 (em ambos os casos). 

(b) y(0) = 0. 

(c) escolher uma constante para t < 0 e outra 
(ou a mesma) para t > 0; aplica-se em t / 0. 

( [2.25P y(t) = (t — C) 2 se t > C e y(t) = 0 caso 
contrário. y(t) = C, com C < 0. 

( 12361 ) y ^ 0. y(t) = (Ce f - l) 3 ; y(t) = 0; 
y(t) = —1. As soluções do primeiro tipo 
cruzam-se com as do segundo em t = — log C ; 
nesses pontos também podemos trocar de um 
tipo para outro. 

( |2.27P (a) ]21og(2/3),+oo[. É possível colar 
com y(t) = 0 à esquerda de t = 21og(2/3), e 
ainda mais à esquerda também é possível colar 
com outras soluções da forma 
y(t) = (Ce t/2 - 2) 4 . 

(b) ]1 — |log3,-(-oo[. É possível fazer colagens 
como em (a). 

(c) ]— 00 ,a[, onde (a — l)e a = 1. (d) R. 

( 12.28D (b) y(t) = u(t) 2 , onde u(t) é solução 
implícita de 2u — 2 log(l + u) = t — C e t > C. 
(c) y(t) = —u(t) 2 , onde u(t) é solução implícita 
de 2u — 21og(l +u) = C — t e t < C. 
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EPISÓDIO 3 


EDOs EM M n E EXPONENCIAIS DE MATRIZES 


Nos episódios anteriores falámos principalmente sobre equações diferenciais de primeira 
ordem, com uma única incógnita. Mas se tivermos um sistema como 


*í = 

x' 2 = -42x 1 + 13x 2 , 


podemos escrevê-lo na forma matricial 


*1 

/ 

0 

1 " 



*2 


-42 

13 


*2 



x' A x 

Considerem a mudança de variável x = Su, com 


S = 


1 1 

7 6 


Nesse caso, a equação x' = Ax fica Su 7 = ASu. Se escrevermos J = S~ 1 AS, temos SJ = AS, 
ou seja, a equação fica Su 7 = SJu ou (multiplicando à esquerda por S -1 , de ambos os lados, 
e simplificando) u 7 = Ju. Chegamos portanto a 


u 


/ 


’ 7 

0 " 


1 


6 

u, 

ou j 

0 


1 


= 7u 1; 


= 6ií 2 , 


As duas funções u 1 (t) e iz 2 (t) ficam completamente independentes, com equações que já 
sabemos resolver. Isto acontece quando a matriz J é diagonal (e portanto a matriz original 
A é diagonalizável). Quando as matrizes não são diagonalizáveis, usamos a forma canónica 
de Jordan. 

Neste episódio, vamos considerar principalmente equações da forma 


Y'(t)=A-Y(t) + B(t), 


onde Y(t) e £(t) são matrizes nxm eAé uma matriz n x n. Se chamarmos yi(t),... ,y m (t) 
às colunas de P(t) e b 1 (t),... ,b m (t) às de B(t), essa equação é equivalente às m equações 

yi(0 =A-yi(t) + b 1 (t), ..., y m (t) =A-y m (t) + b m (t). 

Desta forma, a notação matricial permite-nos representar várias equações simultaneamente. 
Omitiremos por vezes a variável t em y(t) ou Y(t), escrevendo apenas y ou Y. 
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(3.1) Exponencial de matrizes. Vamos chamar exponencial de A (ou de At) à solução M de 

M'(t) =A-M(t), M(0) = ld. 

(O teorema de Picard-Lindelõf garante que uma tal solução existe e é única. Conseguem 
verificar os detalhes?) Então a derivada de 

Y(t) = M(t)-Y 0 é Y / (t) = M / (t)-Y 0 =A-M(t)-Y 0 =A-Y(t), 

mostrando que essa Y(t) é a única solução—de novo por Picard-Lindelõf—de 

Y / (t)=A-Y(t), Y(0) = Y 0 . 

Sabendo isso, podemos verificar que (considerando s constante) 

M(t)-M(s) = M(t+s), 

pois (verifiquem!) ambas as expressões são (a única) solução de 

Y'(t) =A- Y(t), Y(0) = M(s). 

Por essa razão, podemos concluir que 

M(t)-M(-t) = M(0) = Id. 

Ou seja, M(t) é uma matriz invertível e a sua inversa é M(—t). Da mesma forma, podemos 
ver que 

M(t)-A-A-M(t) = 0, 

pois (verifiquem também aqui) ambas as expressões são (a única) solução de 

Y / (t)=A-Y(t), Y(0) = 0. 

Concluímos portanto que 

M(t)-A = A-M(t), 

isto é, a exponencial de A comuta com A. A um nível mais profundo, a matriz M(t) descreve 
o efeito da passagem do tempo. O facto de M(t) comutar com A está relacionado com 
M(t) ■ M(s) = M(t + s) (que poderia traduzir-se por “a passagem de tempo s seguida da 
passagem de tempo t é equivalente à passagem de tempo t +s”). Ambos são cruciais para 
que uma abordagem simplificada às equações em M n seja possível. (Daí focarmo-nos no 
caso em que A é constante.) 

Resumindo tudo isto: usando a notação e At em vez de M(t), acabámos de mostrar que 
e At • Y 0 é a única solução de Y'(t) = A - Y(t) com condição inicial Y(0) = Y 0 . Mostrámos 
(confirmem!) também as identidades 

gAR+s) = e At . g* (e At )-i = e~ At , e e At -A = A-e At = (e At )'. 

(3.2) Exercício. Conseguem mostrar que seAeB comutam, então 

e At ■ B = B ■ e At e e At ■ e Bc = e^ A+B) f ? 

O raciocínio é o mesmo que antes: mostrar que ambos os lados de cada equação são (a 
única) solução de uma mesma equação diferencial. 
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(3.3) A fórmula de variação das constantes para matrizes n x m pode ser obtida com um 
raciocínio semelhante ao que fizemos em dl.22[ ). Um “fator integrante” (se nos permitirmos 
usar tal designação neste caso) para a equação 

é e~ At multiplicado à esquerda, pois a equação torna-se 

e~ At YXt ) - e~ At ■A - Y(t) = e~ At -B(t), 
que é equivalente (verifiquem!) a 

(e~ At -Y(t)y = e~ At 

Se tivermos uma condição inicial em t = t 0 , podemos integrar e obter 

r t 


-At 


no - e- At ° • nt 0 ) = 


e Ar -B(r)dr. 


Tal como no caso n = 1, resolvemos em ordem a Y(t) e chegamos a 

r t 

e A(t-r) . B( r )dr. 


y(t) = e A b- t °)-y(t 0 ) + 

E tal como antes, podemos interpretar isto como 


to 


r t 


Y(t) = passagem *^ 0 • U(to) + 


passagem 1 ^' -B(r)dr. 


J t 0 


(3.4) Um exemplo concreto permite-nos ilustrar tudo isto. A solução de 


y / (0 = 


-8 1 
-4 -12 


y(0 + 


i + t 
2 - 2 1 


y(^oj = 

tn 


é dada por 


y(t) = e At ( 


e At °y 0 


b(t) 
r t 

i 

to 


0 

1 

y 0 


e“ Ar b(r)dr). 


Em (3.20), vamos ver que 


e At = e~ 10t 


1 + 2 t t 

—4t 1 — 2 t 


Substituindo na função integranda, temos 
e _Ar b(r) = e 10r 


CN 

1 

T—1 

1 _ 

—r 


1 + r ‘ 

- .IOr 

1 

CO 

1 

T —1 


' e 10r (1 — 3r) ' 

4 r 

1 + 2 r 


2-2 r 

— c 

2 + 6 r 


e 10r (2 + 6 r) 


Assim, o integral fica 
r t 


f „ior 


e A ' b(r) dr = 


o 


o . 


e 10r (1 — 3r) 
e 10r (2 + 6 r) 


dr = 


J 0 £ e 10r (1 — 3r) dr 
J 0 £ e 10r (2 + 6 r)dr 


(3.5) Exercício. Terminem as contas e obtenham a expressão simplificada (sem integrais) 
para y(t). 
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(3.6) Cálculo de e At . Há várias estratégias para calcular e At . Começamos por observar que 
se Y(t) é uma solução de Y / (t) = A - Y(t), então Y(t) = e At ■ Y(0) permite-nos concluir que 
Y(t) é invertível exatamente quando Y(0) o for (e portanto, quando qualquer outra Y(t 0 ) 
o for). 

Assim, suponhamos que A é uma matriz n x n e que temos n soluções yi(t),... ,y n (t) de 

y / (t)=A-y(t). 

Então a matriz Y(t) com colunas yi(t),...,y n (t) é uma solução de 

Y / (t)=A-Y(t), 

logo é da forma Y(t) = e A ^ t ~ t( " ■ Y(t 0 ). Se os vetores y^t),... ,y n (t), para algum t fixo, 
forem linearmente independentes, a matriz Y(t) será invertível. E vimos há pouco que 
nesse caso a matriz Y(t 0 ) é invertível. Somos assim levados a concluir que 

gA(t-to) = Y( t ). y(t 0 )-i. 

Esta construção de e At corresponde à matriz wronskiana. Se Z(t) for uma solução de 

Z / (t)=A-Z(t), 


podemos escrever 

Z(t ) = e A(t - to) • Z(t 0 ) = Y(t) • Y(to )" 1 ■ Z(t 0 ). 


(3.7) Exponencial de matriz diagonal. Voltemos então a um exemplo do início deste epi¬ 
sódio. Queremos encontrar todas as soluções do sistema 

y = A ■ y, onde A = 

Neste caso, é fácil ver que v x = [ J ] gera os vetores próprios com À = 7 e que v 2 = [ ° ] 
gera os vetores próprios com À = 6. 

Se v é um vetor próprio de A com valor próprio Ã, então y = e At v é uma solução de 
y 7 = A ■ y. Com efeito, temos 

/= ( e At - v y = (e At ) / • v = e At • v _Áy^ = e At ■ Av = A-e Xt v = A-y. 

Av 

No caso da nossa equação, temos então duas soluções (que são linearmente indepen¬ 
dentes pois os próprios vetores v x e v 2 o são): 

y 1 (t) = e 7 t v 1 = e y 2 (t) = e 6 t v 2 = g 6 t • 

Com estes dois vetores, podemos construir a matriz 

nt)=[yi(t) y 2 (0 ] = 0 e 6t • 

Como neste caso Y(0) = Id, vemos que e At = Y(t) • Y(O )” 1 = Y(t). Ou seja, a matriz 
exponencial associada a uma matriz diagonal é também ela uma matriz diagonal: cada 
valor próprio Ã é substituído por e ? t . 
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(3.8) Exercício. Usem a fórmula da variação das constantes em R 2 para resolver 


y(0 = 


7 0 
0 6 


y (0 + 



y(0) = 0. 


Comparem com o que obteriam escrevendo como sistemas de equações a uma variável e 
resolvendo cada equação usando Q1.22D . 


(3.9) Exercício. Vamos provar que se A ■ S = S ■ B, então e At ■ S = S ■ e Bt . Vamos usar 
um raciocínio semelhante ao que fizemos no exercício ( |3.2D e imediatamente antes. Tanto 
e At ■ S — S ■ e Bt como a matriz 0 são (a única) solução de Z'(t) = A ■ Z(t) com Z(0) = 0. 
Conseguem preencher os detalhes que faltam? 


(3.10) Exercício. Diagonalizem as seguintes matrizes e determinem as suas exponenciais. 


(a) 


11 3 
-6 2 



(c) 


-9 10 
-8 9 


(d) 


-9 4 
-5 0 


(Notem que, para ser mais fácil inverter a matriz dos vetores próprios, pode ser boa ideia 
escolhê-los de modo a que a matriz tenha determinante 1.) 


(3.11) Exercício. Resolvam os seguintes sistemas de equações diferenciais com t como va¬ 
riável independente: 

(a) x[ = llx 1 + 3 x 2 , x ' 2 = -6x 1 + 2x 2 , x 1 (0) = 0, x 2 (0) = 0; 

(b) Xj = —Xj - 12x 2 + 31, X 2 = 4x x + 13x 2 - 2t, Xj(0) = 0, x 2 (0) = 0; 

(c) Xj = -9x x + 10x 2 , x' 2 = - 8 x 1 + 9 x 2 , x 1 (0) = 0, x 2 (0) = 1; 

(d) x^ = -9x 1 + 4x 2 + t, x' 2 = —5x : + t, x : (0) = 5, x 2 (0) = 6. 


(3.12) Vejamos mais um exemplo. Se 


A = 


4 0 0 

0-3 0 

0 0 11 


então e l5 e 2 , e e 3 (os vetores da base canónica) são vetores próprios de A, pois A - e x = 4e 1; 
A - e 2 = —3e 2 , e A - e 3 = lle 3 . Por isso, 


e At ■ e x = e 4t e 1; 


e At • e 2 = e 3t e 2 , 


e At ■ e 3 


e m e 3 


e concluímos que 


e At = 


> 4t 0 0 
0 e~ 3t 0 
0 0 e Ut 
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(3.13) Exponencial de bloco de Jordan. Suponham agora que temos 


y = A ■ y, onde A = 

Para esta matriz A, temos apenas um vetor próprio independente: v 1 = [J], Pelo 
raciocínio que fizemos há pouco, sabemos que yi(t) = e^Vj = Vj é uma solução da equação. 
Como a solução também é dada pela exponencial, obtemos 

e At • v x = Vl . 

Temos A • v 2 = 0v 2 + v 1 . Queremos calcular y 2 ( t) = e At ■ v 2 . Mas 
y 2 (t) = (e At • v 2 y = e At • TEv, = e At • v 1 = Vl . 

Vl 

Integrando ambos os lados da equação entre 0 e t, obtemos 

y 2(0 = Y 2 (°) + tv n ou seja, ’ v 2 = v 2 + tvj. 


0 1 
0 0 


Se construirmos uma matriz com colunas yi( 0 ,y 2 (t)> 

H0= [ yi(0 y2(0 ] = 


1 - 

r 1 11 

J - 

i 

T — 1 

O 


temos, mais uma vez, Y(0) = Id. Portanto, e At = Y(t) • T(0 ) _1 = Y(t). 

É razoável perguntar o que acontece quando o valor próprio não é 0. Por exemplo, 
consideremos 

" 7 1 
0 7 


B = 


Mas nesse caso A = B — 7 lá tem valor próprio 0, por isso 

, logo e Bt = e 7t 

Consideremos agora o bloco de Jordan n x n com valor próprio A = 0 


e Bt ■ e ~ 7t = e Bt ■ e~ 71dt = e (B “ 7Id:) 1 = e At = 


1 t 
0 1 


1 t 
0 1 


N = 


0 10 0 0 
0 0 10 0 
0 0 0 1 0 
0 0 0 0 1 
0 0 0 0 0 


Se chamarmos e 1; ... e 5 aos vetores da base canónica, podemos ver que N-e k = e k _ 1 quando 
k > 1 e N ■ ei = 0. Como é que isso ajuda? 

Escrevendo y fc (t) = e Nt ■ e k , já vimos que y^t) = e 2 e y 2 (t) = e 2 + te^ Temos ainda 

yjt(0 = (e Nt • e k y = e Nt • N • e k = e Nt • (0e fc + e*^) = e Nt ■ e fc _j = y k - X {t). 


Assim, temos 

3^(0 = y2( t ) = e 2 + te 1 . 
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Primitivando e usando y 3 (0) = e 3 , chegamos a 

y 3 (0 = e 3 + te 2 + Çe 1 . 

Com o mesmo género de raciocínio, obtemos ainda 

y 4 (0 = e 4 + te 3 + je 2 + Ijer; 
YsCO = e 5 + te 4 + V e 3 + b e 2 + 4 ! e l- 


Construímos agora uma matriz a partir dos y fc (t): 


y(0=[yi(t) y 2 (0 y3(0 y4(0 ys(0 


1 

t 

t 

2 ! 

t 

3! 

t 

4! 

0 

1 

t 

t 2 

2 ! 

t 3 

3! 

0 

0 

1 

t 

£ 

2 ! 

0 

0 

0 

1 

t 

0 

0 

0 

0 

1 


Como 7(0) = Id, concluímos que e Nt = Y(t) • 7(0) 1 = Y(t). 
Se quisermos considerar o caso geral de 

A 1 0 0 0 

0 A 1 0 0 

J= 0 0 A 1 0 

0 0 0 A 1 

0 0 0 0 A 


podemos observar que J — Aid = IV. Portanto, 


£ Jt _ e (AId +N)t _ pAId t . p Nt _ p Xt 


= e 


e = e 


1 t 
0 1 
0 0 
0 0 
0 0 


2 ! 

t 

1 

0 

0 


3! 

t± 

2 ! 

t 

1 

0 


4! 

£ 

3! 

£ 

2 ! 

t 


Por exemplo, se 


então a sua exponencial é 


J = 


7 10 0 0 
0 7 10 0 
0 0 7 1 0 
0 0 0 7 1 
0 0 0 0 7 


e Jt = e 7t 


1 t 
0 1 
0 0 
0 0 
0 0 


£ 

2 

t 

1 

0 

0 


r 

3! 

£ 

2 

t 

1 

0 


4! 

£ 

3 2 

£ 

2 

t 
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(3.14) Outro caso relevante é o das matrizes que se podem decompor como 


C = 


' A 

0 ' 

0 

B 


(onde Ae B são matrizes quadradas e cada 0 é uma matriz nula com as dimensões necessá¬ 
rias para que C seja outra matriz quadrada). Vamos mostrar que 


e 


Ct 


e At 

0 

0 

e Bt 


Lembram-se da definição da exponencial de matrizes? Dissemos que e Ct é a (única) solução 
de M\t) = C ■ M(t) com M(0) = Id. Para provar que e Ct é realmente dada pela fórmula 
acima, basta verificar essas duas condições: 


e At 

0 

y_ 

■ (e At ) / 

0 


A - e At 

0 


" A 

0 " 


e At 

0 

0 

e Bt _ 

J ~ 

0 

(e Bt y ] 


0 

B ■ e Bt _ 


0 

B 


0 

e Bt _ 


' 0 a0 

1 

O 


’ Id 

1 

o 

0 

e S0 


1 

O 

Id J 


(Conseguem verificar que isto chega?) 

(3.15) Exercício. Determinem as exponenciais das seguintes matrizes. 



r o i i 


O 

O 

LO 


O 

O 

T— 1 

1 

(a) 

i 

O 0 

CO h- 

1_ 

; (b) 

0 4 1 

0 0 4 

; (c) 

0-1 1 

0 0-1 


(3.16) Exercício. Resolvam os seguintes sistemas de equações diferenciais com t como va¬ 
riável independente: 

(a) x\ = 3x x + x 2 , x' 2 = 3x 2 + 1, x 1 (0) = 0, x 2 (0) = 0; 

(b) xí = 5x 1; x 2 = 4x 2 + x 3 , X 2 = 4x 3 , x 1 (0) = l, x 2 (0) = 0, x 3 (0) = l; 

(c) x[ = -x 1 + t, x' 2 = -x 2 + x 3 , Xg = x 3 , x 1 ( 0 ) = 0 , x 2 ( 0 ) = 1 x 3 ( 0 ) = 2. 

(3.17) Resta a pergunta: e como fazemos quando a matriz não é diagonal nem um bloco 
de Jordan? Podemos ter a sorte de encontrar soluções suficientes para usar a matriz wrons- 
kiana: e At = Y(t) • Y(0) -1 . Ou podemos tentar encontrar a forma de Jordan, e aproveitar 
o que já sabemos sobre exponenciais de matrizes diagonais e de blocos de Jordan. Feliz¬ 
mente, na prática a questão é bastante menos técnica do que esta longa discussão poderia 
fazer suspeitar. Em todo o caso, vamos dedicar algum tempo a uns quantos exemplos. 


(3.18) A forma de Jordan de uma matriz. Antes de mais, o que é uma forma de Jordan de 
uma matriz? É o que usamos quando gostávamos de ter uma diagonalização, mas não há 
vetores próprios suficientes. Uma matriz de Jordan é uma matriz com blocos de Jordan ao 
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longo da diagonal, e zeros em todas as outras componentes. Assim, 


4 

1 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

2 

1 

0 

0 

0 

0 

2 

1 

0 

0 

0 

0 

2 


são matrizes de Jordan, enquanto que 

5 2 0 0 
-13 0 1 
0 0 0 1 
0 0 0 0 


7 

1 

0 

0 

0 

7 

0 

0 

0 

0 

-2 

0 

0 

0 

0 

3 


2 3 

0 


3 1 

0 

0 2 

0 

e 

0 3 

0 

0 0 

11 


0 -1 

7 


não são. 

Não vamos demonstrar todos os detalhes. O que nos interessa aqui não é a demons¬ 
tração; o que realmente queremos é uma estratégia que nos permita obter uma forma de 
Jordan. Vamos descrever por alto uma tal estratégia, e depois explicá-la com uns quantos 
exemplos. 

Temos então uma certa matriz A. A primeira coisa a fazer é determinar o polinómio 
característico de A. Cada valor próprio A ocorre com uma certa multiplicidade (algébrica) 
nesse polinómio. Por exemplo, no polinómio (A — 7) 4 , o valor A = 7 tem multiplicidade 4. 

Para cada valor próprio A, calculamos a matriz A — Aid. Um vetor próprio generalizado 
de A com valor próprio A é um vetor no núcleo de (A — Ald) fc , para algum k > 1. Dois 
pontos cruciais: quando aumentamos k nunca “perdemos” vetores próprios generalizados, 
e nunca é preciso usar k maior que a multiplicidade algébrica de A. 

Então a estratégia será a seguinte: começamos por determinar a nulidade de A — Aid. 
(Recordem que a nulidade de uma matriz é o número de variáveis livres—ou de colunas 
sem pivot—e é também a dimensão do núcleo. O núcleo de A — A Id consiste nos vetores 
próprios habituais.) Se houver vetores próprios suficientes, determinamo-los e temos o 
problema resolvido. Senão, determinamos a nulidade de (A — Aid) 2 , (A — Aid) 3 , (A — Aid) 4 , 
etc., até que coincida com a multiplicidade algébrica. Quando isso acontecer, significa que 
não vale a pena aumentar o expoente. 

Imaginem que (para uma matriz A) o valor próprio A = 7 tem multiplicidade (algébrica) 
4. Os blocos de Jordan para A = 7 podem ser de vários formatos diferentes: 


7 

0 

0 

0 

0 

7 

0 

0 

0 

0 

7 

0 

0 

0 

0 

7 


7 

1 

0 

0 

0 

7 

0 

0 

0 

0 

7 

0 

0 

0 

0 

7 


7 

1 

0 

0 

0 

7 

0 

0 

0 

0 

7 

1 

0 

0 

0 

7 


7 

1 

0 

0 

0 

7 

1 

0 

0 

0 

7 

0 

0 

0 

0 

7 


etc. 


Se A — 7 Id tiver nulidade 4, significa que conseguimos encontrar quatro vetores pró¬ 
prios independentes e portanto podemos diagonalizar A. Obtemos blocos como na primeira 
matriz mostrada. 

Se A — 7 Id tiver nulidade 3, temos que subir o expoente. Como a nulidade de (A — 7 Id) 2 
é 4, podemos parar. Escolhemos então qualquer vetor v 2 no núcleo de (A — 7Id) 2 . Em 
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seguida, escolhemos v x = (A — 7Id) • v 2 . (Tipicamente, v 1 ^ 0; se não for, temos de voltar 
atrás e escolher outro v 2 .) A consequência dessa escolha é que (A — 7Id) • v 2 = v 1; ou 
A ■ v 2 = 7v 2 + v 1 . Mas recordem que a nulidade de A — 7Id era 3. Este v x é um vetor do 
núcleo de A — 7Id; podemos encontrar mais dois: serão o v 3 e o v 4 . Ora, escrevendo isto 
em termos de matrizes, ficamos com algo como 

A • [ v : v 2 v 3 v 4 ] = [ A • A • v 2 A • v 3 A • v 4 ] = [ 7v x 7v 2 + v 1 7v 3 7v 4 ] 

7 10 0 
0 7 0 0 
0 0 7 0 ’ 

0 0 0 7 

que corresponde à segunda matriz mostrada acima. 

Como imaginam, tentar descrever os outros casos em termos genéricos fica rapidamente 
monótono ou confuso. É melhor passarmos a alguns exemplos concretos. 


= [ Vj v 2 v 3 v 4 ] • 


(3.19) Primeiro exemplo. Começamos com um caso relativamente simples, o da matriz 


A = 


16 0 
0-5 0 

-12 -12 -5 


O polinómio característico de A é (À — 1) (À + 5) 2 . 

Os vetores próprios com X = 1 são obtidos resolvendo o sistema (A — Id) • v = 0: 


0 

6 

0 

0 


1 

0 

-6 

0 

0 

, cujas soluções são os múltiplos de v l = 

0 

-12 

-12 

-6 

0 


-2 


Para X = — 5, podemos verificar que A + 5Id tem nulidade 2. Como coincide com a 
multiplicidade algébrica, podemos calcular imediatamente os vetores próprios, resolvendo 
(A + 5 Id) • v = 0: 


6 6 0 

0 

ir /6 

110 

0 

0 0 0 

0 


0 0 0 

0 

-12 -12 0 

0 

L 3 + 2 L 1 

0 0 0 

0 


As soluções deste sistema são dadas por 


Xi + x 2 = 0 => 

*i 

*2 

= *2 

-1 

1 

+x 3 

0 

0 


*3 


0 


1 


= x 2 v 2 + x 3 v 3 . 


Ou seja, encontrámos dois vetores próprios independentes com X = —5. 
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Dissemos antes que se v é um vetor próprio de A com valor próprio A, então y(t) = e At v 
é uma solução de y/(t) = A • y(t). Como encontrámos três vetores próprios independentes 
( y i, v 2 , e v 3 ), encontrámos também três soluções: 


Yi(t) = e t v 1 = e t 

1 

0 

, Y 2(0 = c St v 2 = e 5t 

1 - 

1—1 i—l 

1 


-2 


0 


e 


y 3 (t) = e 5t v 3 = e 5t 


0 

0 

1 


Construímos uma matriz Y(t) com estas três soluções como colunas: 


nO=[yi(t) Y2(0 Y3(0 ] 


e t —e 5t 0 
0 e~ 5t 0 

—2e t 0 e~ 5t 


Observando que a matriz Y(0) = [ v 3 v 2 v 3 J é invertível (pois já vimos que os três ve¬ 
tores são, de facto, uma base) e recordando o que dissemos sobre a exponencial de matrizes, 
concluímos que e At = Y(t) • Y(0) _1 . Ou seja, obtemos 



O 

LO 

1 

CD 

1 

<D 


1-10 

II 

CD 

0 e~ 5t 0 


0 10 


1 

LO 

1 

CD 

O 

CD 

CN 

1 


-2 0 1 


Há outra forma de prosseguir após a determinação dos vetores próprios. Consiste em 
observar que, se construirmos uma matriz S = [ v 3 v 2 v 3 J com os três vetores próprios 
como colunas, então 

A • S = A • [ Vj v 2 v 3 ] = [ A • v 3 A • v 2 A • v 3 ] 


[ v 3 — 5v 2 — 5 v 3 ] = [ Vj v 2 v 3 ] • 


1 

0 

0 


0 

-5 


0 

0 


= S ■ D. 


D 

(Este é o raciocínio que está por trás da definição tradicional de diagonalização de matrizes.) 
Como A- S = S ■ D, temos ■ S = S ■ e Dt . Mas como D é uma matriz diagonal, já sabemos 
que 


e f 0 0 

0 e~ 5t 0 
0 0 e~ 5t 


Portanto, = S ■ e Dt ■ S 1 , ou 


1-10 


e f 0 0 


1-10 

0 10 


0 e~ 5t 0 


0 10 

-2 0 1 


0 0 e~ 5t 


-2 0 1 


Se fizerem o primeiro produto de matrizes, podem ver que obtemos a mesma expressão que 
tínhamos antes. (E portanto, podem escolher qual dos dois métodos preferem usar.) 
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(3.20) Segundo exemplo. Quando a matriz é diagonalizável, a situação é sempre seme¬ 
lhante ao que fizemos com a matriz A. Mas o que se passa se a matriz não for diagonalizável? 
Consideremos a matriz 


B = 


-8 1 
-4 -12 


cujo polinómio característico é (A +10) 2 . Temos portanto um único valor próprio, À = —10, 
com multiplicidade algébrica 2. Como a matriz 


B-ÃId = B + 10Id = 


2 1 
-4 -2 


tem nulidade 1 (lembram-se como calcular a nulidade?), temos de considerar (B + 10 Id) 2 , 
cuja nulidade é 2 (pois trata-se da matriz nula). 

Escolhemos então qualquer v 2 no núcleo de (B + 10Id) 2 . Por exemplo, v 2 = [5] (para 
quê complicar?). Em seguida definimos v x = (B + 10 Id) • v 2 = [ _ 2 ] . 

Essas escolhas garantem que (B + 10 Id) • v 1 = {B + 10 Id ) 2 • v 2 = 0. Por conseguinte, 
B ■v 1 = —lOvj e B ■ v 2 = —10v 2 + v 1 (verifiquem!). Usando esses dois vetores, obtemos 
duas soluções de y'(t) = B ■ y(t): 

y 1 (t) = e- 10 t v 1 = e- 10t [ * 


y 2 (t) = e 1 ° t (v 2 + tv 1 ) = e 


-10t 


t 

1-21 


Se construirmos a matriz T(t) = y^t) y 2 (t) ] , obtemos 
e Bt = T(t) • T(O ) -1 = e _10t 


1 t 
-2 1 — 2 t 



1 

o 

T— 1 


1 

T— 1 

CN 

1 


-i -i 


Alternativamente, depois de termos escolhido v, e v 2 , podemos obter a forma de Jordan 
para B. Pondo S = [ v 1 v 2 J , temos 


8-S=B-^v 1 v 2 J = [^B-v 1 B-v 2 J 

= [ -lOv! 10v 2 + v a ] = [ v a v 2 ] 


—10 


0 -10 


= S ■ J. 


A exponencial de J é conhecida, por isso podemos dizer 


e Bt = S ■ e Jt ■ S -1 = 


1 0 

-2 1 


10t 


1 t 
0 1 


1 0 

-2 1 


-i -i 


Se fizerem as contas, podem confirmar que chegamos de facto ao mesmo resultado que no 
parágrafo anterior. 


(3.21) Terceiro exemplo. Experimentemos agora com a matriz 


C = 


-8 13 

0-2 0 
-12 2 4 
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cujo polinómio característico (verifiquem) é (À + 2) 3 . Temos portanto um único valor pró¬ 
prio, A = — 2, com multiplicidade (algébrica) 3. 

Como a nulidade de 


C + 2 Id = 


-6 13 
0 0 0 
-12 2 6 


é 2 (que é menor que a multiplicidade algébrica, 3), temos de considerar (C + 2 Id ) 2 que 
tem nulidade 3 (pois é a matriz nula). 

Queremos escolher um v 2 no núcleo de (C + 2 Id ) 2 mas fora do núcleo de C + 2Id—de 
contrário v 1 = (C + 2Id) • v 2 seria 0, e teríamos de voltar atrás para escolher outro v 2 . Uma 
das (imensas) escolhas possíveis é 



0 


1 

V 2 = 

1 

e Vj = (C + 2 Id) • v 2 = 

0 


0 


2 


Com isto, encontrámos um dos dois vetores próprios independentes (são dois, porque a 
nulidade de C + 2Id era 2). Precisamos de encontrar o outro (que será o v 3 ). Para tal, 
precisamos das soluções de (C + 2 Id) • v = 0. Através do método de Gauss, 


-6 13 

0 


-6 13 

0 

0 0 0 

0 


0 0 0 

0 

-12 2 6 

0 

L 3 -2L, 

0 0 0 

0 


vemos que as soluções correspondem a — óxq + x 2 + 3x 3 = 0. Entre as soluções, temos que 
escolher uma independente de Vj e v 2 . Pode ser, por exemplo. 


v 3 = 


0 

-3 

1 


O resultado de tudo isto é que encontrámos três vetores linearmente independentes. 
Como (C + 2Id) • Vj = 0, temos C • v x = — 2v x . Como (C + 2 Id) • v 2 = v 1; temos C • v 2 = 
—2v 2 + Vj. Como (C + 2 Id) • v 3 = 0, temos C ■ v 3 = — 2v 3 . 

A estes três vetores podemos associar soluções de 5 /(t) = C • y(t): 


CM 

1 

O) 

II 

> 

CM 

1 

II 

1 

0 

, y 2 (t) = e 2 t (v 2 + tv 1 ) = e 2t 

t 

1 


2 


21 


e y 3 (l) = 6 _ 2 t v 3 = e~ 2t 
Como no exemplo anterior, construímos Y(t) = [ yi(t) y 2 (t) y 3 (t) J e obtemos 


e ct = y(t) • y(o) _1 = 

e~ 2t te~ 2t 0 

0 e~ 2c —3e~ 2t 


1 0 

0 1 

1 - 

O CO 

1 


2e~ 2t 2 te~ 2t e~ 2c 


2 0 

1 


-1 
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Alternativamente, podemos construir explicitamente uma forma de Jordan para a matriz 
C. Com efeito, pondo S = [ v x v 2 v 3 J , temos 

C • S = C • [ v a v 2 v 3 ] = [ C • Vj C • v 2 C • v 3 ] 



-2 

1 

0 

2v 1 -2v 2 + v! — 2 v 3 ] = [ v 1 v 2 v 3 ] 

0 

-2 

0 


0 

0 

-2 


j 


Ora, como a matriz J é formada por dois blocos de Jordan, temos 


e Jt = 

ei o —2 ] f 

o o 


1 

1 

to 

1 

0 

i i 

4-» i—1 


O o 


e 2t 

0 

te 2t 
e~ 2t 

0 

0 


1 

O 

o 

e [-2]t 


O 

o 

e 2t 


1 

O 

0 

e -2t 


Temos então e Ct = S-e Jt ■S 1 . Se fizerem o primeiro produto, podem ver que esta expressão 
para e Ct realmente coincide com a que obtivemos acima. 


(3.22) Quarto exemplo. Mais um exemplo, agora com a matriz 


F = 


1 

3 


-12 


0 0 
-3 1 

0 -3 


O seu polinómio característico é (A — 1) (A + 3) 2 . 

Os vetores próprios com A = 1 são os vetores no núcleo de F — Id. Resolvendo o sistema 


0 

0 

0 

0 


1 

3 

-4 

1 

0 

, obtemos os múltiplos de v 3 = 

0 

-12 

0 

-4 

0 


-3 


Como a nulidade de 


F + 3 Id = 


4 0 0 
3 0 1 
-12 0 0 


é 1 (menor que a multiplicidade algébrica 2), temos que considerar 


(F + 3Id) 2 = 


16 0 0 
0 0 0 
-48 0 0 


Esta matriz tem nulidade 2. Vamos então escolher qualquer v 3 no seu núcleo, exigindo 
apenas que v 2 = (F + 3 Id) • v 3 seja não-nulo. Por exemplo, tomemos 


v 3 = 

0 

0 

e v 2 = (F + 3 Id) • v 3 = 

0 

1 


1 


0 
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(Se, por coincidência, obtivéssemos v 2 = 0, teríamos que voltar atrás para escolher outro 
v 3 .) 

Encontrámos assim três vetores linearmente independentes. Como (F — Id) • v 1 = 0, 
temos F • Vj =v 1 . Como (F + 3 Id) • v 2 = 0, temos F • v 2 = —3v 2 . E como (F + 3 Id) ■ v 3 = v 2 , 
temos F • v 3 = —3v 3 + v 2 . 

A estes três vetores podemos associar soluções dey / (t) = F-y(t): 



T—1 


O 

_1 

O) 

II 

O) 

II 

0 

00 

1 

II 

CM 

> 

00 

1 

O) 

II 

CN 

1 


1 

CO 

1 

_1 


-1 

o 

_1 


e y 3 (0 = e 3t (v 3 + tv 2 ) = e 3t 


0 

t 

1 


Como nos outros exemplos, construímos y(t) = [ Yi(t) y 2 (t) y 3 (t) ] e obtemos 



Oi 

O 

O 


10 0 

e Ft = y(t) • y(o) _1 = 

0 e“ 3t te~ 3t 


0 10 


—3e c 0 e“ 3t 


-3 0 1 


Alternativamente, podemos construir explicitamente uma forma de Jordan para a matriz 
F. Pondo S = [ v 3 v 2 v 3 J , temos 

F • S = F • [ Vj v 2 v 3 ] = [ F • Vj F • v 2 F • v 3 ] 



1 

0 

0 

[ Vi —3v 2 -3v 3 + v 2 ] = [v 1 v 2 v 3 ] ■ 

0 

-3 

1 


0 

0 

-3 


j 


Como a matriz J é formada por dois blocos de Jordan, temos 




o 

o 



0 0 

e Jt = 

1 

O o 

eii-íl 

— 

1 

O o 

e~ 3t 

" 1 t 

0 1 



Temos também e Ft = S ■ e Jt ■ S 1 . 


e f 0 0 

0 e“ 3t te~ 3t 
0 0 e“ 3t 


(3.23) Quinto exemplo. Para um último exemplo, consideramos a matriz 


G = 


3 0 0 

3 0 1 

10 -9 6 


cujo polinómio característico (verifiquem) é (A — 3) 3 . Vemos que G tem um único valor 
próprio, A = 3, com multiplicidade (algébrica) 3. 
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Como a nulidade de 


G — 3 Id = 


0 0 0 
3-3 1 
10 -9 3 


é 1 (confirmam?), teremos de considerar (G — 3 Id) 2 . Mas esta matriz tem nulidade 2 
(concordam?), que é ainda inferior à multiplicidade algébrica. Por isso, teremos que ir até 
(G — 3 Id) 3 , que é a matriz nula. 

Assim, escolhemos um qualquer v 3 e definimos v 2 = (G — 3 Id) • v 3 e v 3 = (G — 3 Id) • v 2 , 
exigindo apenas que resulte v 3 ^ 0. Por exemplo, tomemos 



0 


0 


0 


0 


0 

V 3 = 

0 

1 

, v 2 = (G — 3 Id) • 

0 

1 

— 

1 

3 

, e v 3 = (G — 3 Id) • 

1 

3 

— 

0 

0 


Azar—chegámos a v 3 = 0. Tentemos de novo: 


T— 1 

1 _ 


T— 1 


O 


O 


O 

_ 1 

0 

0 

, v 2 = (G—3Id> 

0 

0 

— 

3 

10 

, e v 1 = (G—3Id)- 

3 

10 

— 

1 

3 


Encontrámos três vetores linearmente independentes. Como (G — 3Id) • v 3 =0, temos 
G • v 3 = 3v 3 . Como (G — 3Id) • v 2 = v 1; temos G • v 2 = 3v 2 + v 2 . E como (G — 3Id) • v 3 = v 2 , 
temos G • v 3 = 3v 3 + v 2 . 

A esses vetores podemos associar estas soluções de y^t) = G • y(t): 



0 


0 

y 1 (t) = e 3 t v 1 = e 3t 

1 

, y 2 ( t ) = e 3 t (v 2 + tv 1 ) = e 3t 

3 + t 


3 


10 + 3t 


e ys (0 = e 3t (v 3 + tv 2 + ^v 1 ) = 


„ 31 


3t + J 
10t + 3| 


Como nos exemplos anteriores, construímos uma matriz Y(t) a partir destas três soluções, 


no = 


0 

„3t 


0 


„3t 


(3 + t)e 3t (3t + y) e 3f 


3e 3t (10 + 3t)e 3í (I0t + 3|-)e 3t 


e concluímos que e Gt = T(t) • T( 0 ) 1 . 

Ou então construímos a matriz S = [ v 1 v 2 v 3 J e calculamos 


G ■ S = G ■ [ 


vi v 2 v 3 = 


] = [ G • v 3 G • v 2 G • v 3 ] 


[ 3v 3 3v 2 +Vj 3v 3 + v 2 ] = [ v 3 v 2 v 3 ] • 


3 10 
0 3 1 
0 0 3 


= S ■ J. 
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Então, como e Gt = S ■ e Jt ■ S 1 e J é um bloco de Jordan, teremos 


.Gt _ 


= S ■ e 


3 1 


1 

0 1 
0 0 


f 2 
t 

1 


(3.24) Exercício. Determinem as formas de Jordan e as exponenciais das seguintes matrizes: 


(a) 

"-9 9" 
-4 3 

5 


(b) 

8 

-9 

1 " 
2 

} 

(c) 

i 

-3 

4 

-2 

5 

-7 

-2 

5 

-7 


-1 -1 

2 



0 

2 

-3 


-1 

1 

-2 

(d) 

-4 0 

-5 


; (e) 

4 

-2 

-2 

; (f) 

-3 

-1 

6 


-4 2 

-7 



4 

-4 

0 


-1 

-1 

4 


1 2 

2 



-2 

7 

-1 





(g) 

-3 6 

3 

í 

(h) 

-3 

7 

-1 






1 -1 

2 



-1 

1 

1 






Os valores próprios da matriz (c) são 0 e —1. Os valores próprios da matriz (d) são —2 e —3. 
Os valores próprios da matriz (e) são ±2. As matrizes (g) e (h) só têm um valor próprio. 


(3.25) Exercício. Resolvam os seguintes sistemas de equações diferenciais com t como va¬ 
riável independente: 

(a) x^ = -9x a + 9x 2 , x ' 2 = -4x a + 3x 2 , x a (0) = 7, x 2 (0) = 5; 

(b) x[ = 8 x a + x 2 + 1, x' 2 = —9x 1 + 2x 2 - 3, x 1 (0) = -3, x 2 (0) = 10; 

(c) Xj = —2x 1 + 7x 2 - x 3 - 2, x 2 = —3x x + 7x 2 — x 3 — 1, x' 3 = — x 3 + x 2 + x 3 + 1, 

x 1 (0) = l, x 2 ( 0 ) = 0 , x 3 (0) = — 1 ; 

(d) x\ = —2x 1 + 7x 2 *- x 3 + 1, x' 2 = -3x a + 7x 2 - x 3 , x' 3 = -x 1 + x 2 + x 3 - 1, 

x 1 (0) = -2, x 2 (0) = —1, x 3 (0) = 1; 


(e) x[ = -2x 1 + 7x 2 - x 3 + 1, 
x 1 {2) = -2, x 2 (2) = —1, 


x 2 = 3 x 3 + 7x 2 x 3 , 

x 3 (2) = 1. 


X 3 = -x 3 +x 2 + x 3 - 1, 


(3.26) Coeficientes não constantes. É razoável perguntar o que se passa quando os coe¬ 
ficientes da equação não são constantes, como em 

y / (t)=A(t)Y(t)+B(t), Y(t 0 ) = y 0 . 

Vamos admitir que exista algum análogo W(t) da exponencial de matrizes, satisfazendo 
W'(t) = A(t) W(t). A solução seria da forma Y(t) = W(t)Z(t)? Se sim, teríamos 

y / (t) = w'(t ) z(t) + w(t) z\t) = A(t) w(t) z(t) + w(t) z'(t) = A(t) y(t) + w(t) z\t). 

Por outro lado, Y'(t) = A(t) Y(t) + B(t), o que nos leva a concluir que B(t) = W(t) Z / (t). 
Assim, desde que W(t) seja invertível, temos Z 7 (t) = W(t ) _1 B(t) e conseguimos obter Z(t) 


47/99 




















João Pedro Boavida, Episódios de Equações Diferenciais, 2015 


por primitivação, usando ainda Z(t 0 ) = W(t 0 ) 1 Z(t 0 ) (de onde vem isso?): 

Y(t) = W(t) (w(t 0 )~ 1 Y 0 + f W(r) _ 1 B(r)dr). 

J t 0 

(De acordo?) É uma simples variante da fórmula da variação das constantes, desta vez com 
a W(t) em vez da exponencial. 

Para construir W(t) precisamos de várias soluções (uma para cada coluna) de w'(t) = 
A(t)w(t), e precisamos que sejam linearmente independentes em t 0 (para assegurar que 
W(t) é invertível numa vizinhança de t 0 ). Ou seja, construímos W(t) a partir de soluções 
do problema homogéneo associado à mesma matriz A(t). 

Vejamos um exemplo. Considerem o sistema 

í (i-t 2 )yí = -tji + y2 + t(i-t 2 ), ji(o) = o, 

\ (i - t 2 )y ' 2 = yi~ ty 2 + (i -1 2 ), j2(o) = i- 

O sistema homogéneo associado é 

í (1 — t 2 )w / 1 = —tw 1 + w 2 , 

\(1 -t 2 )w' 2 = w 1 -tw 2 , 


e podemos verificar que [ j ] e [ j ] são duas soluções linearmente independentes (desde 
que 1 — t 2 ^ 0). Definimos então W(t) = [ t í ] ■ 

O sistema original pode ser reescrito (se aceitarmos uma divisão por 1 — t 2 ) na forma 


y / (0 = 


" -t 1 ' 

1 

t 

1 -t 

1 -U y(t) + 

1 


(3.27) Exercício. Qual a solução y(t)? Qual o maior intervalo em que é solução da equação 
original? Qual o intervalo em que W(t) é invertível? 


(3.28) Exercício. Considerem a equação 


y'(t) = 


o o 
2 1 0 


y (0 + 


3t 

0 


y(o) = 


i 

-i 


Sabendo que (1, t 2 — 1) e (1, t 2 + 1) são soluções do sistema homogéneo, determinem y(t). 

(3.29) Exercício. Considerem o sistema 

f zí (0 = tz 1 (t)-z 2 (t) + l, 


z' 2 ít) = (1 + t 2 )Zl(t) - tz 2 (t ) + t, 


z i(0) = z 2 (0) = 0. 

Sabendo que (t, t z — 1) e (1, t) são soluções do sistema homogéneo, determinem z(t). 


(3.30) Exercício. Considerem o sistema 

í w'(t) = 2 w 1 (t) - e l w 2 (t), 
y w' 2 (t) = 2 e -t w 1 (t) — 2 w 2 (t) + 1 , 

Sabendo que (e £ , 1) e (1,2e _t ) são soluções do sistema homogéneo, determinem w(t). 


w 1 (0) = w 2 (0) = 1. 
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(3.31) Exercício *. Vamos considerar a equação geral 5 /(t) = A(t)y(t) +b(t) definida em 
a < t < b, onde cada A(t) é uma matriz n x n com componentes a í; (t). Vamos usar o 
teorema de Picard-Lindelõf, com a função F : t,x >-» A(t) x + b(t). 

(a) Mostrem que FC^yq) - F(t,y 2 ) = A(t) (y x - y 2 ). 

(b) Suponham que cada função a l j é limitada em a < t < b. Mostrem que existe C > 0 
tal que, para a<t <b ey 1 ,y 2 ei", temos llFCt^) - F(t,y 2 )|| < C ||y x -y 2 ||. 

(c) Se existir C como na alínea anterior, mostrem que existe uma única solução y satis¬ 
fazendo a condição inicial y(t 0 ) = Yo> onde a < t 0 < b. 

(d) Mantendo as hipóteses das alíneas anteriores, mostrem que podemos escolher n 
condições iniciais linearmente independentes. 

(3.32) Exercício. Seja agora a equação homogénea w 7 (t) = A(t) w(t) em a < t < b, onde 
cada A(t) é uma matriz n x n, cujas componentes são limitadas em a < t < b (como no 
exercício anterior). Consideramos ainda n soluções w l3 ..., w n dessa equação, bem como a 
matriz W(t) com colunas w^t),..., w„(t). 

(a) Fixado um vetor v e R n , mostrem que t —> W(t)v é solução da equação homogénea. 

(b) Se W(t) v = 0 para algum a < t < b, mostrem que W(t) v = 0 para qualquer outro 
a < t < b. 

(c) Se {w 1 (t),..., w n (t)} é linearmente independente para algum a < t < b, mostrem 
que o é para qualquer outro a < t < b. 

(d) Mostrem como construir, a partir de n vetores linearmente independentes, soluções 
w l3 ..., w n de modo que a matriz W(t) correspondente seja invertível para todo o a < t < b. 
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(3.33) Respostas de exercícios. 


031) (— + —- —e~ Wt + —e~ 10t 

tiU-l v 100 10 100 5 


2 _ t , 43 -10t 

50 5 ^ 50 E 


6 t -10t 
5 


)■ 


( —-— - H—— 

V 49 7 4Q ’ 


49 


£ )- 


(d) 


(e) 


-2 0 0 
0-3 1 

0 0-3 


(l+ 2 t)e“ 3t 

e~ 3t —e~ 2t 

(t—l)e~ 3 í +e~ 

_4 te -3t 

3e~ 2t —2e~ 3c 

(3—2t)e" 3 í —3e 

-4te _3t 

2e~ 2t —2e~ 3t 

(3—2t)e _ 3 í —2e 

‘ 2 0 

0 ' 



0-2 1 



O 

o 

-2 . 


[ 2e 8t —e 5t e s ‘-e 5t 1 

" (l+ 2 t)e“ 2t 

e 2 r —(l+ 2 t)e~ 2f 

(l + t)e - 2 í —e 2t ' 

L — 2 e 8 t + 2 e 5í -e 8 t + 2 e 5 ‘ J ' 

4te -2t 

e 2 ‘-4te- 2 ‘ 

(l+ 2 t)e“ 2 t —e 2r 

Ie 5r -6e 7t +6e 5t 1 

4te _2t 

-4te~ 2 ' 

(l+ 2 t)e -2t 


4e n -3e b 


2e 7t -2e 5c 
—4e t 
—4e t 

—4e - 4 t +5e~ 5í 4e" 4t - 
—5e - 4 t +5e~ 5t ^ _4t - 


e t +5e~ t 5e t -5e“ c 1 
e t +4e~ t 5e t —4e~ t J ' 


(b) 

(c) 

(d) 


dÕTfl ) (a) ( 0 , 0 ). 

(c) segunda coluna da exponencial. 

fcD ÍÍ-Í + Í4. -?§_ í-lj-A + -26_) 

W V 5 25 ^ ^ 25e St 5 5 25 ^ ^ 25c st J 


,5t J- 

to 

o 

o 

_1 


1 -t t 

—2 t 

1 (f) 

0 0 1 

? 

—3t 3+3t—2e 2t 6e 2t —6—6í 

• 

o 

o 

o 


-t 1+t—e 2t 3e 2t —2—2t 

J 

"3 1 0" 


' 1—2t 2t 2t ' 


-*1:1. (g) 

0 3 0 

j 

— 3t l+3t 3t 

e 3t . 

-4e~ 5 J 

_0 0 3 _ 


t -t 1-t. 



'2 1 0" 


l-4t-2t 2 7t+3t 2 -t- 

(h) 

0 2 1 

5 

— 3 t-t 2 l+5t+ 

3 1 2 -t- 

2 

f , — + -t - —e 7t ). 

_0 0 2 


t 2 -t t-| 

2 l-t + 

5 49 7 49 J 





l [3~25l ) (a) (7 + 3t,5 + 2t)e“ 3t . 
(b) (te 5t - ^e 5t - |, ^e 5t - 3 te b 


(c) (1 - 3te 2t - t 2 e 2t , 1 - ±e 2t - 2te 


(HHU) (a) [; ;]e 3 


+ !)• 

2 t 



■e 5t 

0 

0 


" 1 

1 — 

o 

o 

(b) 

0 

e 4t 

te 4 ‘ 

• (c) 

0 

1 t 


_ 0 

0 

e 4r . 


_0 

0 1 _ 


-t 2 e 2t 

2 ’ 


e 2t + h 2 e 2t ). 


te 


l r 2 2t 3 , 5 2t , 

4 ’ 8 ' H 




3 _ 5 .21 
• 8 8 
jte 2t + i 

4 4 


+ Tt 2 e 2t 


)■ 


ÜH) W (i-^ 3t + ke 3t ,-i + ie 3t ). 


(b) (e 5t , te 4t 'e 4t ). 

(c) (t — 1 + e _t , e _t + 2 te _t , 2 e _t ). 

ra> (a)[- 0 3 _ 1 3 ];[ 1 :í t i^]e- 

(b) 

(c) 


r 5 j i- r 

!+ 3t f 1 .5t 

Lo 5 J’ L 

—9t 1 —3t J ‘ 

1 

1 

0 

0 


' 2—e~ t 2e~‘-2 2eT t -2 ' 

0 0 1 

3 

2e~ f —2—t 5+t—4e~ £ 4+t-4e~ f 

0 

0 

0 


_ 3+t—3e _t 6e~ t — 6 — t 6 e _í —5— t _ 


(d) (-1 ■ 

lte 2t 

4 

(e) igual, mas com t substituído por t — 2. 
( f3T27l ) y(t) = (t + t 2 ,1 + t). R. ]— 1,1[. 
( f^28l ) y( 0 = (1 + |t 2 , -1 + t 2 + |t 4 )- 
( 13391 ) z(t) = (t,t 2 ). 

( [37301 ) w(t) = (—1 + 2e c — te 1 ,3— t — 2e _t ). 
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EPISÓDIO 4 


EQUAÇÕES DIFERENCIAIS DE ORDEM SUPERIOR 


As equações diferenciais de ordem superior podem ser resolvidas com as técnicas que discu¬ 
timos no episódio [ÕJ Por exemplo, f" — 13f' + 42/ = 0 (a variável independente continua 
a ser t) pode ser simulada usando x 1 = f e x 2 = f' ■ O sistema fica então 


tN 

* 

II 

II 

b-T 

*i 

/ 

0 1 “ 


*i 

\ , ,, , ou 

= 

-42 13 


{ X 2 = f = -42/ + 13 /' = —42x 1 + 13x 2 , 

. * 2 . 



L *2 J 



A 


Esta é exatamente a equação discutida no início desse episódio. O polinómio característico 
da matriz A é — A(13 —A)+42, ou (A—7) (A—6). Quais os vetores próprios correspondentes? 
Digamos queAv 1 = 7\\ e Av 2 = 6v 2 . Assim, 


A[vi v 2 ] y(t) = 

■V 

S 


m 
r (o 


= |~ M - Vi M ■ v 2 ] = [ 7v 1 6v 2 ] = S 


7 0 
0 6 


e podemos concluir que 
e* =S-e Dt -S~ 1 =S- 


,7t 


0 e 


6 1 


■ s~ 


x(t) = S ■ 


,7t 


0 e 


6t 


■S _1 -x(0). 


Em particular, vemos que a primeira coordenada da solução (correspondente a Xj = /) é 
uma combinação linear de e 7c e e 6t . 

Isto encaixa num padrão mais geral: se usarmos D para indicar a derivação, podemos 
escrever f" — 13/ '+ 42/ = (D 2 —13D+42)/ e fatorizar D 2 — 13D+42 = (D — 7) (D — 6). Por 
outro lado, podemos ver que (D — 7) aniquila e 7t , ou seja, que (D — 7) (e 7t ) = (e 7t y — 7e 7t = 
0. Da mesma forma, (D — 6) aniquila e 6t . Vemos então que as funções aniquiladas por 
D 2 — 13D+42 (entre as quais está a solução / da nossa equação) são uma combinação linear 
das funções aniquiladas pelos seus fatores (que neste caso não têm raízes em comum). O 
método dos aniquüadores sistematiza esta estratégia, e será o nosso primeiro assunto neste 
episódio. 


(4.1) Exercício. Escrevam a equação matricial correspondente a uma equação como 

f" + af" + bf' + cf = 0. 

Qual o polinómio característico da matriz que encontram? (Calculem o determinante 
usando a regra de Laplace em relação às primeiras colunas das matrizes.) 
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(4.2) Vamos olhar para a equação f"+f— 2/ = g, com condições iniciais/(O) = /'(O) = 3. 
Escrevendo-a na forma matricial, obtemos 


' / " 

/ 

’ 0 1 ' 

' / ' 


’ 0 " 


' /( 0 ) ' 


’ 3 ‘ 

. f'. 


2 -1 

. f' . 

+ 

. g _ 




3 




y(o) 


Pela fórmula da variação das constantes, sabemos que 


y(t) = e At 




rt 


e^y(0) 

solução da eq. 
homogénea 



e A 'b(r)dr. 


J o 

-V-- 

solução particular 


(4.3) Exercício. Calculem a exponencial da matriz, e usem-na para resolver a equação no 
caso g(t) = 6e f + 2. 

(4.4) Só nos interessa a primeira linha. Vemos que, se não tivermos condições iniciais, 
temos de escolher um vetor constante em M 2 (ou seja, a primeira linha envolve tantas 
constantes quanto a ordem da equação) e que há sempre uma solução particular. 

Porém, também vemos que, já que as soluções envolvem exponenciais (que podem ser 
determinadas a partir dos valores próprios da matriz, ou mesmo diretamente a partir da 
equação—como?), seria desejável evitar a fórmula da variação das constantes. 


(4.5) Derivadas e aniquiladores. Começamos por identificar o núcleo de (D — À). Dizer 
(D — X)f = 0 é dizer f — Xf = 0. Esta equação não é exata, mas e~ Xt é um fator integrante. 
Assim, podemos escrever 

f'-Xf = 0 <=} e~ Xt f-e~ Xt Xf =0 <=} (e _At /) / = 0. 

Ou seja, e~ Xt f é uma constante C, ou / = Ce Xt . 

Imaginem agora que temos qualquer função / (aniquilada ou não por (D - 
queremos calcular g = (D — X)f = f 1 — Xf. Usando o mesmo fator integrante e 
de contas, continuamos a ter uma equação linear), obtemos 

e~ Xt g = e~ Xt f' — Xe~ Xt f = (e~ Xt ff ^ g = e Xt (e“ At /)'. 

Concluímos assim que (D — X)f = e Xt (e _At /) / . 

(4.6) Exercício. Mostrem que (D — X) k f = e Xt {e~ Xt (para k inteiro positivo). Com 
base nesse resultado, mostrem que (D — X) k f = 0 (i.e., / está no núcleo de (D — X) k , i.e., 
(D — X) k aniquila /) precisamente se e~ Xt f é um polinómio P(t) de grau menor que k, ou 
seja, se /(t) = P(t) e Xt . 

(4.7) E o que acontece se o nosso operador tiver valores próprios distintos? Em geral, pode¬ 
mos escrever a equação original na forma ABf = 0, onde A e B são operadores envolvendo 
apenas a derivada D (e portanto Ae B comutam) e sem valores próprios em comum. Nesse 


Ã)) e que 
Xt (no fim 
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Episódio 4: Equações diferenciais de ordem superior 


caso, todas as soluções de ABf = 0 são da forma f = a + b, onde Aa = 0 e Bb = 0 (isto não 
é óbvio; conseguem demonstrar?). 

A ideia ficará mais clara com um exemplo. Imaginem que queremos encontrar as solu¬ 
ções de (D + 2)(D — 1 ) 2 D f = 0. Como os fatores (D + 2), (D — l) 2 , e D não têm valores 
próprios comuns (os valores são, respetivamente, — 2, 1, e 0), sabemos que / é uma combi¬ 
nação linear de elementos nos núcleos de (D + 2), (D — l) 2 , e D. Ora, o núcleo de (D + 2) 
consiste nos múltiplos de e~ 2t . O núcleo de (D — l) 2 consiste em polinómios de grau menor 
que 2 multiplicados por e l . E o núcleo de D consiste nas constantes. Assim, vemos que / é 
da forma 

/ = ^i e _^ + C 2 e t + C 3 t e^ + C 4 

núcleo de (O + 2) núcleo de (D - l) 2 núcleo de D 


(4.8) Podemos querer o contrário: dada uma função, encontrar um operador que a aniquile 
(ou seja, encontrar uma aniquüador ). Para tal, começamos por juntar as parcelas associadas 
ao mesmo valor próprio (ou seja, com a mesma exponencial), e determinar, para cada valor 
próprio, qual o grau necessário. 

Por exemplo, quem é um aniquilador para e c — t 5 e 5t + 3t 2 e l — t 3 e 4t + 10e 5t —4 te 1 ? Se 
juntarmos as parcelas associadas a cada valor próprio, podemos ler diretamente: 

(1 + 3t 2 — 4t) e f + (—t 5 + 10)e 5t + (—t 3 ) e 4t 

aniquilado por (D - l) 3 aniquilado por (D - 5) 6 aniquilado por (D — 4) 4 

Assim, (D — l) 3 (D — 5 ) 6 (D — 4) 4 (ou qualquer outro operador com os mesmos fatores) é 
um aniquilador para essa função. 

Um ponto importante é que só podemos aplicar este método depois de termos simpli¬ 
ficado cada parcela de modo a ter apenas uma exponencial e um polinómio. Por exemplo, 
se tivermos um produto de exponenciais, temos de simplificá-lo; se tivermos um produto 
de exponenciais e funções trigonométricas, temos de converter as funções trigonométricas 
para exponenciais complexas e simplificar. 

(4.9) Exercício. Determinem o aniquilador mais simples de cada uma das seguintes funções: 


(a) e 4t ; 

(b) 3e 4t ; 

(c) e 3t ; 

(d) e 3t + 3e 4t ; 

(e) e lt ; 

(f) e~ ic ; 

e it + e~ it 


Cg) 2 ; 

(h) cos t; 



2it -2 it 


(i) te 3t+2it 

; © te 3t ~ 2it ; 

Qf C b 

(B «* 2 . ; 

(1) te 3í cos(2t) 

(m) te 3t - 

t 2 e 3t ; (n) te 3t + t 2 e 3t sin(2t); 

(o) sin(2t) + t 3 cos(2t); 

(p) (sin(2t)) 2 . 


(4.10) Vale a pena tomar nota deste padrão: 

(D - a) 2 + b 2 = (D — a) 2 - (bi) 2 = (D - a - bi)(D - a + bi ) 

aniquila combinações lineares de e^ a+bl ' ,t e e ( a-bl )£, q Ue coincidem (porquê?) com as com¬ 
binações de e at cos(bt) e e at sin(bt). 
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Da mesma forma, ((D — a) 2 + b 2 ) k = (D — a — bi) k (D — a + bi) k aniquila produtos 
de polinómios de grau menor que k por essas exponenciais, ou seja (porquê?), aniquila 
combinações lineares de parcelas das formas t ( e at cos(bt) e t f: e at sin(bt) (onde í < k). 

(4.11) Exercício. Indiquem três aniquiladores distintos para 

e 2t + a/3 te 2t + (cos t) 2 + t 2 cos(2t) + e c (cos(2t) + t 2 + t 3 sin(3t)). 


(4.12) O método dos aniquiladores. Vejamos então como usar o método dos aniquilado- 
res para resolver a equação 

f"+f'~2f = 6e t + 2, /(O) =/ / (0) = 3. 

Começamos por escrever a equação na forma 

(D 2 + D-2)f =6e f + 2 <=} (D + 2)(D- 1)/ =6e f + 2. 


Neste momento, sabemos resolver imediatamente a equação homogénea correspon¬ 
dente: 

(D + 2)(D-1)/ =0. 

As suas soluções são da forma C 1 e~ 2t + C 2 e t (a solução geral da equação homogénea). 
Nitidamente, as equações homogéneas são mais fáceis de resolver. 

A equação original é não homogénea, pois tem 6e f + 2 no lado direito. Esta função é 
aniquilada por (D — 1)D. Ao aplicarmos este operador a ambos os membros da equação 
original, obtemos 

(D 2 + D — 2)f = 6e f + 2 => (D - 1)D (D 2 + D - 2)f = (D - 1)D (6e f + 2) = 0. 

Notem que os polinómios em D comutam sempre entre si (conseguem demonstrar?). Por 
isso, a nova equação, depois de simplificada e fatorizada, fica 

(D + 2)(D-l) 2 D/ = 0, 


cujas soluções são, como vimos em Q4.7D, 


f — Cie 2t + C 2 e t + C 3 te t + C 4 . 

V -v-' "--^ 

solução homogénea solução particular 


Mas notem que nem todos os valores de C 3 e C 4 são compatíveis com a equação original. 
Na verdade, qualquer outro lado direito (em vez de 6e c H- 2) que tenha o mesmo aniquilador 
dará origem a soluções desta forma. 

Para determinar C 3 e C 4 , substituímos f* = C 3 te t + C 4 na equação original (não vale 
a pena incluirmos parcelas da solução homogénea, pois já sabemos que estão no núcleo de 
(D + 2) (D — 1), ou seja, que depois de as substituir na equação original e simplificar, a sua 
contribuição é 0). Temos então 


(n / = C 3 e t + C 3 te t e (/T = C 3 e t + C 3 e c + C 3 t e c = 2C 3 e t + C 3 te t . 
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Substituindo na equação original — = 6e l + 2, obtemos 

(2 C 3 e l + C 3 t e f ) + (C 3 e f + C 3 t e f ) - 2(C 3 t e l + C 4 ) = 6e f + 2, 

ou simplesmente 3 C 3 e t — 2C 4 = 6e f + 2. Como e‘ e 1 são linearmente independentes, 
podemos concluir que os seus coeficientes não se misturam. Assim, obtemos 3C 3 = 6 e 
—2C 4 = 2, ou C 3 = 2 e C 4 = -1. A solução geral (i.e., ignorando as condições iniciais) da 
equação original é então 

/ =C 1 e“ 2t + C 2 e t + 2te t -l. 

Falta-nos apenas determinar Q e C 2 , e para isso vamos usar as condições iniciais. Como 
f = —2C 1 e~ 2t + C 2 e c + 2e t + 2 te\ 

temos / (O) = C 1 + C 2 — 1 = 3 e / / (0) = —2Cj + C 2 + 2 = 3. A solução desse sistema é C 3 = 1 
e C 2 = 3, pelo que concluímos que 

/ = e~ 2t + 3e f + 2te £ — 1. 


(4.13) Vejamos outro exemplo: 

f"-2f' + 10f = -10, /(0) = 0, f\ 0) = —2. 

O operador no lado esquerdo é D 2 — 2D + 10 = (D — l) 2 + 3 2 , e um aniquilador para a 
função do lado direito é D. Assim, vamos querer resolver 

((D — l) 2 + 3 2 )D f =0. 

Mas as soluções dessa equação são 

/ = Ci e c cos(3t) + C 2 e c sin(3t)+ C 3 . 

V -V- 

solução homogénea /* 

Como ( f*Y = (= 0, substituindo f* na equação original obtemos 10C 3 = —10, ou 
C 3 = — 1. Logo, a solução geral da equação original é 

f = Cie 1 cos(3t) + C 2 e l sin(3t) — 1 


e a sua derivada é 

f' = Ci e l cos(3t) — 3Cj e l sin(3t) + C 2 e l sin(3t) + 3C 2 e c cos(3t), 
pelo que obtemos /(O) = Q — 1 = 0 e /'(O) = C 3 + 3C 2 = —2. Ou seja, Cj = 1, C 2 = —1, e 

/ = e f cos(3t) — e f sin(3t) — 1. 

(4.14) Exercício. Determinem a solução geral de cada uma das seguintes equações: 

(a) f" - 8 f + 15/ = 0; 

(b) f" — 8/ / + 16/ = 0; 

(c) f" — 8f' + 17/ = 0; 
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(d) f" - 8f + 20/ = 0; 

(e) f'" — 3f" + 3f' — / = 0. 

(4.15) Exercício. Determinem a solução geral de cada uma das seguintes equações: 

(a) f" - 8f + 15/ = 3e 3t + 8e f ; 

(b) f" — 8 f + 16/ = 6e 4t — 8 sin t — cos t; 

(c) f" - 8f + 17/ = —10te t + 4e 4t - lle f ; 

(d) f" — 8f' + 20/ = 12e c cos t — 14e c sin t; 

(.é)f'"~ i 3f" + 3f'-f = 6e t . 


(4.16) Exercício. Determinem a solução de cada uma das seguintes equações: 

(a) f" — 8f' + 15/ = 2e 5t , /(0) = 1, f\ 0) = 4; 

(b) /" - 8f + 16/ = 2e st , /(0) = 1, /'(0) = 4; 

(c) f" — 8/ / + 17/ = 2e 5t , /(0) = 1, / / (0) = 4; 

(d) / // -8/ / + 17/ = -3e 4t cos(2t), /(0) = 1, / / (0) = 4; 

(e) /"-8/ / + 20/= -3e 4t cos(2t), /(0) = 1, / / (0) = 4; 

(f) / /// -3/ // + 3/ / -/=6e t , /(0) = / // (0) = -l, / / (0)=1. 

(4.17) Exercício. Determinem a solução de cada uma das seguintes equações: 

(a) f" — 2/ = 4, /(0) = / / (0) = 2; 

(b) f"' + 8/ = 12e _2t , /(0) = -2, / / (0) = -l, /"(0) = 0. 


(4.18) Parte não homogénea definida por trocos. Para resolver uma equação como 

f"-2f'+10f = g, f (0) = 0, /'(0) = -2, 

onde g é uma função definida por troços, temos na verdade de resolver várias equações 
(uma para cada troço). Por exemplo, digamos que 

g(t) = —10, se t < n, e g(t) = 0, se t > n. 

Então a nossa solução terá a forma 

/(0=/i(t), se t < 71, e /(t)=/ 2 (t)> se t > 71, 

onde /j é uma solução de 

/ // - 2 / / + 10 / = - 10 , /( 0 ) = 0 , / / ( 0 ) = — 2 , 
no troço t < n e f 2 é uma solução de 

f" -2f' + 10f = 0 


no troço t > n. Uma vez que a função / tem de ser contínua e ter derivada contínua 
(de outra forma, f" não existiria e a equação original fl4.18[ ) não faria sentido), temos 
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necessariamente _/j(7i) = / 2 (7t) e //(ti) = /^(íi), que nos dá as condições iniciais para o 
segundo troço. 

O primeiro troço coincide com o exemplo ( |4.13D , restrito a t < n. Já vimos (verifiquem!) 
que a solução geral de f" — 2f' + 10/ = —10 é /(t) = C 1 e l cos(3t) + C 2 e t sin(3t) — 1. 
Concretamente, no caso específico em que / (O) = 0 e /'(O) = —2, vimos que a solução é 
fi(t) = e f cos(3t) —sin(3t) —1. Calculando os valores em t = n, obtemos fi(n) = —e n — l 
e fí(n) = 2e n , e por conseguinte / 2 (tt) = — e n — 1 e t) = 2e". 

No mesmo exemplo também vimos que a solução geral de f" — 2 f + 10/ = 0 é dada 
por /(t) = C 1 e f cos(3t) + C 2 e l sin(3t) (onde as constantes podem ou não ser as mesmas 
que na solução anterior). Usando as condições iniciais para f 2 , vemos que C 1 = 1 + e~ n e 
C 2 = — 1 — |e -71 . (Verifiquem!) 

Obtemos assim a solução final 


m= 


e l cos(3t) — e f sin(3t) — 1, se t < n, 

(1 + e _rr ) e l cos(3t) — (1 + |e _7r ) e f sin(3t), se t > n. 


(4.19) Exercício. Considerem ainda 

f" ~ 2f' + 10/ = g, 

e continuem a aproveitar as conclusões do exemplo ( |4.13| ) . 

(a) Qual a solução f 1 para g = 0 e condições iniciais /-/O) = 0 e /j(0) = —2? Mostrem 
que fi(n) = 0 e //(ti) = 2e rt . 

(b) Qual a solução / 2 para g = —10 e condições iniciais / 2 (7t) = 0 e /^(n) = 26^? 

(c) Qual a solução / quando 


g(t) = 0, se t < 7i, e g(t) = —10, se t > ti 
e as condições iniciais são / (O) = 0 e / 7 (0) = —2? 


(4.20) Exercício. Considerem agora 


f" — 8/ / + 15/ = g, / (0) = 1, / / (0) = 4. 


Determinem a sua solução em cada um destes casos: 


(a) g(t) = 


2e 5t , se t < 1, 
0, se t > 1; 


(b) g(t) = 


0, se t < 1, 

2e 5t , se t > 1. 


(4.21) Coeficientes não constantes. Como em ( |3.26D , também aqui é razoável perguntar 
o que se passa quando a equação não tem coeficientes constantes. Olhemos de novo para a 
0. reescrita ligeiramente: 


f' 


2x + l 


f' = 1-4* 2 . 
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Na forma matricial, e chamando y ao vetor, ficamos com 


y'(x) = 


0 

0 


2x+l J 


yO) 


AM 


0 

1 -4x 2 
~ bM 


Pela fórmula da variação das constantes (vamos considerar t 0 = 0), sabemos que 
y(x) = W(x) fw(0) _1 y(0) + í W(r) _1 b(r)drl, 

J o 

onde as colunas de W(t) são soluções do problema homogéneo w'(x) = A(x)w(x), ou seja, 
são obtidas a partir de soluções de f" — ^—f' = 0. 


(4.22) Exercício. Mostrem que fi(x) = 1 e / 2 (x) = x 2 + x são soluções da equação homo¬ 
génea. Qual a matriz W(x) correspondente? Calculem 

= W(0)“ 1 y(0)+ í W(r) _1 b(r)dr 

J o 

e usem o resultado para mostrar que as soluções da equação original (não homogénea) são 
da forma /(x) = C 1 (x)/ 1 (x) + C 2 (x)/ 2 (x). (Daí ser a fórmula da variação das constantes.) 


Ci(x) 

C 2 (x) 


(4.23) Exercício. Considerem agora a equação 

x 2 f" ~xf'+f = 3x, /(1) = / / (1) = 0. 

Sabendo que / 2 (x) = x e / 2 (x) = x logx são soluções da equação homogénea, qual a 
matriz W(x) correspondente? Qual a solução da equação? O que mudava se mudássemos 
as condições iniciais? 


(4.24) Exercício. Resolvam 

(cosx)/" + (sinx)/ 7 = (cosx) 3 , /(O) = 1, f'( 0) = 4, 

sabendo que _/j(x) = 1 e / 2 (x) = sinx são soluções da equação homogénea. (De que outras 
formas conseguem resolver esta equação?) 


(4.25) Exercício. Resolvam 

+V - 4/= 6 *> /(l) = -5, / / (l) = -8, r(l) = 0 , 

x x z x- 3 

sabendo que /j(x) = x/2, / 2 (x) = x 2 , e / 3 (x) = x 3 /2 são soluções da equação homogénea. 
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Episódio 4: Equações diferenciais de ordem superior 


(4.26) Respostas de exercícios. 


r i 

/ 

' 0 1 0 


■ / - 

f 

= 

0 0 1 


/' 



_ —c — b —a _ 




B 

+ clX‘ ^ + bX + c 


2 e £ +^ 


B [ 

mesma discutida em Q4.12 ) 


a equaçao e a 


B (a) D - 4. (b) D - 4. (c) D — 3. 

(d) (D-4)(D-3). (e) D-i. (f) D + i. 
(g)(D-i)(D + i). (h)(D-i)(D + i). 

(i) (D — 3 — 2i) 2 . 0) (D - 3 + 2i) 2 . 

(k) (D-3- 2 if (D-3 + 2 if. 

(l) (J? — 3 — 2i) 2 (D - 3 + 2i) 2 . (m) (D - 3) 3 . 

(n) (D — 3) 2 (D — 3 — 2i) 3 (D - 3 + 2i) 3 . 

(o) (D - 2i) 4 (D + 2i) 4 . (p) D (D - 4í) (D + 4i). 


( 14.11D qualquer múltiplo de (D — 2) 2 D (D 2 + 
2 2 ) 3 (D - l) 3 ((D - l) 2 + 2 2 ) ((D - l) 2 + 3 2 ) 4 . 


( |4. 16D (a) te 5£ + e 3t . (b) 2e S£ — e 4£ — 2te 4t . 
(c) e 5t — e 4t sin t. (d) e 4£ cos(2t). 

(e) e 4£ cos(2t) — ^te 4£ sin(2t). 

(f) t 3 e £ -2t 2 e £ + 2te £ -e £ . 


( R~T7l) (a)-2+(2+ ^ 5 )e l/2£ +(2-^)e-' /2£ . 
(b) te~ 2t — 2e t cos(-/3 t). 

H4.19D (a)-|e £ sin(3t). 

(b) 

—e t_í£ cos(3t) + |e £_7£ sin(3t) — |e £ sin(3t) — 1. 

{ — |e £ sin(3t), se t < ti, 

—e t-7r cos(3t)+|e í-7r sin(3t)—|e £ sin(3t)—1, sen<t. 


B3> (a)/(t)={ 


te 5£ +e 3£ , 

|e 5£ +(l— |e 2 )e 3£ . 


se t < 1, 
se t > 1. 


(b)/(t)={ 


|e 5t + |e 3t , 

-e 5t + te 5t +(f + |e 2 )e 3t , 


se t < 1, 
se t > 1. 


d4l4l ) (a) Cj e 3£ + C 2 e 5£ . (b) Q e 4£ + C 2 te 4£ . 

(c) C 1 e 4£ cos t + C 2 e 4£ sin t. 

(d) Q e 4£ cos(2t) + C 2 e 4£ sin(2t). 

(e) Cj e £ + C 2 te £ + C 3 t 2 e £ . 

( P~l5l ) (a) C 1 e 3£ + C 2 e 5£ -|te 3£ + e £ . 

(b) Cj e 4t + C 2 te 4£ + 3t 2 e 4£ - sin t - ^ cos t. 

(c) Cj e 4£ cos t + C 2 e 4£ sin t — te £ — ^e £ + 4e 4£ . 

(d) C 1 e 4£ cos(2t)+C 2 e 4£ sin(2t)+|- cos t— sin t. 

(e) Q e £ + C 2 te 1 + C 3 t 2 e £ + t 3 e £ . 


(±23 "'«=[; 

C 1 (x)=/(0)-|x 2 + |x 3 +ix 4 ; 

C 2 (x)=/ / (0) + x-x 2 . 

( I4.23P W(x) = [ j iog° S +r];/(x)= fx(logx) 2 . 
( |4.24P /(x) = i(sinx) 2 + 1 + 4sinx. 

( [4.25D /(x) = x 4 — 6x 2 , em x > 0. 
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EPISÓDIO 5 


TRANSFORMADA DE LAPLACE 


O nosso tópico neste episódio é a transformada de Laplace. O seu principal efeito é trans¬ 
formar derivadas em produtos por variáveis. Ou seja, transformar equações diferenciais em 
equações algébricas. Por exemplo, f" — 13f' + 42/ = 0 transforma-se em 

(s 2 F(s) - 5/(0) - /'(O)) - 13 (s F(s) - / (0)) + 42,F(s) = 0 

(é tradicional chamar F(s) à transformada de /(t)). Daqui é fácil obter F(s), o que signi¬ 
fica que, se tivermos uma forma prática de inverter a transformada de Laplace, temos um 
método alternativo de resolução de equações. 


(5.1) Funções definidas por troços. Como vamos lidar com funções definidas por troços, 
começamos por introduzir as funções de Heaviside. Definimos 


H 


a 


: t *-> 


1, se t > a, 
0, se t < a. 


(Notem que não especificamos Ff a (a), porque não tem qualquer papel no que vamos fazer.) 
Esta função permite-nos “cortar” o gráfico de uma função / e ficar apenas com a parte do 
gráfico à direita de a, pois 


(H a /)(t) = 



se t > a, 
se t < a. 


(Mais uma vez, não especificamos o valor da função na fronteira entre troços.) 

A mesma notação pode ser generalizada para descrever funções com mais troços, pois 
se a < b então 


(H a -H b )(t) = 


1, se a < t < b, 

0, se t < a ou b < t 
(verifiquem). Isto permite-nos escrever a função 


g : t< 


e , se t < —1, 

cos t + 1, se — 1 < t < 2, 


„3t 


e 3t + e 5t . 


se 2 < t < 3, 
se 3 < t 


como 


e f (1 - H_ x ) + (cos t + 1) (H_j - H 2 ) + e 3t (fí 2 - H 3 ) + (e 3t + e 5t )H 3 

= e f (l — H_ 1 ) + (cost + 1)(FÍ_ 1 -H 2 ) + e 3t H 2 + e 5t H 3 . 
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Episódio 5: Transformada de Laplace 


Reparem que (1 — H_ 1 ) corresponde a manter o gráfico da função desde —oo até —1. Repa¬ 
rem também que na segunda versão e 3t se mantém de 2 em diante (tanto no troço 2 < t < 3 
como no troço 3 < t) e que e 5t se mantém de 3 em diante. 

Da mesma forma, podemos reconstituir uma função por troços a partir da expressão 
com funções de Heaviside. Por exemplo, se 

Kt) = e 3t H 2 (t) + e 4t H s (0 + (e 5t - e 3t )H 6 (t), 


vemos que os pontos que separam os diferentes troços são 2, 5, e 6. A primeira parcela só é 
tida em conta à direita de 2, a segunda parcela só é tida em conta à direita de 5, e a terceira 
parcela só é tida em conta à direita de 6. Ficamos então com 


Kt)={ 


o, 


, 3 1 


e 3t + e 4t . 


,3t 


+ e 4t + (e 


51 


se t < 2, 
se 2 < t < 5, 
se 5 < t < 6, 
e 3t ), se 6 < t. 


(5.2) Exercício. Usem funções de Heaviside para escrever estas funções: 


(a) 

(d) 


í et ’ 

se t < 1, 

(b) • 


se t < 1, 


e 3t , se —2 < t < 1, 

[e t + e 2t 

, se 1 < t; 


se 1 < t; 

(c) ■ 

e 2t , se 1 < t; 

(t 2 

se 0 < t < 3, 


11 2 

se 0 < t < 1, 


t — t 2 , se 1 < t < 2, 

COS t, 

se 3 < t < 5, 

(e) < 

e f , 

se 1 < t < e. 

CO - 

t 2 + e l , se 2 < t < 4, 

^sin t, 

se 5 < t; 



se e < t; 


1 — e f , se 4 < t < 6. 


Escrevam estas funções na forma tradicional (descrição por troços) 


(5.3) Exercício. 

(a) H n + H 2 -3H 5 ; 

(c) H n + H 2 - 3H 5 + 5 H 1 - 2; 
(e) H 0 (t) + cos(3t)H_!(t); 


(b) H n + H 2 -3H 5 + 5H 1 ; 

(d) 

(0 H 0 (t) + cos(3t)H_ 1 (t) + e c . 


(5.4) Exercício. Calculem estes integrais: 


"4 


" 4 

(a) H 2 (t)dt; 

Jo 

f 4 

(b) 

H 2 (t)tdt; 

Jo 

r 4 

(c) H_ 2 (t)dt; 

Jo 

(d) 

H_ 2 (t) + H 2 (t)t+H 5 (t)dt 
Jo 


(5.5) A transformada de Laplace. Vamos considerar apenas funções / cujo primeiro troço 
é nulo (ou seja, todas as parcelas envolvem uma função de Heaviside) e satisfazendo uma 
condição do género |/(t)| < Ce at (onde a,Cel são fixos para cada função, mas podem ser 
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diferentes para funções diferentes). Definimos a transformada de Laplace de uma tal função 
por 


xifKs) = 




onde s é uma variável complexa. Frequentemente vamos usar a abreviatura 


r oo 


X c {fKs) = X{H c fKs) = 


e- st H c (t)/(t)dt = 


e~ st f(t)dt. 


Mas notem que esta abreviatura não é de modo nenhum “oficial”. Vamos usá-la só porque 
simplifica as coisas, mas não contem encontrá-la com o mesmo significado noutros sítios. 
O mais habitual é usar apenas funções / definidas em ]0, +oo[, e portanto a transformada 
de Laplace que encontram habitualmente é aquilo a que estamos a chamar S£q {/}. Vamos 
também usar essa convenção: se falarmos em “transformada de Laplace de /” sem indicar 
um valor de c e se / não envolver funções de Heaviside (ou seja, se não houver nenhum 
valor de c implícito), então vamos assumir que se trata de c = 0. 

Se fizermos a mudança de variável (qual?), podemos observar que 


^cí/(0Ks) = 


e~ sc f(t)dt = 


-s(t+c) 


/(t + c)dt 


e- st Ht + c)dt = e- cs X 0 {f(t + c)Ks), 


ou, alternativamente, 

X c {fít-c)Ks) = e- cs X 0 {f(t)Ks) 
(conseguem traduzir de uma versão para a outra?). 


(5.6) Começamos por calcular 


n oo 


X 0 {e at } = 


r oo 


= lim 

T—>oo 


-st+at dt = 


Q—(.s—a) t 


0 


e" (M)t dt= lim 

T —>oo 

T g—(s—a)T 

= — lim - 

t= o T—>oo s — a 


-( s~a)t 


dt 


Notem que a primitiva (feita no início da segunda linha) só é válida ses^ae que o limite 
final só existe se Res > a. (Conseguem verificar que o integral original, mesmo com a 
primitiva correta, também não converge ses / a?) Temos também 

g—(s—a)c 

S£ c {e at j = e - cs ^ 0 {e a(t+c) } = e ~ cs e ac ^ 0 {e at } = 


s — a 

„at 


Algo semelhante se passa com funções / satisfazendo |/(t)| < Ce at (ou seja, funções / 
dominadas por alguma exponencial). Usando s 0 = Res: 




7(t)dt 


< 


1 


1/ (01 dt 


< C 


e~s 0 1 e at át = c ^f c {e at }(s 0 ) = C 


-(s 0 -a)c 


s 0 -a 


62/99 



Copyright © 2015, João Pedro Boavida. Versão atual disponível em http : //web . técnico. ulisboa.pt/joao .boavida/2015/ACED/. 


Episódio 5: Transformada de Laplace 


o que mostra que o integral impróprio definindo E£ a {f} está bem definido pelo menos no 
semiplano direito Res > a. (Veem porquê?) 

Embora o integral só esteja bem definido num semiplano direito, para as funções / que 
vamos considerar, a função E£ a {f } pode ser estendida a todo o plano complexo (exceto num 
número finito de polos); todos os polos ficam fora desse semiplano direito. 


(5.7) Exercício. Determinem as transformadas de Laplace das seguintes funções: 
(a) e 5t ; °- 5t - r^zir^St. ^ ~ 2 it 

(e) 10e 4t + e 5t ; (f) 3e 


(b) 3e 5t ; 

(c) H 2 (t)e 5t ; 

(d) e 

(f) 3e 4t - e 6t ; 

Cg) H 0 (t)e 4í -H 2 (t)e 5t ; 

(h) J 


e 4t ,se 0 < t < 2, 


e 5t ,se 2 < t. 


Recordem a convenção que estamos a usar: se não houver indicação em contrário, assumi¬ 
mos que as funções são não-nulas apenas em ]0,+oo[ (e neste caso, E£{f} = ££ 0 {f}). 


(5.8) Exercício. Quais as funções com estas transformadas de Laplace? 



(e) 


e -2( s -D 

s-1 ’ 


20 

(b) —r; 
s - 3 

(c) 

e ~ 2s e 2 


CD —- + 

s — 1 s 

(g) 


20 4 

s - 3 s - 5’ 
3 

(s-l)(s + 2) ; 


(d) 


e 


-2 s 


s 


(h) 


e~ s +1 
1*^1 ' 


(5.9) Conhecendo a transformada de Laplace da exponencial, podemos facilmente obter a 
transformada de Laplace dos senos e cossenos, usando 

Y 0 {e at e ibt } = & 0 {e ía+bi)t } =---. 

s — a — bi 


Ficamos então com 

X 0 {e at cosíbt)} = xJe at elht+ 2 e lbt | = lif 0 { e (° + bõt} + l i?o { e (a-l , Ot} 

1 1 s —a + bi+s —a —bi 

2(s — a — bi) 2(s — a + bi) 2 (s — a — bi)(s — a + bi) 

0 ~ a) 

(s — a) 2 + b 2 

Da mesma forma, 

J? 0 {e at sin(bt)} = ^e at £lbt ~ C ^ j = ^ 0 {e (a+60t } ~ ^ 0 {e (fl ~ b0t } 

1 1 (s — a + bi) — (s — a — bi) 

2i(s — a — bi) 2i(s — a + bi) 2i(s — a — bi) (s — a + bi) 

b 

(s — a) 2 + b 2 

(Podemos usar o método do final de Q5.5D para calcular S£ c {e at cos(bt)} e J£ c {e at sin(bt)}. 
Mas como podem constatar rapidamente, o resultado não é nada de especial...) 
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(5.10) Exercício. Quais as transformadas de Laplace destas funções? Recordem que se 
não houver indicação em contrário, assumimos que as funções são não-nulas apenas em 
]0,+oo[ (e neste caso, S£{f} = i? 0 {/}). 

(a) e 3t sin(5t); (b) e 3t (cos(5t)) 2 ; (c) H n (t)e 3t sin(5t); (d) H 7r (t)e 3t (sin(5t)) 2 . 


(5.11) Exercício. Quais as transformadas de Laplace das funções do exercício ( |5.3D ? 


(5.12) Exercício. Quais as funções cujas transformadas são as seguintes? Usem funções de 
Heaviside ou descrevam os troços individualmente. (Pode ser útil rever o final de ( |5.5D .) 


(a) 

(d) 


s — 3 

s 2 — 6s + 25 ’ 
s - 4 

s 2 - 6s + 25 ’ 


(g) 


3 

(s 2 — 2s + 2) (s 2 — 2s + 5) ’ 


(b) 

(e) 

(h) 


3 

s 2 — 6s + 25 ’ 

1 

s 2 — 2s + 2 ’ 

1 

~õ -1" 

s — 2s 2 


_2s 



(c) 

(0 

(i) 


s 

s 2 — 6s + 25 ’ 

(s 2 — 2s 2) — s 
s (s 2 — 2s + 2) 


s 2 — 2s 2 


(5.13) Continuamos a nossa exploração das propriedades da transformada de Laplace, com 

r 00 rT 


x c {fXt)Ks) = 


e ~ st f'(t) dt = lim 

T —>oo 


e~ st f\t) dt. 


Usando integração por partes, obtemos 
rT 


e~ f (t)dt = 


-i t 


e_s /(t) 


rT 


+ 


_ t=C 


se~ st f(t)dt, 


o que nos leva (como?) a 


X c {f , (t)Ks) = -e- cs fíc)+sX c {fKsl 


Já agora, é razoável perguntar qual a derivada (em ordem a s) de S£ a {f }(s). Assumindo 
que podemos derivar a função integranda em ordem a s, ficamos com 


5 ^ ,= á 


e /(t)dt = 


de~ sc 

ds 


■/(t)dt = 


r oo 

— te~ st f (t)dt = —if c {t/(t)}. 
J C 


Porém, é preciso algum cuidado: não só estamos a assumir que podemos trocar a derivada 
em ordem a s com o integral em ordem a t, mas o integral em ordem até um integral 
impróprio, pelo que estamos na verdade a trocar a derivada do limite de um integral pelo 
limite de um integral da derivada. 

Prossigamos então com mais cuidado. Temos (verifiquem que os passos são seguros) 


= lim 

ds h-> o 

= lim 
h^O 


-^cí/Ks + h)--^ {/K 5 ) 


= lim lim 

h—>0 T —>oo 


fVc^/Mdt- f T e~ st f (t) dt 


lim 

T —>oo 


r T g -{s+h)t _ e ~st 


-/(t)dt = lim lim 

/i^>OT->oo 


r T e~ ht — 1 


-e- s 7(t)dt. 
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Episódio 5: Transformada de Laplace 


Para podermos continuar, precisávamos que o limite com h —* 0 estivesse dentro do integral. 
Mas como 


lim 


e~ ht - 1 


= -t 


h-> o h 

(trata-se da derivada de e~ ht em ordem a h ), sabemos que para qualquer e > 0 existe 5 > 0 
tal que 

e~ ht - 1 


\h\<5 


+ t 


< e. 


Recordem agora que estamos a assumir que |/(t)| < Ce at para constantes a,C el. Então, 
para |h| < 5, temos 


lim 

T —>oo 


r T e~ ht - 1 


7 -e~ f(t)dt - lim (-t)e“ s 7(t)dt 

h T^oo I 


T 


lim 

T —>oo 

n 

J 

c 


+ t e~ st f (t)dt 


< lim 

T —>oo 

C T 

J 

c 


e~ ht - 1 


+ t 


e ~ st |/(0| dt 


rT 


< Ce lim 

T —>oo 


e~ st e at dt = CeS£ c {e at } = Ce- 


~(s—á)c 


s — a 


Daqui concluímos (porquê?) que 


lim lim 

h ->0 r->oo 


r T 


+ t e st f(t)dt = 0, 


o que nos leva (como?) a 


lim lim 

íi—>0 T—>oo 


r T e~ ht - 1 


rT 


, -e _ /(t)dt = lim 

h T-> oo 


-te- s 7(0dt. 


que é o que precisávamos para continuar o nosso cálculo. Voltando a ele: 

rT -ht i 


= lim lim 

ds h—>OT->oo 


rT 


= lim 

T —>oo 


e“ s 7(t)dt 


te s 7(t)dt = -i? c {t/(t)}(s). 


(5.14) Exercício. Usem este resultado para calcular S£ c {t k e at }. 

(5.15) Resumindo tudo o que dissemos: 




ds 


S£ 0 {e at cos(bt)} = 


s — a 


(s - a) 2 + b 2 ’ 


X c {f'}=sl£ c {f}-e- cs f{c), 
k\ 


X 0 {t k e at } = 


C£ Q {e at sin(bt)} = 


(s — a) fc+1 ’ 
b 

(s — a ) 2 + b 2 
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(5.16) Exercício. Quais as transformadas de Laplace das funções do exercício fl4.9D? 


(5.17) Exercício. Qual a transformada inversa de cada uma destas funções? 


(a) 


(d) 


(g) 


s 2 + l 


2s 


(s 2 + l) 2 ’ 
s — 1 

(s-l) 6; 


(b) 

(e) 

(h) 


s 2 -l 


s — 1 


(s 2 — 2s + 10) 2 ’ 
s 


(s-1) 6 ’ 


(c) 


(f) 


(i) 


—2s 


(s 2 — l) 2 ’ 

1 

(5-i) e; 

„3s 


(s-1) 


6 ' 


(5.18) Resolução de EDOs usando transformada de Laplace. Mas o que tem a transfor¬ 
mada de Laplace a ver com equações diferenciais? Por exemplo, conseguimos resolver a 
equação 

f' = 3f, /(0) = 12 

usando a transformada de Laplace? Sim. 

É tradicional usar maiúsculas, como F, para a transformada de Laplace de funções em 
minúsculas, como /. Ou seja, é comum escrever F(s) = i? 0 {/(t)}(s) (e aqui vamos sempre 
usar a transformada com c = 0). Aplicando a transformada a ambos os membros da equação 
diferencial, ficamos com 

&{/'} = &{3f}. 

Mas 


^íf'} = sE£{f} — f (0) = s F(s) — 12 e S£{3f] = 3S£{f] = 3F(s), 
o que nos leva à equação 

sF(s)- 12 = 3F(s). 

Reparem que esta equação já não envolve derivadas: transformámos uma equação diferen¬ 
cial numa equação algébrica. De facto, é fácil resolvê-la: 

12 

sF(s) -3F(s) = 12 <=> F(s) =-. 

s - 3 

Nem sequer precisamos de nos preocupar com o caso s = 3: estamos já no mundo das 
funções meromorfas, e estamos apenas a constatar que F tem um polo em s = 3. 

Temos assim 

X 0 {f}=X{H 0 f} = --. 

s - 3 

Assumindo que a transformada de Laplace é injetiva, concluímos que 

H 0 (t)/(t) = 12H 0 (t)e 3t , ou / (t) = 12e 3t . 

Este raciocínio é potencialmente subtil. A solução / que procuramos é uma função con¬ 
tínua com derivada contínua (porquê?). Mas a transformada de Laplace só opera sobre 
funções com o primeiro troço nulo. Mais: como temos a condição inicial em t = 0, e como 
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S£ c {f'} = J£{H c f'} envolve /(c), convém-nos cortar a função em t = 0. Assim, o que real¬ 
mente usamos é {/} = ££{H 0 f} e = S£{H Q f'}. Assim, no final quem obtemos 

(invertendo a transformada de Laplace) é H 0 f. 


(5.19) Exercício. Resolvam estas equações diferenciais: 

(a) f - 3/ = 2e f , /(0) = 0; 

(b) f' - 3/ = 2e c , /(0) = 1; 

(c) f'-3f = e 3t , /(O) = 2; 

(d) f'~3f = —3, /(0) = 3; 

(e) f' — 3f = e 3t — 2e t , /(O) = 4. 


(5.20) Vejamos outro exemplo. Vamos resolver outra vez a equação ( |4.12D 
f"+f'-2f = 6e t + 2, f (0) = / / (0) = 3. 

Temos 

X{f'}=sF{s)~f{0)=sF{s)-3, 

F£{f"\ = s -/'(O) = s (s F(s) - 3) - 3 = s 2 F(s) - 3s - 3, 


portanto a equação original passa a 


(s 2 F(s) - 3s - 3) + (s F(s) - 3) -2 F(s) = 


s — 1 


2 

s 


S£{f"} 

ou, depois de simplificar, 


S£{f} 




(s 2 + s — 2s) F(s) — 3s + 6 + 


6 


^{6e‘} S£{2\ 


+ -. 


S — 1 S 

Resolvendo em ordem a F(s) e fatorizando o denominador, chegamos a 

6 , 2 


3(s + 2) + —- + = 3 

F(s) = - 5 =- r + 


6 


+ 


(s + 2)(s-l) s-1 (s-l) 2 (s + 2) s (s + 2)(s — 1) 

Para continuar, precisamos de decompor as duas últimas frações em frações parciais. Temos 

6 2 _ 2(s + 2)-2(s-1) 2s — 2(s — 1) 

(s — l) 2 (s + 2) s (s + 2)(s — 1) (s — 1 ) 2 (s + 2) s(s + 2)(s —1) 

2 2 2 2 
+ ■ 


(5 — 1) (s — 1) (5 + 2) (5 — 1) (5 + 2) 5 (s + 2) 


(s + 2) — s 


1 1 
- - + 


Assim, chegamos a 
F(s) = 


e concluímos que 


s-1 (s-1) 2 s s + 2 


(s-1) 2 s(s + 2) (s-1) 2 s s + 2 


—I-— g{3e t + 2te f — 1 + e 2t } 


/(t) = 3e c + 2te f — 1 + e 2í . 
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(5.21) Exercício. Escolham condições iniciais em t = 0 para as equações de fl4.14D e fl4.15D . 
Usem a transformada de Laplace para resolver essas equações usando essas condições iniciais. 
Que adaptações é preciso fazer se não conhecermos as condições iniciais? 


(5.22) A transformada inversa. Por esta altura, já devem ter reparado que, mesmo já 
tendo F(s), nem sempre é imediato determinar a função /(s) original. Por vezes, é mais 
rápido usar a transformada inversa. 

Se o integral que define F(s) = £t? c {f}(s) = i?{H c /}(s) (onde / é alguma função domi¬ 
nada por e at ) convergir em Res > a, então a transformada inversa é 


r b+iT 


^~ L {FKt) = (H c fXt)=— lim 

2m T—><x 


e st F(s) ds, 


b-iT 


onde b é um real em Res > a e o integral é feito ao longo de um segmento de reta. (Não 
vamos tentar provar este resultado.) É prático escrever este integral com uma notação 
ligeiramente diferente: 

^ r b+ioo 

(H c /)(t) = — e 5t F(s)ds. 

2m Jb-ioo 

Reparem que a expressão de S^~ 1 {F} não depende do ponto c onde começa o pri¬ 
meiro troço de H c f . Porém, se estivermos a considerar alguma função / que não está 
“cortada” para trás de nenhum c, temos de a “cortar” (ou seja, temos de escolher algum 
c—tipicamente o ponto onde temos condições iniciais—e calcular i? c {/}), e a função que 
obtemos a partir deste novo integral é a função “cortada”. 

Notem também que a exponencial neste integral (que define a transformada inversa 
i? -1 {F}) tem o sinal oposto da exponencial no integral que define a transformada direta 
(isto é, a transformada original Isto é um fenómeno habitual com transformadas 

integrais (essencialmente, porque todas são casos particulares de uma transformada mais 
geral que tem essa propriedade). 

Outro ponto importante é que a reta ao longo da qual calculamos o integral (a reta 
vertical Res = b) está contida na região Res > a (isto é, b > a), onde a é a constante da 
exponencial e at que domina/. Em particular, F não tem quaisquer singularidades à direita 
dessa reta. 

Vamos usar o teorema dos resíduos para calcular o integral da transformada inversa. 
Deste ponto em diante vamos considerar apenas funções F satisfazendo 


|F(s)| < C 


quando s —* oo, 


onde C, d e R e k é um inteiro positivo. (Como só nos interessa o comportamento quando 
s —* oo, não se justifica ajustar esta condição de modo a dar-lhe sentido em s = 0.) Notem 
que todas as transformadas com que temos lidado satisfazem esta condição. 

Vamos concluir que, se t < d e a reta Res = b estiver à direita de todos os polos de F(s), 
então 

-y ç b+ioo 

—: e st F(s)ds = 0. 
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Por outro lado, se t > d (e a reta R es = b continuar à direita de todos os polos), então 


1 

2ni 


b+ioo 


e st F(s)ds = soma dos resíduos des 


b—ioo 


e st F(s). 


(5.23) Exercício. Para o caso t < d, considerem o caminho y que contorna o retângulo de 
vértices b — iX 2 , X — iX 2 , X + iX 2 , b + iX 2 . 

(a) Usem o teorema dos resíduos para mostrar que e st F(.s) ds = 0. 

(b) Usem estimativas de integrais para mostrar que os integrais ao longo dos segmentos 
horizontais e do segmento à direita convergem para 0 quando X —* +oo. (Não se esqueçam 
que t < d, nem das hipóteses sobre F(s).) 

(c) O que podem concluir sobre /(t) = S£~ X {F } quando t < d? 




(5.24) Exercício. Para o caso t > d, considerem o caminho y que contorna o retângulo de 
vértices -X - iX 2 , b - iX 2 , b + iX 2 , -X + iX 2 . 

(a) Usem o teorema dos resíduos para calcular J e st F(s) ds. 

(b) Usem estimativas de integrais para mostrar que os integrais ao longo dos segmen¬ 
tos horizontais e do segmento à esquerda convergem para 0 quando X —* +oo. (Não se 
esqueçam que t > d, nem das hipóteses sobre F(s).) 

(c) O que podem concluir sobre /(t) = S£~ X {F } quando t > d? 

(5.25) Vejamos se isto nos ajuda com o exemplo fl5.20P . Tínhamos 

J2>{H 0 /}=F(s) = +-- 6 + 2 --, 

s_—(s — iy (s + 2) s (s + 2)(s — 1) 

S £,oUeq 
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queremos determinar /. Para a primeira parcela nem precisamos da transformada inversa, 
pois é imediato que a função original é 


^" 1 {jn}=«o(t)(3e t ). 
A função original para a segunda parcela é 


JT 


6 


(s - 1 ) 2 (s + 2) J 2ni 


r 2+ioo 


6e 


st 


2—ioo 


(s — l) 2 (s + 2) 


ds. 


Se t < 0, sabemos que o resultado é 0. Se t > 0, precisamos dos resíduos. Temos um polo 
simples em s = — 2 e um polo duplo em s = 1. Os seus resíduos são dados pela fórmula de 
Cauchy: 

6 e st í 6e st \ 6e~ 2t _ 2 _ 2t 

£% (s - l) 2 (5 + 2) - ( (TT? J s= _ 2 - (=3? “ 3 C ’ 

6 e st Í6e st \ í 6te st {s+ 2) - 6e st \ , 2 t 

Res --- 5 --- = - = ----r- = 2 te 1 -e f . 

s=t (s -1) 2 (s + 2) \^s + 2y s=1 (s + 2) 2 J s =i ^ 

Ou seja, a função original é H 0 (t)(|e _2t + 2te t — |e f ). 

Finalmente, a função original para a terceira parcela é 


í 2 1 

= J_ 

* 2+í °° 2e s t 

\s(s + 2)(s - 1) J 

2ni % 

2—ioo S +2) ( 5 


Como o integral só é não-nulo se t > 0 e como todas as singularidades são polos simples, 
ficamos com 


H 0 ■ (Res+ Res +Res) = H 0 (t) 

s=0 s=—2 s= 

= H 0 (t) 


2 e 


,St 


- (s + 2) (s 

1 

1 3 6 


1 ) J s =0 
2 

+ 3 6 ' 


+ 


2e st \ / 2e st a A 

s (s — 1) Js=-2 V5 (5 + s=i J 


Juntando as três parcelas, chegamos a 

(fío/Xt) = H 0 (t) (3e f ) +H 0 (t) (|e- 2t + 2t e f - |e f ) +H„(t) (-1 + |e“ 2t + |e f ) 
= H 0 (t) (3e f + e“ 2t + 2te t - 1). 


(5.26) Vejamos um último exemplo, desta vez com uma função definida por troços. A 
equação 

f" — 2f' + 10/ = —10H n , f (0) = 0, /'(0) = -2 

pode ser resolvida usando a transformada de Laplace. Se F = Y 0 {f}, temos 

2’ 0 {f'}=sF(s)-f(.0) = sF{s), 

&o{/") = s 2 Hs ) - 5/(0) - /'(0) = 5 2 F(s) + 2, 

— 10e~ ns 

&oi~10 H n } = -. 

s 
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Assim, ao aplicarmos a transformada de Laplace à equação original obtemos 


(s 2 F(s ) + 2) - 2(s FÇs )) + 10F(s) = 


— 10e _7IS 


s 


que nos leva a 

2 10e~ ns 

^- —= --. 

s 2 - 2s + 10 s{s 2 -2s + 10) 

Para determinar H 0 f, vamos separar F(s) em parcelas, de acordo com as exponenciais 
envolvidas (pois a exponenciais com constantes diferentes, correspondem troços diferentes 
onde a função que procuramos é não-nula). Notamos ainda que 


s 2 - 2s + 10 = (s - 1 - 3i) (s - 1 + 3i) = (s - l) 2 + 3 2 . 


A função original é então dada por 


(fí 0 /)(t) = ^ 


-i 


-2e 


-Ot 


(s — l) 2 + 3 2 


■se 


-i 


-10e” 


s (s — 1 — 3i)(s — 1 + 3i) 


A primeira parcela é não-nula apenas quando t > 0. De facto, podemos reconhecê-la 
diretamente: 


se 


-i 


(s — l ) 2 + 3 2 
A segunda parcela é 

, f — 10e _7IS 

se~ 


2 

= — S£~ 
3 




s (s — 1 — 3i)(s — 1 + 3i) 


1 

2ni 


f 2+ioo 


2—ioo 


s (s — 1 «3i)(s — 1 + 3É) 


■ ds 


(reparem que todos os polos estão à esquerda da reta vertical Res = 2, por isso podemos 
usá-la no integral). Se t < n, este integral é 0. De contrário, somamos todos os resíduos. 
Posto de outra forma, 


se 


= H 


-i 


10e" 


s (s — 1 — 3É) (s — 1 + 3i) 


= H„(t)íRes+ Res + Res ) 
s=0 s=l+3i s=l—3i 


ff -10e^~^ s \ f -lOe^-^ 5 a f -lOe^-^ 5 a > 

\ ms - 1 — 3i)(s — 1 + 3i)2s=o"*" Vs (5 — 1 + 3í)2s=i+3í"*" Vs (s — 1 — 3í) J s =i- 3 í, 


-(s - 1 — 3i)(s 
= H n (t ) 


10g( 1+3 d( t-ír ) 


: (s — 1 — 3i)J s=l 


(1 + 3í) (1+ 3i - 1 + 3i) (1 — 3i)(l — 3i — 1 — 3i) 


Nem sempre é útil, mas neste caso vamos simplificar um pouco mais. Vamos usar u = t — n 
e 10 = (1 — 3i) (1 + 3i). Ficamos então com 


lOgD+SiKt-") lo e (l-3i)(t-7r) 

(l + 3i)(l + 3i —l + 3i) _ (l-3i)(l-3Í-1-3Í) 

(1 — 3Í) g( 1+3í )u (1 + 3i)e^~ 3i ^ u g(l+3i)u _g(l—3i)u g(l+3i)u _|_ g(l—3i)u 

-- — --- — - = -l-;- + - 

6 i —6 i 6 i 2 

1 1 t 

= — 1 — -e u sin(3u) + e u cos(3u) = — 1 — -e f n sin(3t — 37 t) + e c n cos(3t — 3n) 
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Juntando tudo, obtemos (finalmente) 

(H 0 /)(t) = - |fí 0 (t) e z sin(31) + H n {t) (-1 - |e t_7T sin(3t -3n) + e l ~ n cos(3t - 3ít)) 

= - |H 0 (t) e c sin(31) + H n (t) (-1 + |e t_7T sin(3t) - e t_Tt cos(3t)). 

(5.27) Exercício. Escrevam esta função H 0 f sem usar funções de Heaviside (isto é, na forma 
habitual, por troços). Façam o mesmo com /. Comparem este método com o que usámos 
no exercício fl4.19D . 


(5.28) Exercício. Conseguem ver a relação entre a mudança de variável que fizemos ainda 
agora (iz = t — rí) e a mudança de variável que fizemos no final de ( |5.5| )? A partir daí, con¬ 
seguem encontrar uma forma um pouco mais “arrumada” de obter a transformada inversa 
de uma parcela envolvendo exponenciais? Qual a relação entre o parâmetro na exponencial 
e a função de Heaviside que ocorre na resposta final? 


(5.29) Exercício. Usem a transformada de Laplace para resolver as equações de Q4.16| ) 


(5.30) Exercício. Usem a transformada de Laplace para resolver as equações de Q4.17| ). 


(5.31) Exercício. Usem a transformada de Laplace para resolver as equações de (|4.20| 
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(5.32) Respostas de exercícios. 


( f5T2l ) (a) e l + e 2t H^t). (b) e £ + (e 2t - 

(c) e 3t /í_ 2 (t) + (e 2t — e 3t )íí 1 (t). (d) 

t 2 /í 0 (t) + (cos t— t 2 )/í 3 (t) + (sin t — cos t)iTc(t). 

(e) t 2 /í 0 (t) + (e t - t 2 )íí 1 Ct) + (t e - e f )iT e (t). 

(f) (t - + (2t 2 - t + eOHaCt) + (1- 

2e t -t 2 )/í 4 (t) + (e t -l)/í 6 (t). 


( 1531 ) (a) 


0, se t < 2, 

1, se 2 < t < 7T, /'in 

2, se 7T < t < 5, 

—1, se 5 < t. 


(C) 

(e) 

(f) 


/• —2, se t < 1, 

I 3, se 1 < t < 2, 

4, se 2 < t < n, (d) 

! 5, se n < t < 5, 

^ 2, se 5 < t. 
f 0, se t < —1, 

cos(3 1 ), se —1 < t < 0, 

[ l+cos(3t), se 0 < t. 
e 1 , se t < —1, 

e^cosCSt), se— 1 < t < 0. 
e f +l+cos(3t), se 0 < t. 



• 0, se t < 1, 

5, se 1 < t < 2, 

6, se 2 < t < ti, 

7, se n < t < 5, 

• 4, se 5 < t. 
se t < 1, 

se 1 < t < 2, 
se 2 < t < 3, 
se 3 < t. 


( 133 ) (a) 2. (b) 6. (c) 4. (d) 10. 


B («0 à- Cb) á. M T?- W á- 

z' N 10 , 1 3 1 r N 1 e - 2 ( S - 5 ) 

( e ) ~ (0 — jZç,- Cg) — 7^s~' 


(h) 


l-e-*-V , 
5-4 


g -2(s-5) 
5-5 ‘ 


( [5781) (a) e 3 t H 0 (t). (b) 20 e 3 t Tf 0 (t). 

(c) (20e 3t — 4e St )íí 0 (t). (d) H 2 (t). (e) e f íí 2 (t). 
(f) e t_ 2 H 2 (t) + e 2 íí 0 (t). (g) (e £ — e _ 2 t )H 0 (t). 
(h) sí^^H l( t)+^H 0 (t). 


d53Õl ) (a) — 


(c) 


_ Se -^( S -3) 

(s-3) 2 +5 2 • 


(s-3) 2 +5- 


■ (b) 


1 


■ + 


5—3 


2(s—3) 2((s—3) 2 +10 2 ) * 


(d) 


V.s-3 


s—3 


(s—3) 2 +10 2 ) • 


d3TTT ) (a) 


:s +e —3e~ 


(b) 


e - ra+e -z s _ 3f .-b S+ 5 e - 


. e +e — 3e +5e — 2 


(c) usando íí 0 : 

(d) É ^^ + ^; (e) i + 

(f) usando H 0 : -d- + - + 


+ 


2(5—31) 2(5+31) ‘ 


s 2 +3 2 ' 


( 1313 (a) e 3t cos(4t)íí 0 (t). 

(b) |e 3t sinÇ4t)TÍ 0 Ct). 

(c) e 3t (cos(4t) + | sin(4t)) íí 0 (t). 

(d) e 3t (cos(4t) - i sin(4t)) Tf 0 (t). 


(e) e £ sin t Tf 0 (t). (f) (1 — e £ sin t) iT 0 (t). 

(g) (e £ sin t — e £ sin( 2 t))íí 0 (t). 

(h) e l sin t H 0 (t) + H_ 2 (t). 

(i) e t+2 sin(t + 2 )H_ 2 (t). 


2(s-i) 

i 


2(s+l) 


(S-3-2Í) 2 


O 


m Ca) i{HÍ(c) J;.(d) J; + i. 
(e) —. (f) — r (g) j + 2(s+i) • 

(h)^ + ^. (i) 

(k) 

® 2(s—3—2i) 2 571 15ZTÕ2 ■ C m ) 

(n) + — 


(s-3+21) 2 ‘ 


2i(s—3—2i) 2 2i(s—3—2i) 2 ’ 


2(s—3—2i) 2 

1 


(5-3) 2 (5-3) 2 • 


(s—3) 2 í(s—3—2i) 3 

( O ) + ttZz + 


i(s—3+20 


s 2 +2 2 (s—2i) 4 (s+2i) 4 


Cp) 5: 


2s 2(s 2 +4 2 ) * 


( 13171 ) (a) cos t H 0 (t). (b) 

(c) is^ísl H 0 (t). (d) tsint H 0 (t). 

(e) Pfe t sin(3t)íí 0 (t). (f) ^tVH 0 (0- 
(g) itVffo(t). (h) (^ 4 + it s )e*H 0 (t). 

(i) £(t + 3) 5 e t+ 3 H_ 3 (tj. 

Q5.19D (a) e 3t — e c . (b) 2e 3£ — e £ . 

(c) 2e 3t + te 3t . (d) 1 + 2e 3t . (e) e £ + 3e 3t + te 3t . 


([327]) 

ru; — se t < o, 

J — |e f sin(3t), se 0 < t < n, 

^ — |e f sin(3t)—l+|e í-7r sin(3t)—e í-7r cos(3t), se n < t. 

{ — |e t sin(3t), se t < n, 

— |e f sin(3t)—l+|e í-ír sin(3t)—e t-7C cos(3t), se n < t. 

( [323 (a) te 5t + e 3t . (b) 2e St - e 4t - 2te 4t . 
(c) e 5t — e 4t sin t. (d) e 4t cos(2t). 

(e) e 4t cos(2t)— ^te 4t sin(2t). 

(f) t 3 e c —2t 2 e t + 2te t — e £ . 


( 13301 ) (a) — 2 +( 2 +d f )e l/ 2 t +( 2 —^)e-' /2t . 
(b) te~ 2t — 2e £ cos(-/3 t). 


d33[l ) (a) { 

(b) {> 


te +e , 
e 3t + 3 g 5t_ l e 3t+2 


se t < 1, 
se t > 1. 


-e^ + ^e^+te 3 ' 


se t < 1, 
, se t > 1. 
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EPISÓDIO 6 


AS SÉRIES DE FOURIER 


Imaginem que queremos resolver a equação 

du _ d 2 u 
d t dx 2 

Como é que alguém poderia resolver uma equação desse género? Na falta de outras ideias, 
uma opção é tentar encontrar soluções por tentativa e erro. Por exemplo, podemos subs¬ 
tituir iz(t,x) = T(t) -X(x) e ver o que acontece. Chegamos (como?) a T'(t ) -X(x) = 
T(t) -X"(x). Separando tudo o que tenha t para a esquerda e tudo o que tenha x para a 
direita, obtemos 

T\t) _ X"(x) 
nÕ ~~ X(x) ' 

Se essas duas expressões são iguais, então não podem depender nem de t (porque não há t 
no lado direito) nem de x (porque não há x no lado esquerdo). Ou seja, são iguais a uma 
constante A. Por exemplo, a segunda equação fica equivalente a 

(D 2 - X)X = 0 

(usando D para indicar a derivada em ordem a x). Suponhamos ainda que 0 < x < 1 e que 
iz(t, 0 ) = u(t, 1 ) = 0 (a condição de fronteira) , o que, se T(t) ^ 0 , implica X( 0 ) =X( 1 ) = 0 . 
Para termos X(x) 7 ^ 0, precisamos de A < 0 (isto não é bem assim com outras condições de 
fronteira), e nesse caso as soluções desta última equação diferencial são 

X(x) = C 1 cos(^x) + C 2 sin^yf— Ãx). 

Mas para termos X(0) =X(1) = 0, temos necessariamente Q = 0 e V—A = nk (para algum 
inteiro fc). Ou seja, neste momento já sabemos que u(t,x) = T(t) sin(Tifcx) (absorvemos a 
constante C 2 na definição de T(t)), que T'(t) = AT(t), e que A = —n 2 k 2 . 

Usando agora D para representar a derivada em ordem a t, vemos que (D — A) T = 0, 
cuja solução é T(t) = e At T(0). Ou seja, a solução final é um múltiplo de e~ K k c sin(Trfcx). 
Mais geralmente, seria razoável assumir que a solução geral seja da forma 

u(t, x) = ^ b k e~ n k 1 sin(nkx). 
k> 0 

Este é o nosso primeiro exemplo de uma série de Fourier. 

Neste episódio, falamos um pouco das séries de Fourier e das suas variantes (as séries de 
senos e as séries de cossenos). No próximo, vamos ver como usá-las para resolver algumas 
equações diferenciais parciais (EDPs)—assim chamadas por envolverem derivadas parciais. 
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(6.1) Funções periódicas. Nas próximas páginas vamos trabalhar apenas com funções pe¬ 
riódicas seccionalmente C 00 . Mais concretamente, vamos fixar um intervalo [A,B] de com¬ 
primento l = B — A, e considerar apenas funções / : M —» C seccionalmente C 00 e com 
período l. 

Dizer que / tem período t é dizer que /(x + () = /(x) para todo otel. Reparem que 
se a função tem período i, também tem período 21 ou 3£ ou —l. Ou seja, da forma como 
estamos a definir ter período l, não há nenhuma sugestão de que l seja o número mais 
pequeno para o qual /(x + 1 ) = /(x). (Conseguem imaginar alguma função com período 
l para todo o l > 0? Conseguem imaginar duas funções com período l para todo oíeQ? 
Linearmente independentes?) 

E o que queremos dizer com uma função ser seccionalmente C 00 ? No nosso caso, apenas 
que existe um número finito de pontos intermédios (que podem ser diferentes de função 
para função) com A = x 0 < < x 2 < ■ ■ ■ < x íc _ 1 < x k = B e que em cada subintervalo 

]x;,X; +1 [ a função é a restrição de uma função C 00 definida num intervalo ligeiramente 
maior. (Pôr isto desta forma evita algumas complicações com limites laterais das derivadas 
da função.) 

Alternativamente, podemos começar com funções seccionalmente C 00 definidas no in¬ 
tervalo ]A,B[ (porquê aberto?) e estendê-las aM de forma a terem período í = B — A. De 
facto, o intervalo exato que escolhemos não é especialmente relevante (porquê?). Porém, 
nos problemas com que vamos lidar pode ser mais prático lidar com intervalos de certos 
formatos. 

(6.2) Produto interno de funções periódicas. Concentramo-nos então nas funções sec¬ 
cionalmente C°° de período l. Escolhemos qualquer intervalo [A,B] de comprimento í = 
B — A, e definimos o seguinte produto interno: 


(/,g) = / 0 )g(x)dx. 

JA 

Como habitualmente, a norma correspondente é definida por ||/|| = 

(6.3) Exercício. Provem as seguintes propriedades deste produto (/, g, h são funções secci¬ 
onalmente C 00 definidas em [A,B], c é um complexo): 

(a) (f,g) = ( g,f) 

(b) (cf,g) =c(f,g) e (f,cg)=c{f,g)] 

(c) (f + g,h) = (f,h ) + (g,h) e (f,g + h) = (f,g ) + (f,h); 

(d) 11/1| > 0; 

(e) se 11 / 1 | = 0 , então / = 0 . 


(6.4) Exercício. Mostrem que se / e g têm período t e ambos os intervalos [A, B] e [A',B'] 
têm comprimento t, então 


r B rB' 

ZO)g(x)dx= /(x)g(x)dx. 


JA 


JA! 
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(6.5) Base ortogonal para funções periódicas. Vamos agora concentrar-nos em funções 
de período 2n. Não vamos tentar demonstrar tudo (ou explicar rigorosamente o que as 
diferentes afirmações significam), mas a ideia é que temos uma base ortogonal para estas 
funções. 

Onde vamos buscar uma tal base? Se o operador D (derivação em ordem a x) fosse 
diagonalizável (de novo, sem tentar definir rigorosamente as palavras que usamos em con¬ 
textos novos), então esperaríamos ter uma base de vetores próprios. Quem são os vetores 
próprios de D ? A função / é um vetor próprio de D com valor próprio À se Df = f — Xf. 
Conhecemos as soluções dessa equação: são da forma e Ax /(0). Mas são periódicas? 

(6.6) Exercício. Mostrem que /(x) = e Ax /(0) define uma função de período 2n só se 
/(O) = 0 ou X = ni, com teS. 

(6.7) Exercício. Dados n, m € Z, mostrem que as funções e mx e e mix são ortogonais (esta¬ 
mos a usar qualquer intervalo de comprimento 27t, por exemplo, [0,27t]) desde que nj^m. 
Qual a norma de e mx ? 

(6.8) Exercício. Se / e g são funções C°° de período 271 (continuamos a usar o intervalo 
[0,27t]), mostrem que: 

(a) /(271) = /(0) e g (27t) = g(0); 

(b) -j^g(x) = ^g(x) (não é a derivada complexa, pois a variável não é complexa); 

(c) {f,D 2 g) = {D 2 f ,g) (onde D representa ^); 

(d) —{f,D 2 f) > 0. 

(6.9) Exercício. Mostrem que e mx (com n e Z) é um vetor próprio de D 2 . Qual o valor 
próprio? 

(6.10) Exercício. Fixem n, m e Z com n ^ m. Tendo em conta que e mx é ortogonal a e mix , 
mostrem que cos(nx) e sin(nx) são ortogonais a cos(mx) e a sin(mx). (Continuamos a 
trabalhar no intervalo [0,2tt].) 

(6.11) Exercício. Fixem um inteiro n > 0. Usando exponenciais complexas ou outra técnica 
que prefiram, mostrem que cos(nx) é ortogonal a sin(nx). Determinem a norma de am¬ 
bas as funções. Mostrem ainda que ambas são ortogonais à função constante igual a 1, e 
determinem a norma desta última função. (Continuamos a trabalhar em [0, 2 tt] .) 


(6.12) Séries de Fourier. Em termos práticos, a conclusão de tudo isto é que vamos usar 
uma base ortogonal com as funções 1, cos(nx) e sin(nx) (usando todos os inteiros n > 0). 
Nesse caso, uma função / de período 2n pode ser escrita como combinação linear dessas 
funções: 

/(x) = a 0 + s a n cos(nx) + ^ b n sin(nx). 

n>0 n>0 
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A uma tal soma chamamos série de Fourier para /. (Notem que isto pode não ser bem uma 
igualdade; vamos esclarecer este ponto um pouco mais à frente.) 

E quem são os coeficientes da série de Fourier? Como as diferentes funções da base são 
ortogonais entre si, vamos ter: 

(/, 1) = (a 0 ,1), (/, cos(nx)) = ( a n cos(nx), cos(nx)), (/, sin(nx)) = ( b n sin(nx), sin(nx)) 

(verifiquem!), o que nos leva a (aqui, n > 0) 

_ (f, 1) _ (/, cos(nx)) _ (/, sin(nx)) 

fl ° lllll 2 ’ ün ||cos(nx)|| 2 ’ n ||sin(nx)|| 2 

(concordam?). 

No caso de um intervalo mais geral, da forma [ A,B ] e com comprimento l = B — A, 
podemos usar as mudanças de variável x = A + e u = y (x — A) = y x — para 
converter entre u e [0, 2n] e x e [A, fí ]. Nesse caso, os limites de integração devem ser 
ajustados (passar de 0 a 2n para de A a B). Podemos então usar uma base ortogonal com 
as funções 1, cos(nu) = cosfyp^ - —y—) e sin(ntí) = sintyp^ - yj^) (tal como antes, 
usando todos os inteiros n > 0). Porém, isto é desnecessariamente complicado; podemos 
perfeitamente eliminar os obtendo uma base ligeiramente diferente: 1, cos(y^) e 
sin( —). Para cada n, as novas funções são combinação linear das antigas e vice-versa 
(conseguem ver porquê?), mas a nova versão é mais prática. 

Assim, dada uma função de período í (usando o intervalo [A,B] de comprimento í = 
B — A), escrevemos 

/(x) = a 0 + 2> n c os*P+2> n sin*p 

n> 0 n>0 

e temos, como antes (e ainda com n > 0), 

(/, 1) (/, cos^) (/, sin^F) 

a °“ ||1|| 2 ’ ün ~ llcos^pil 2 ’ bn ~ llsin^pil 2 ' 

Ou, tendo em conta, como vimos em fl6.11| ), que ||1|| 2 = í e ||cos yr^H 2 = ||sin ^y^H 2 = |, 

a 0 = \{fA), a n = |(/,cos^p), b n = |(/,sin yy). 

Neste ponto, vale a pena fazer um aviso: alguns autores usam uma convenção ligeira¬ 
mente diferente, em que a 0 tem a mesma fórmula que a n (n > 0), e na série de Fourier de 
/ usam y em vez de a 0 . Com qualquer das convenções, o resultado final é 

fw= ,1) + cos—)cos — + 

n> 0 n>0 

(De novo: isto não é bem uma igualdade e vamos esclarecer os detalhes mais à frente.) 

(6.13) Vejamos um exemplo. Vamos considerar a função / de período 1 definida no inter¬ 
valo ]0,1[ por/(x) = x. (Qual é a sua definição no intervalo ] 1,2[? E no intervalo ]0,2[? 
Não dissemos qual o valor em x = 1. Faz diferença?) 
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A sua série de Fourier é dada por 


a 0 + 2 a n cos(27tnx) + ^ b n sin(27tnx), 

n >0 n >0 


e os respetivos coeficientes (de novo, n > 0) são 


r 1 


a o = (/>!) = 


x dx, 


a n = 2 (/, cos(27tnx)) = 


b n = 2 (/, sin(27tnx)) = 


r 1 
o 

r 1 


2x cos(27tnx)dx, 


2x sin(27tnx)dx. 


Basta-nos calcular estes integrais. Para o primeiro, temos simplesmente 


r 1 


x dx = 


r x 2_l1 


Para o segundo, podemos usar integração por partes: 
r 1 r_1_'v -i i r i 


2x cos(27tnx) dx = 


x sin(27tnx) 


nn 


J o 


sin(27tnx) f 1 sin(27tnx) 


■ dx = 


nn 


dx 


nn 


Este último integral é na verdade (1, dEÍllEi}) (j sto s 5 f un dona porque o integral é ao longo 
de todo o intervalo!), portanto o seu valor é 0 (porquê?) e a n = 0 quando n > 0. Para o 
terceiro coeficiente usamos de novo integração por partes: 


r i 


2x sin(27tnx) dx 


x cos(27tnx) 


-i i r 1 


nn 


J o Jo 


cos(2írnx) 


dx 


nn 


O último integral é (1, _ çosÇ^px)), q Ue também é 0 (porquê?), e ficamos apenas com a outra 
parcela: 


r i 


2x sin(27tnx) dx = 


x cos(27tnx) 


nn 


- i 

. o 


cos(27tn) 1 
nn nn 


Concluímos que a série de Fourier é 


1-2 


sin(27tnx) 


n> o 


nn 
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(6.14) Vejamos um exemplo parecido, mas agora com uma função ímpar: a função / de 
período 2 definida no intervalo ] — 1,1[ por /(x)= x. (Qual é a sua definição no intervalo 
]1,3[? E no intervalo ]0,2[?) A sua série de Fourier é 


a 0 + E a n cos(Ttnx) + ^ b n sin(7tnx) 

n>0 n>0 


e os respetivos coeficientes (com n > 0) são 




r i 


xdx. 


a n = (/,cos(7tnx)) = 
b n = (/, sin(7tnx)) = 


-i 

r i 

-i 


x cos(nnx) dx, 
x sin(7tnx) dx. 


Para o primeiro integral, temos 

r i 

-i 


x dx = 


r x 2 


= 0 . 


-i 


Também podíamos observar que a função é ímpar e o intervalo é simétrico, logo o integral 
é 0. Aliás, essa observação permite-nos concluir imediatamente que 

* i 

x cos(7inx)dx = 0. 

J-1 


Para o terceiro integral precisamos de integração por partes: 


" 1 

x sin(7tnx)dx = 

cos(7inx) 
—x- 

‘-í 

J -1 

nn 

-i J 


cos(7tnx) 


dx 


nn 


Este último integral é (1, — cos ^ nx l ) (comparem com o exemplo anterior, e reparem que 
a notação do produto interno se aplica agora sobre um intervalo diferente), ou seja, 0 
(porquê?). A outra parcela é 


~ 1 

cos(nnx) 

x sin(7tnx) dx = 

—x- 

J-i 

nn 


cos(7tn) ^cos(—7tn) 


nn 


nn 


2 tir. 

nn 
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A nossa série de Fourier fica portanto 


n>0 


(~D n 

nn 


sin(Ttnx). 


(6.15) Exercício. Seja / uma função ímpar de período 21. Mostrem que os coeficientes a 0 e 
a n (com n > 0) da sua série de Fourier obtida a partir do intervalo [—1,1] são todos nulos. 
Mostrem ainda que os coeficientes b n podem ser calculados usando 


ri 


K = j 


/O) sin^dx. 


21 


ou seja, usando a base relevante para o intervalo [—1,1], mas calculando os produtos in¬ 
ternos (tanto o integral no numerador, como o que dá origem ao denominador) como se o 
intervalo fosse [0,£], 


(6.16) Exercício. Seja / uma função par de período 2 í. Mostrem que os coeficientes b n 
(com n > 0) da sua série de Fourier obtida a partir do intervalo [—1,1] são todos nulos. 
Mostrem ainda que a 0 e a n (com n > 0) são 



/0)dx, 


ri 


/O) cos^dx. 


(6.17) Séries de senos. Por vezes é útil ter uma série exclusivamente com senos para al¬ 
guma função /. No exercício Q6. 1 5D vimos que se / for impar, a sua série de Fourier tem 
apenas senos (ou seja, os coeficientes associados aos cossenos são nulos). Isso mesmo foi 
ilustrado no exemplo ( |6.14| ). E como os senos são funções ímpares, a soma de uma série 
só com senos também é uma função ímpar. Mas às vezes precisamos de algo que parece 
contraditório: uma série só com senos mesmo quando / não é ímpar. 

Comparem os exemplos ( |6.13| ) e Q6.14| ). No primeiro caso, a função não era ímpar, 
calculámos a série de Fourier em relação ao intervalo [0,1], e a série não tinha só senos. 
No segundo caso, a função era ímpar, calculámos a série em relação ao intervalo [ — 1,1], 
e a série tinha só senos. Naturalmente obtivemos séries diferentes, mas reparem que no 
intervalo [0,1] representam a mesma função (comparem os gráficos). 

Dessa forma contornamos a aparente contradição: quando precisarmos de ter uma série 
só com senos, também só precisaremos que ela corresponda à função num intervalo da 
forma [0, i]. Portanto, podemos fazer um prolongamento ímpar da função ao intervalo 
[—1,1] e calcular a série de Fourier (para uma função de período 21) em relação ao novo 
intervalo, [—1,1]. De acordo com o exercício Q6.15D , essa série tem apenas senos; daí 
chamarmos-lhe uma série de senos. 


(6.18) Exercício. Calculem as séries de senos destas funções: 

(a) / : ]0, n[ —» K : x >-» sinx; 

(b) / : ]0, n[ —* M : x >-» sin(2x); 
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(c) / : ]0,7i [ —* R : x <-* cosx; 

(d) / : ]0,7i[ —» M : x <-> cos(2x); 

(e) / : ]0, tt[ —> M : x —> 1; 

(f) / : ]0,7i[ —* M : x >-» (cosx) 2 ; 

(g) / : ]0,2[ —> M : x —> x 2 . 


(6.19) Exercício. Considerem uma função / ímpar e de período 2Í, satisfazendo /(x + £) = 
—f{x). Mostrem que na sua série de senos em relação ao intervalo [—£,£] só aparecem 
coeficientes de índice ímpar, ou seja, b n = 0 quando n é par. (Façam uma mudança de 
variável x >-» x + í no integral e depois usem propriedades de / e do seno.) 


(6.20) Séries de cossenos. Da mesma forma, por vezes é útil ter uma série só com cossenos 
para uma função / originalmente definida em [0,£]. Como vimos no exercício Q6.161 ), uma 
tal série corresponde a uma função par. Assim, fazemos um prolongamento par de / a 
[—£,£] e calculamos a série de Fourier (para uma função de período 2£) em relação ao 
novo intervalo, [—£,£]. De acordo com esse exercício, a série tem apenas cossenos; daí 
chamarmos-lhe uma série de cossenos. (Notem que uma série de cossenos pode incluir o 
termo constante, que corresponde ao cosseno com n = 0.) 

Considerem de novo a função / do exemplo fl6.13D definida por /(x) = x no intervalo 
]0,1 [ e de período 1. Esta função não é par. Porém, se só estivermos interessados no que 
se passa no intervalo ]0,1 [, podemos usar o prolongamento par g : x —* |x| ao intervalo 
] — 1,1 [ e depois prolongar periodicamente, com período 2. Esta função é par, e portanto a 
sua série de Fourier em relação ao intervalo [ — 1,1] é uma série de cossenos. 

(6.21) Exercício. Calculem a série de Fourier de g em relação ao intervalo [ — 1,1] e confir¬ 
mem que é realmente a série de cossenos para /. 



(6.22) Exercício. Considerem agora a função de período 2 definida por /(x) = x no inter¬ 
valo ]0,2[. Determinem a série de Fourier, a série de senos, e a série de cossenos dessa 
função. Comparem com as séries dos exemplos ( |6. 13] ) e fl6.14[ ) e do exercício anterior. 


(6.23) Exercício. Calculem as séries de cossenos destas funções: 

(a) / : ]0,7t[ —» K : x —» 1; 

(b) / : ]0,7t[ —> M : x —> cosx; 

(c) / : ]0, n[ — » E : x >-» cos(2x); 
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(d) / : ]0, n\_ — * R : x —* (cosx) 2 ; 

(e) / : ]0, n[ —» M : x <-* (sinx) 2 ; 

(f) / : ]0, ti[ —> R : x >—*• cos(3x). 


(6.24) Exercício. Considerem uma função / par e de período 2 í, satisfazendo /(x + í ) = 
—/(x). Mostrem que na sua série de cossenos em relação ao intervalo [—£,£] só aparecem 
coeficientes de índice ímpar (atenção: ímpar!), ou seja, a n = 0 quando n é par. 


(6.25) Valor de uma série de Fourier num ponto. Uma pergunta natural é: a função 
obtida a partir da série de Fourier coincide com a função original / ? 

Se / for C°°, sim (mas não vamos tentar demonstrá-lo). De contrário, é preciso algumas 
precauções. Para as funções que temos usado (limitadas, periódicas, seccionalmente C 00 ) 
e nos pontos onde / é contínua, a soma da série coincide com o valor de /. O valor da 
série num ponto de descontinuidade de / é a média dos limites (à esquerda e à direita) de 
/ nesse ponto. (É tradicional exigir apenas que a função seja seccionalmente C 1 , mas nos 
exemplos que vamos explorar isso não acrescenta nada.) 

Ilustramos com a função / do exemplo fl6.14| ) : uma função de período 2, definida em 
] —1,1 [ por/(x) = x, e com série de Fourier 

,(-D n 


Z- 

n> 0 


nn 


sin(Trnx). 


Em pontos no interior do intervalo [ — 1,1], a função original é C 00 , por isso a soma 
coincide com o valor da função. Por exemplo, em x = |, ficamos com 


/(\) =| = 2 sin(Tini). 

n> 0 


nn 


Se n for par, o seno é 0. Se n for ímpar, temos n = 2k + 1 (com k > 0), e então 


f_l ~\2k+l 

|=2 ~ 2 -r^ T ^ SÍn O 

k> 0 


2fc+l 


n (2 k + 1) 


>=z 


sin(fc7i + §) = 2 


fc >0 


n (2 k + 1) 


k> o 


n (2 k + 1) 


(- 1 )" 


Desta forma, multiplicando tudo por jj, obtemos 


?=£ 


k> 0 


(-D* 

2 fc +1 


(6.26) Para perceber o que se passa num ponto de descontinuidade, voltemos ao exemplo 
( 16.131 ): uma função / de período 1 definida em ]0,1[ por /(x) = x, e com série de Fourier 


1 

2 


E 

n> 0 


sin(27rnx) 

nn 


Esta função tem descontinuidades nos inteiros (porquê? consultem o gráfico). Como a 
função é periódica, basta avaliarmos o que se passa em x = 0 (pois nos outros inteiros 
será igual). Reparem agora que o valor da série em x = 0 é a média do limite à esquerda 
lim x ^ 0 - /(x) e do limite à direita lim x ^ o +/(*)• 
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(6.27) Exercício *. Mantemos a função / do parágrafo anterior e consideramos uma fun¬ 
ção g de período 1, limitada, e seccionalmente C 00 . Suponham que lim x ^ 0 - g(x) = a e 
lim x ^ 0 + g(x) = b, com aj^b. Definam h : x >-» g(x) + (b — a)/(x) em x & TL. 

(a) Determinem lim x ^ 0 - h(x) e lim x ^ 0 + h(x), e mostrem que é possível prolongar h a 
x = 0 por continuidade. Qual o valor em x = 0 da função prolongada? 

(b) Qual a soma da série de fr(x) em x = 0? Qual a soma da série de (b — a)/(x) em 
x = 0? Qual é então o valor da soma da série de g(x) em x = 0? 

(c) Conseguem generalizar esta conclusão para outros pontos de descontinuidade (além 
de x = 0)? E para outros intervalos (além de [0,1])? 

(6.28) Exercício. Considerem a função de período 1 definida em ]0,1[ por h(x) = x 2 . 

(a) Determinem a série de Fourier de h(x). 

(b) Substituindo x = determinem o valor de Eg=i 

(c) Substituindo x = 1, determinem o valor de Eig=i 72 - 

(6.29) Exercício. Considerem agora a função de período 1 definida em ]0,1[ por fr(x) = x 4 . 

(a) Determinem a série de Fourier de fr(x). 

(b) Substituindo x = ^, determinem o valor de V 00 , ^ V . 

(c) Substituindo x = 1, determinem o valor de E^Li 74 - 


(6.30) Respostas de exercícios. 


( fC7l ) V2n. 

H6.9D -n 2 . 

H6.11P ypri\ </n; V2n. 


H 6 . 18 P (a) sinx. (b) sin(2x). 

(c) E sin(nx). (d) E ^ sin(nx). 


n>l 
n par 


n> 0 
n ímpar 


Ce) E 77 sin(nx). (f) E 5 n "_ 4 2 n \ sin(nx). 


„ .... „ n(n 3 — 4n) ' 

n>0 n>0 v ' 

n ímpar n ímpar 

16+8(7T 2 n 2 —2)C—l) n „• nnx 
n> 0 ' n 3 n 3 S1U 2 * 


(g) E 


(|6.2ip i - En >0 -ricosínnx). 

1 n ímpar n n 


d6.23|) (a) 1. (b) cosx. (c) cos(2x). 

(d) | | cos(2x). (e) | | cos(2x). (f) cos(3x). 


46271 ) (a) b. b. b. (b) b. 

46281 ) (a) 

| + E„> 0 7^7 cos( 27 inx) - E n>0 77 sin( 27 inx). 

(b) -g. (0 f 


W |+E„>o(t7ê- 7v) C0S ( 27Tnx ) + 

E n >o(7^ - 77) sin( 27 inx). (b) -^f. (c) g. 
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EPISÓDIO 7 


ALGUMAS EDPs 


Neste episódio, resolvemos alguns exemplos de equações diferenciais parciais, ou seja, equa¬ 
ções que envolvem derivadas parciais. 

Começamos com alguns exemplos de equações de primeira ordem. Imaginem que têm 
uma função u : R 2 —* R : x, y >-» u(x,y) e que y : t >-» y(t) é um caminho em R 2 . Abusando 
um pouco a notação, escrevemos y(t) = (x(t), y(t)). Omitindo t sempre que não haja risco 
de confusão, temos 


é( u(r(tM ) 


éC-M 


dx 

dt 


du 

dx 


(*,y) + 


dy 

dt 


du 

dy 


(x,y). 


Isto sugere que, se tivermos uma equação 


du du 

a(x,y)- ^-(x,y) + b(x,y)- — 0,y)=/(x,y), 
dx dy 


procuremos, usando as técnicas do episódio [3] uma solução para 


dx 

dt 




dy 

dt 


b(x,y) 


e construamos com ela um caminho y(t) = (x(t),y(t)) (uma curva característica'). Substi¬ 
tuindo esse caminho na equação, ficamos com 

3 u 3 u 

aW0j(0)'7-Wt)j(0) + í>Wt)j(0) , 7-W0j(0) = /(x(0.y(t)) 1 

'-v-' dx ' ---' dy 

dx dy 

dt dt 


ou seja, tendo em conta o que vimos acima, 

^(u(r(0)) =/0(t),y(t)). 

Esta é uma equação de primeira ordem, apenas na variável t, que podemos resolver com 
as técnicas do episódio [I] Se conhecermos o valor de u num ponto do caminho, podemos 
usá-lo como condição inicial e obter os valores de u ao longo de y, desde que a equação 
tenha uma solução única, e desde que não haja outros caminhos a cruzar-se com este e que 
indiquem valores diferentes (usamos as ideias do episódio [2] para esclarecer estes detalhes). 

Também veremos exemplos de equações de segunda ordem. Para estas, as técnicas 
do episódio [6] permitem-nos escrever a função procurada como uma série de senos ou de 
cossenos (conforme os casos), e a partir da EDP obter equações de ordem superior para os 
coeficientes da série, equações essas que resolveremos como no episódio [4] ou no [5] 
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(7.1) O método das características. Considerem a equação 

du 2 

— = ~yu, u(0,y) = y . 
dx 

(Neste e nos próximos exemplos e exercícios, vamos trabalhar com funções u(x,y).) Depen¬ 
dendo do contexto, à segunda condição chamamos condição de fronteira (quando vemos x 
como variável espacial) ou condição inicial (quando vemos x como variável temporal, ou a 
condição como condição inicial de uma EDO). 

Como a equação não envolve derivadas parciais em relação a y, podemos tratar y como 
constante e observar que nesse caso = — yu implica u = C e~ yx (ainda se lembram como 
resolver estas equações?), onde C é uma constante em relação a x mas pode depender de 
y. Como nos é dito que iz(0, y) = y 2 , substituímos x = 0 em iz(x,y) = C e~ yx e obtemos 
iz(0, y) = C. Concluímos assim que C = y 2 e que u(x, y) = C e~ yx = y 2 e~ xy . 


--- 




(o,y n ) 

(D To) ^ 





-- 



Vejamos isto à luz do que dissemos há pouco sobre as curvas características. Se escolher¬ 
mos uma variável independente t e escrevermos x(t) = t e y(t) = y 0 , onde y 0 é constante, 
ficamos com ^ = 1 e ^ = 0. Dessa forma, 

^(“ wa J(t ») = cã(“ (t>y o ) ) = 

é o lado esquerdo da equação. Se escrevermos g(t) = iz(x(t), y(t)), constatamos que 

g / (t) = ^ (“OCO, y(0)) = ^WOj(t)) = -y(0 • ri(x(t),y(t)) = —y 0 g(t). 

Isto é, g'(t) = —y 0 g(t). A condição inicial desta EDO é g(0) = u(0,y 0 ) = yg. Ora, a 
solução dessa equação com essa condição inicial é g(t) = y^ e~ yot (verifiquem). Portanto 
u(t,y 0 ) = g(t) = Jo e ~ yot - Emx = t ey = y 0 , temos u(x,y) = y 2 e~ xy . 

(7.2) Vejamos outro exemplo: 

du 

— = u-cosx + e y, u(0,y) = e y . 

5x 

Também aqui, usamos as características x(t) = t e y(t) = y 0 . 

Em vez de usar g(t) = u(x(t), y(t)) como acima, vamos simplesmente escrever u (que 
será u(x,y) ou u(x(t),y(t)) conforme os contextos) e usar e (quando estejamos a 
referir-nos às derivadas de x,y >-» u(x,y)) ou ^ (quando estejamos a referir-nos à derivada 
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de t >-» u(x(t),y(t)), ou seja, a g^t)). Isto permite-nos aliviar a notação substancialmente, 
sem trazer riscos sérios de confusão. Temos assim 

d u dx du dy du du 

dt dt dx dt dy dx 

e portanto, substituindo x = t e y = y 0 , a equação original é equivalente a 


u' = u- cos t + e sült y 0 , u(0,y 0 ) = e^, 

onde u' = u'(t) = ^ e y 0 é uma constante. Esta equação é linear e podemos resolvê-la como 
em (|1.22D. Escolhendo /i(t) = e _smt como fator integrante, obtemos consecutivamente 


e- smt u' 


cos t • e sm c u = y 0 


(e -sint u)' = y 0 


- sm t,, _ 


u = To I + C, 


onde C é uma constante (em relação a t). Desta forma, u = y 0 te smt + Ce smt . Como x(0) = 
0 e y(0) = y 0 , se substituirmos t = 0 nessa fórmula obtemos u(0,y 0 ) = C (verifiquem!). 
Comparando com a condição inicial u(0,y 0 ) = e y °, concluímos que C = e y °. Portanto, 


u(x(t), y(t)) = u(t, y 0 ) = y 0 te sinc + e*> +sinit . 


Em x = t e y = y 0 , temos u(x, y) = xye smx + e - y+smx . 


(7.3) Nada obriga a que a condição inicial seja dada em pontos da forma (0,y 0 ). Vejamos 
o que se passa com 

= 1 + sin(2y), u(-y 2 , y) = -y 2 sin(2y). 

dx 

Como queremos que ff = ff ' ff + ff • fy seja |^, escolhemos ff = 1 e ff = 0. Isto 
leva-nos a x(t) = x 0 + t e y(t) = y 0 . A equação original fica 

du ? ? 

— = 1 + sin(2y 0 ), i í(-y 0 ,y 0 ) = -y 0 sin(2y 0 ). 

Como a condição inicial é dada em pontos da forma (— y^, y 0 ), vamos escolher (x(0), y(0)) = 
(x 0 ,yo) = (—yg,yo), ou seja, x(t) = t — y^. Desta forma, em t = 0 ficamos com 

u(x( 0), y (0)) = u(x 0 , y 0 ) = u(-y 2 , y 2 ) = -y 2 sin(2y 0 ). 

A solução de ff = 1 + sin(2y 0 ) é u = t + t sin(2y 0 ) + C, onde C é constante em relação 
a t (mas pode depender de y 0 ). Para determinar o valor de C, substituímos t = 0 na 
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expressão que acabámos de encontrar para u, obtendo u(x(0),y(0)) = u(x 0 ,y 0 ) = C. Como 
a condição inicial era u(x 0 ,y 0 ) = —ygSÍn(2y 0 ), concluímos que C = —ygSÍn(2y 0 ) e que 
ti = t + t sin( 2 y 0 ) - y 2 sin( 2 y 0 ) = t + (t - y 2 ) sin( 2 y 0 ). 

Mas queremos escrever u só em termos dei ey (isto é, sem usar t ou y 0 ). Recordemos 
as parametrizações x{t) = t—y^ e y{t) = y 0 . Se as resolvermos em ordem a tey 0 obtemos 
yo = y et = x + Jo =x "*” y 2 - Dessa forma, 

u(x,y) = t + (t-y 2 ) sin( 2 y 0 ) = x + y 2 + x sin( 2 y). 

(7.4) Exercício. Encontrem a solução u(x,y) de cada uma destas equações. 

(a) ff =y2 - sinx ’ u (°>y) = 2; 

(b) = 2 e 2x + 2x, u( 0 ,y) = 1 + y; 

(c) = 2 u + cosy, u( 0 ,y) = cosy; 

(d) f^ = 8 x-y, u(x,x) = 4x 2 ; 

(e) yg = e* + 3y, u(0,y) = l/y (com y > 0); 

(0 2^x =ysin^x, ii(7t 2 ,y) = y. 

(7.5) Consideremos agora 

du du , /0 

3— + 2— = ií + ye y/ , u(x, 0 ) = x. 
dx dy 

Neste caso, queremos que ^ = f/ ' §7 + 77 ' seja 3|^ + 2|y, pelo que precisamos 
de ^ = 3 e ^ = 2. Além disso, queremos que o ponto inicial (x(0), y(0)) seja da forma 
(x 0 ,0). Para tal, escolhemos x(t) = x 0 + 3t e y(t) = 2t. 



Usando estas curvas características e tendo em conta que ^ = 3|^ + 2|y e que y = 2t, 
a equação original pode ser expressa assim: 

du 

— = u + 2 te, u(x 0 , 0 ) = x 0 . 
dt 

Esta equação é linear e e -t é um fator integrante. Resolvendo-a, obtemos u = t 2 e t + Ce c . 
Em t = 0, temos u(x(0), y(0)) = C. Pela condição inicial, u(x(0),y(0)) = u(x 0 ,0) = x 0 , 
logo C = x 0 . Substituindo este valor, obtemos u = t 2 e f + x 0 e c = e l (t 2 + x 0 ). 

Falta-nos exprimir u apenas em termos de x e y. Resolvendo x = x 0 + 3tey = 2t em 
ordem a t e a x 0 , vemos que t = |y e que x 0 = x — 31 = x — |y. Assim, 

u(x,y) = e t (x 0 + t 2 ) = e y/2 (x - |y + |y 2 ). 
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(7.6) Exercício. Determinem as soluções u(x,y) da equação 

du dl l 

3— + 2— = /(x,y,u), u(x,0) = x, 

dx dy 

onde / é cada uma das seguintes funções: 

(a) /(x,y, ti) = y; (b) /(x,y, u) = 4u; (c) /(x, y, u) = uy; (d) /(x,y,u) = x-y 2 . 

(7.7) Se as curvas características forem parametrizadas por x(t) = t + x 0 e por soluções 
y(x) de equações y' = F(x,y) em ordem a x, teremos = 1 e ^ = F(x,y). Assim, 

du dx du dy du du du 

d7 = dTãt + dr^ = ãI +Ffcy) 'ã? 

Por exemplo, para a equação 

du du 

— ~2y— = e x + y + u, u(0,y) = y. 
dx dy 

Para as características serem soluções de^- = le^j = —2y passando por pontos da forma 
(0,y), escolhemos x(t) = t e y(t) = e~ 2t y 0 . (De onde vem este y(t)?) Escrevendo t e y 0 
em ordem axej, temos também t = x e y 0 = e 2t y = e 2x y. 

Já a equação diferencial fica 

du 

— = e +y + u = e +e y 0 + u, uÇO ,y 0 ) = y 0 . 

A equação é linear, e usando o fator integrante e -t obtemos e~ l u = t — |e _3t y 0 + C. 
Substituindo em t = 0 e usando a condição inicial, obtemos u(0,y 0 ) = — ^y 0 + C = y 0 . 
Logo, C = -y 0 . Resolvendo em ordem a u, e depois substituindo as expressões para t e y 0 
em ordem axey, chegamos então a 

u = te t - |e- 2t y 0 + Ce f = te f - |y 0 e" 2t + f y 0 = xe* - f y + f ye 3 *. 

(Veem o que aconteceu nesta última passagem?) 

(7.8) Exercício. Encontrem as soluções u(x, y) das seguintes equações. 

(a) §f-2yg = 2x-2y, u(0,y) = l+2y; 

(t>) È ~ 2y §^ = y ~ 2x y> u (°>y) = j; 

(c) Ê + x §f = y ’ u (°>y) = °; 

(d) §f + xg = x-2y, u(0,y) = l; 

(e) g + xg = 2 ux, u(0,y) = e y \ 

(f) |^ + cosx|y = x-cosx, u(0,y) = y; 

(g) ff + cosx g = 1, u(0,y) = y 2 ; 

(h) + c °sx = u, u(0,y) = -2y; 

® fl + ^ly = cosx > u (0,y) = cosy. 
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(7.9) Exercício * Vamos identificar condições nas quais a equação 

du du 

— + F(x,y)- — = /(x,y,u) 
ox ay 

(com alguma condição inicial) tem solução única u : El —* E : x, y >-» u(x,y), onde El é 
um aberto conexo por arcos. (Não vamos tentar determinar se a função u assim obtida é 
diferenciável.) Em cada alínea, mantemos as hipóteses das alíneas anteriores. 

(a) Suponham que F é contínua em El e localmente Lipschitz em relação a y. Mostrem 
que, para cada (x 0 ,y 0 ) e El, há exatamente uma solução de 


V(t) = i, 

/(t) = F(x(t),y(t)), 


x(0) = x 0 , 

y(o) = y 0 , 


com todos os seus pontos em El e que não pode ser prolongada para um intervalo maior 
sem sair de El. 

(b) Mostrem que se duas curvas como as das alínea anterior (mas com pontos iniciais 
distintos) se intersetam, então são a mesma curva (mas com uma parametrização diferente). 

(c) Suponham que / é contínua, e localmente Lipschitz em relação a u. Dada uma curva 
(x(t),y(t)) como na alínea (a), mostrem que a equação correspondente 

— =/(x(t),y(t),u) 

tem uma única solução para cada valor inicial tz(x 0 , y 0 ). 

(d) Queremos especificar condições de fronteira ao longo de um curva em El. O que é 
preciso exigir dessa curva para garantir que as curvas características não se intersetam? 

(e) O que é preciso exigir do conjunto El para garantir que cada ponto de El está nalguma 
curva característica? Nessas condições, como obter o valor de u em cada ponto de El? 


(7.10) Também podemos ter situações em que nenhum dos coeficientes é constante. Veja¬ 
mos o que se passa com 

du du ~ 

íx + y) — + y — = x + 2xy + 3y 2 , u(x,l) = 2x. 
ox oy 

Como queremos que ^ = ~j 7'§7 + fj"'f^ = ( x + .y)§ 7 +y§y e que cada característica 
passe por um ponto da forma (x 0 ,1), procuramos soluções de 

jx' = x + y, Jx(0) = x 0 , 

(/= \y(o) = i. 

Escrevendo na forma matricial, como fizemos no episódio [3] temos 


X 

/ 

' i i ‘ 

X 


' x(0) ' 


*0 

_ y _ 


0 1 

_ y _ 

y 

. y(°) . 


1 


cuja solução (usando a exponencial da forma de Jordan) é 


x(t) 


e t te í 

*0 


e l x 0 + te 1 

. y(0 . 


0 e f 

1 


e l 


89/99 




















João Pedro Boavida, Episódios de Equações Diferenciais, 2015 


Tendo em conta que ^ = (x + y) + y |y, a equação original fíca agora 

— = (e f x 0 + te c ) + 2(e t x 0 + te t )e t + 3e 2t , u(x 0 ,1) = 2x 0 . 

dt 

Integrando por partes, obtemos u = e l x 0 + te c — e l + e 2t x 0 + te 2t + e 2t + C. Substituindo 
t = 0 e usando a condição inicial u(x 0 ,1) = 2x 0 , obtemos ainda 2x 0 = 2x 0 + C, ou seja, 
C = 0. Concluímos assim que u = e l x 0 + te f — e l + e 2t x 0 + te 2t + e 2t = x — y + xy + y 2 
(de onde vem a última igualdade?). 


(7.11) Exercício. Considerem 

du du 

x— + y— = u, u(x,-l) = x. 
dx óy 

(a) Determinem a parametrização da curva característica com x(0) = x 0 e y(0) = — 1. 
Qual o maior intervalo onde a parametrização está definida? 

(b) Que pontos estão em alguma dessas curvas? As curvas intersetam-se? 

(c) Determinem u no maior conjunto possível. 


(7.12) Exercício. Considerem agora 



du 

x—=fíx,y) 

dy 


com condição inicial iz(x, 0) = 0 válida em x > 0. 

(a) Determinem a parametrização da curva característica com x(0) = x 0 e y(0) = 0, 
para x 0 > 0. Qual o maior intervalo onde a parametrização está definida? 

(b) Supondo que/(x,y) = x 2 — y 2 , determinem u em M 2 \ {(0,0)}. 

(c) Supondo agora que /(x, y) = 1, determinem u em M 2 \ {(x, 0) : x < 0}. E possível 
estender u a mais pontos? 


(7.13) A equação de continuidade. Dizer que uma certa quantidade é conservada, é dizer 
que, para cada região V, 
rçr 

densidade = — I I fluxo através da fronteira 

lv JJõV 

(dV é o bordo de V). Se chamarmos p à densidade (da quantidade em causa), e usarmos 
a regra de Leibniz (para trocar a derivada para dentro do integral), o lado esquerdo dessa 
equação fica 


d 

dt 


d 

dt 


densidade = — 
dt 


P = 


dp 

d t 


Se chamarmos F ao vetor que descreve a direcção de movimento (e quantidade de material 
movimentado) em cada ponto, e N ao vetor normal a dV apontando para fora da região, e 
depois usarmos o teorema da divergência, o lado direito da fórmula inicial fica 


fluxo através da fronteira = 


dV 


F • N = 


dv 


divF. 
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Ora, se duas funções contínuas têm sempre os mesmos integrais, independentemente da 
região de integração V, é porque são necessariamente iguais. Obtemos portanto a chamada 
equação de continuidade: 

d p 

— = — divF. 
dt 

(Notem que aqui estamos a usar a divergência de forma ligeiramente diferente do habitual, 
incluindo apenas derivadas em relação às variáveis espaciais.) 


(7.14) A equação do calor pode ser obtida a partir da equação de continuidade, se acei¬ 
tarmos algumas hipóteses (com cuja justificação—e condições de validade—não nos vamos 
preocupar). O nosso objetivo é encontrar uma equação diferencial para a temperatura 
u{t,x,y,z ) em cada ponto (x,y, z) de um objeto e em cada instante t. A equação, é claro, 
só se aplicará no interior do objeto; teremos de especificar (de alguma outra forma) o que 
se passa na fronteira do objeto. 

Como a energia é conservada, usamos a correspondente equação da continuidade. Para 
isso, assumimos que a densidade de energia (térmica) em cada ponto é proporcional à 
temperatura (absoluta) nesse ponto (ou seja, p = au, para alguma constante a > 0) e 
que o fluxo de energia é proporcional ao gradiente da temperatura, mas no sentido oposto 
(ou seja, F = — bgradu, para alguma constante b > 0). Combinando tudo, a equação da 
continuidade fica 

du r d 2 u d 2 u d 2 u\ 

a- = bib, S c^u=b[— 2 + — 2 + — 2 ). 

Para nós, os valores de a e b são irrelevantes, e por isso combinamo-los no mesmo lado da 
equação. Vamos considerar apenas o caso unidimensional. (Como usamos apenas derivadas, 
não precisamos que u seja a temperatura absoluta .) Obtemos portanto 

du _ d 2 u 
dt dx 2 

Mas como resolvemos uma equação deste género? 


(7.15) Resolução de EDPs usando séries de Fourier. Queremos encontrar a função u(t, x) 
com domínio t>0e0 <x<| satisfazendo 

d ii ^ 

— = 4—-x - 17u, u(t, 0) = u(t, f ) = 0, iz(0, x) = 4sin(3x) - 7sin(6x). 

dt dx z J 

(As condições no centro chamamos as condições de fronteira, à última condição chamamos 
a condição inicial. Para evitar ter de definir as derivadas parciais na fronteira, trabalhamos 
apenas com funções definidas numa vizinhança do domínio.) 

Vamos usar a seguinte estratégia: fixando o valor de t, queremos que a função /(x) = 
u(t,x) satisfaça /(O) = iz(t, 0) = 0 e /(|) = u(t, |) = 0. Como vamos precisar de calcular 
a derivada = |^. convém que / seja diferenciável (sem saltos nem cantos). Podemos 
começar por fazer um prolongamento ímpar de / ao intervalo [ — (como /(O) = 0, 
isto assegura que / seja diferenciável em x = 0). Se fizermos um prolongamento periódico 
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(de período ^), também asseguramos que / é diferenciável em x = |. Desta forma, / é 
uma função ímpar de período pelo que tem uma série de senos com o seguinte formato: 


/ o)=XX sin 

k>0 


2nkx 

2n/3 


b k sin(3 kx). 

k>0 



Mas é importante notar que este raciocínio se aplica independentemente para cada t. 
Desta forma, os coeficientes podem perfeitamente depender de t. Ou seja, estamos a lidar 
com uma solução da forma 

u(t,x) = E b k {t) sin(3fcx). 

k> o 

Para encontrar os b k (t), calculamos as derivadas parciais desta série (vamos admitir, 
sem justificar, que isto é válido) e substituímo-las na equação original. Como 

x)=y. b k M sin(3 kx), ^(t, x) = V b k (t) (-3 2 k 2 ) sin(3 kx) 

dt ti 5x ti 

(porquê?), a equação original fica 

^ b' k (t) sin(3fcx) = 4^ b k (t) (—3 2 k 2 ) sin(3fcx) — 17^ b k {t) sin(3fcx). 
k> o fc>o fc>o 

Recordando que os senos formam uma base para as funções que estamos a considerar (fun¬ 
ções ímpares de período =y), reconhecemos que os seus coeficientes não se misturam. Por 
isso, para cada k > 0, temos 

b' k (t ) = -4 • 3 2 k 2 b k (t ) - 17b fc (t) = (-36fc 2 - 17) b fc (t). 

Podíamos usar aniquiladores para resolver esta equação, mas este caso é tão simples que 
sabemos imediatamente que a solução é b k (t) = e^~ 36k ~ 17 ^ c b k (0). Assim, a solução da 
equação original é 

u(t, x) = ^ g(-36fc 2 -i7) t bksin ( 3fcx )_ 
k> 0 

Para determinar os valores de b k ( 0), substituímos t = 0 na solução geral, obtendo 

u(0,x) = E b fc (0) sin(3fcx). 

k>0 

Neste ponto, pode acontecer uma de duas coisas: ou a condição inicial já vem expressa 
neste tipo de série (caso em que lemos os coeficientes diretamente), ou não vem (caso em 
que temos de calcular os coeficientes da condição inicial—prolongada da mesma forma que 
os u(t,x)—numa série deste tipo, usando os produtos internos). 

No nosso caso, a condição inicial iz(0, x) = 4 sin(3x) — 7 sin(6x) já vinha expressa como 
uma série do tipo procurado. Assim, podemos ler diretamente que fq(0) = 4, b 2 (0) = —7, 
e que nenhum outro k ocorre. Ou seja, a solução final da nossa equação é 

u(0, x) = 4e ( - 36 “ 17)t sin(3x) - 7e ( - 36 ' 22 “ 17)t sin(6x). 
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(7.16) Se as condições de fronteira fossem |^(t, 0) = |^(t, |) = 0 e, como antes, /(x) = 
u(t,x) para algum t fixo, então o prolongamento compatível seria um prolongamento par. 
Assim, teria uma série de cossenos 

/(x) = u(t,x) = E a fc (t) cos(3fcx). 
k>0 



(7.17) Se as condições de fronteira fossem iz(t,0) = |^(t, |) = 0 e, como antes, /(x) = 
iz(t, x) para algum t fixo, então o prolongamento compatível seria um prolongamento ímpar 
(para evitar ter um canto em / (O) = 0). Porém, não podemos fazer um prolongamento de 
período pois corremos o risco de ter saltos no outro extremo do intervalo. 



Para evitar este problema, fazemos um novo prolongamento, até ao intervalo [ — ], 
e depois consideramos uma função de período Dessa forma, a série fica 


/(x) = u(t,x) = ^]b fc (t) 
fc> o 


sin 


2nkx 

4n/3 


^—i 3kx 

2j sin—. 

fc>o Z 



(7.18) Exercício. Como ilustrado na figura, a função / tem a seguinte simetria: /(x + = 

—/(x). Usando o exercício Q6.19D , mostrem que t) = 0 se k é par. Verifiquem que o 
mesmo se passa independentemente do comprimento do intervalo original. Concluam que 
para este tipo de prolongamento a série de senos inclui apenas k ímpar. 

(7.19) Se as condições de fronteira fossem |^(t, 0) = u(t, |) = 0 e, como antes, /(x) = 

u(t,x) para algum t fixo, então o prolongamento compatível seria um prolongamento par 
(para evitar ter um salto em / (O) = 0). Tal como para o caso anterior, não podemos fazer 
um prolongamento de período desta vez porque corremos o risco de ter cantos no 
outro extremo do intervalo. Fazemos então um novo prolongamento a [ — e depois 
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consideramos uma função de período ff. Seguindo um raciocínio semelhante ao de antes, 
chegamos a 

^—i 3kx 

/(x) = u(t,x) = a fc(t) cos—• 

fc >0 z 

k ímpar 


(7.20) Exercício. Para esta função, /(x + ff ) = —/(x). Usem o exercício Q6.24P para mos¬ 
trar que a fc (t) = 0 se k é par. Verifiquem que o mesmo se passa independentemente do 
comprimento do intervalo original. Concluam que para este tipo de prolongamento a série 
de cossenos inclui apenas k ímpar. 

(7.21) Exercício. Considerem a EDP 

du d 2 u 

"ÕT "õ 2 
Ot OX Á 

onde iz(t,x) está definida para t > 0 e 0 < x < 2n. Encontrem a solução satisfazendo 
as condições de fronteira e inicial listadas em cada uma das alíneas (ou seja, cada alínea 
corresponde a uma solução distinta). 

(a) iz(t,0) = u{t,2n) = 0, u(0,x) = sinx; 

(b) u(t,0) = iz(t,27i) = 0, u(0,x) = 3sinx + sin(2x); 

(c) u(t, 0) = u(t, 27i) = 0, u(0,x) = sin(x/2); 

(d) u(t,0) = u(t,27i) = 0, u(0,x) = cosx; 

(e) u(t,0) = u(t, 27i) = 0, u(0,x) = x; 

(f) ff (t, 0) = ff (t, 2 ti) = 0, u(0, x) = 3 cos(2x); 

Cg) ff(t,0)= ff(t,27i) = 0, u(0,x) = sinx; 

(h) |f(t,0)=ff(t,27r) = 0, u(0,x) = x; 

(i) ff (t, 0) = ff (t, 2n) = 0, u(0, x) = H 0 (x) - fí^x); 

(j) u(t, 0) = ff (t, 2n) = 0, u(0, x) = sin(2x); 

(k) u(t,0) = ff (t, 271) = 0, u(0,x) = 3sin(2x) - 5cos(4x); 

(l) u(t,0) = ff (t,2n) = 0, u(0,x) = x; 

(m) iz(t, 0) = ff (t, 2n) = 0, u(0,x) = l; 

(n) ff (t, 0) = iz(t, 2n) = 0, u(0, x) = 1; 

(o) ff (t, 0) = iz(t, 27t) = 0, u(0,x) = cos(2x) + 5 cos(4x); 

(p) ff (t, 0) = iz(t, 2n) = 0, u(0, x) = cos(3x). 
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(7.22) Exercício. Consideramos agora uma variante com coeficientes não constantes: 

du , d 2 u 

—~(t,x) = t 3 —- + 5u, 
dt dx z 


com í>OeO<i<|e com condições de fronteira |^(t, 0) = |^(t, -|) = 0. Qual a 
solução geral, escrita como série de cossenos? Qual a solução com condição inicial u(0, x ) = 
n — e cos(2x)? 


(7.23) Exercício. Quais as soluções, com domínio t>0e0<x<2ne com condições de 
fronteira u(t, 0) = u(t, 2 n) = 0 e condição inicial u(0, x) = sin |, das seguintes equações? 


du d 2 u 

(a) trr = 8-r— 2 +(smt)u; 
dt dx Z 



4 õ 2 u 
t + 3 dx 2 


—2 tu; 


du d 2 u 

(c) — = — r+u-smx. 
dt dx z 


(7.24) A equação das ondas aparece em contextos diferentes. Em dimensão (espacial) 
1, com variáveis ter, podemos pensar que u descreve o deslocamento (por exemplo, de 
pontos num material elástico) em relação à posição de equilíbrio, e que o momento é a 
quantidade conservada. Para um material de densidade constante, é razoável dizer que o 
momento é a |^, onde a > 0 é uma constante, e aproximar a variação na distância entre 
pontos do material por |^. Como a transferência de momento corresponde a uma força e 
a força elástica é proporcional à diferença entre essa distância e o seu valor de equilíbrio, 
escrevemos F = — onde k > 0. Portanto, de acordo com a equação da continuidade 


' 3 
d du 


(I7.13D, teríamos = c 2 4~^r, onde c 2 = - é uma constante positiva. 

-■ dt dt dx dx a r 


A forma adequada de modelar o problema em dimensão mais alta depende do pró¬ 
prio problema. Para os nossos fins, vamos simplesmente admitir que a nossa análise em 
dimensão 1 se mantenha válida ao longo de cada curva tangente a F = — kgradiz, e que 
chegaríamos sempre a algo como (exemplificando em dimensão 3) 


d 2 U y y y f d 2 U d 2 U 

T = c divgrad u = c lap u = c ^ 

dt z v dx z dy 


d 2 u a 

Jz 2 )' 


Os próximos exercícios exploram variantes desta equação em dimensão 1 (ou 2, contando 
também a variável t). Como vamos ver, c pode ser vista como a velocidade das ondas. 


(7.25) Exercício. Queremos encontrar a solução u(t,x), com t>0e0<x<n, da EDP 

d 2 u d 2 u 

— - y = —- y d" 3lí 

dt 2 dx 2 

com condições de fronteira u(t, 0) = u(t, 7t) = 0 e condições iniciais iz(0,x) = sinx — 
3sin(2x) e ||(0,x) = 4sin(2x). 

(a) Comecem por mostrar que a série adequada é u(t, x) = Xifc>o sin(fcx). 

(b) Mostrem que a equação satisfeita por b k (t) é ( D 2 — 3 + k 2 ) b k = 0 (usando D para 
indicar a derivada em ordem a t). 
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(c) Para cada k > 0, determinem b k { 0) e b' k { 0). Em particular, decidam quais os valores 
de k que não são relevantes para o resto do problema. 

(d) Para cada valor relevante de k, determinem a solução b k da equação da alínea (b) 
com as condições iniciais indicadas na alínea (c). Notem que valores diferentes de k levam 
a tipos diferentes de solução. 

(e) Juntem toda a informação das alíneas anteriores e obtenham uma expressão para a 
solução iz(t,x). 

(7.26) Exercício. Encontrem a solução u(t,x), com t>0e0<x<7t, da EDP 

d 2 u d 2 u 
——r = ——~ + 4u 

dt 2 dx 2 

com condições de fronteira iz(t, 0) = u(t, n) = 0 e condições iniciais iz(0, x) = 2sin(2x) — 
3 sin(3x) + 4sin(4x) e ff (0, x) = 3 sin(2x). 

(7.27) Exercício. Considerem a EDP 

d 2 u _ d 2 u 
d t 2 dx 2 ’ 

onde iz(t,x) está definida para t>0e0<x< 2n. Em cada alínea, encontrem a solução 
satisfazendo as condições indicadas. 

(a) u(t,0) = u(t, 2n) = 0, u(0,x) = 3sinx + sin(2x), ff(0,x) = 0; 

(b) u(t,0) = u{t,2n) = 0, u(0, x) = 0, ff(0, x) = sin(x/2); 

(c) u(t,0) = u(t,27t) = 0, u(0,x) = 3sinx + sin(2x), = sin(x/2); 

(d) iz(t, 0) = u(t,2n) = 0, u(0, x) = x, |^(0, x) = 0; 

(e) g(t ) 0) = g(t,2n) = 0 ) u(0,x) = 0, ff (0,x) = 3cos(2x); 

(f) |f(t,0)=ff(t,27i) = 0, u(0, x) = 0, §f(0,x) = 7t; 

(g) ff(t,0)= ff(t,27i) = 0, u(0,x) = cosx, ff (0,x) = x; 

(h) u(t, 0) = ff (t, 27t) = 0, u(0,x) = sin(3x/4), ff (0,x) = f sin(x/4); 

(i) u(t, 0) = ff (t, 20 = 0, u(0,x) = 0, ff (0,x) = 1; 

0') ff(t,0) = u(t,2O = 0, u(0,x) = —6cos(x/4), ff(0,x) = 0; 

(k) ff(t,0) = u(t,2O = 0, iz(0,x) = 0, ff(0,x) = —6cos(x/4); 

(l) ff (t,0) = u(t, 20 = 0, u(0,x) = ff (0,x) = —6cos(x/4). 

(7.28) Exercício. Encontrem a solução u(t,x), com t>0e0<x<f, da EDP 

d 2 u d 2 u 

—r = - 3 cosx sm(2t), 

ot A ox z 

com condições de fronteira u(t, 0) = ff (t, f) = 0 e condições iniciais u(0,x) = sin(3x) e 
ff (0, x) = 2 cosx. 

(7.29) Exercício. Encontrem a solução geral u(t,x), com t>0e0<x<l, da EDP 

d 2 u 2 d 2 u 
dt 2 dx 2 
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(c > 0 é uma constante) com condições de fronteira iz(t, 0) = iz(t, 1) = 0. 

(a) Considerem primeiro condições iniciais da forma iz(0,x) = Asin(fx) e |^(0, x) = 
B sin(fx). Quais os valores razoáveis para t? 

(b) Fixem iz(0,x) = Asin(fx) e |^(0,x) = 0. Qual a solução? Escrevam-na com base 
em sin(f ct + tx) e sin(£ct — tx). Qual diriam que é a velocidade dessas ondas? 

(c) Fixem agora tz(0, x) = 0 e ^j(0,x) = B sin(fx). Qual a solução? Escrevam-na com 
base em cos(fct + tx) e cos(£ ct — tx). Qual diriam que é a velocidade dessas ondas? 

(d) Qual a solução geral da EDI) escrita em série de senos? 


(7.30) Exercício. Considerem de novo uma solução iz(t,x) de 

d 2 u _ 2 d 2 u 
~dt 2=C fa 2 ’ 

onde c > 0 é uma constante. Vamos usar a mudança de variável a = x — ct e b = x + ct. 

(a) Escrevam t e x em função d e a e b. Usem essas expressões para escrever u em 
função deaeb. 

(b) Calculem e mostrem que é igual a 0. Concluam que é uma função apenas 
de a, e que u = f(a) + g(b) = /(x - ct) + g(x + ct). 

(c) Mostrem que se iz(t, 0) = 0 para todo o t > 0, então g(x) = —/(—x). Mostrem que 
se iz(t,0) = iz(t,() = 0 para todo o t, então / e g têm período 2 1. Conseguem encontrar 
outras propriedades de / e g quando u satisfaz os vários tipos de condições de fronteira? 

(d) Se h for uma primitiva de ^f^(0,x) com fr(0) = g(0) — /(0), mostrem que 

u(0,x) — h(x) u(0, x) + h(x) 

f 0) =-^- e g(x) =---. 

(e) Concluam que 


ií( 0, x + ct) + iz(0,x — ct) 1 

u(t,x) = -1- 

2 2c 


/-x+ct 


x—ct 


ÕU 

dt 


(0,s)ds. 


(7.31) A equação de Laplace aparece frequentemente quando lidamos com potenciais de 
campos. Um caso típico é o campo eletrostático, em que, na ausência de carga, o potencial 
é uma função harmónica, ou seja, é uma solução de 

d 2 u i d 2 u i d 2 u ^ 

Jx 2+ J^ 2 + ~dz 2= °' 

Em duas variáveis, sabemos que a parte real e a imaginária de uma função holomorfa são 
harmónicas. Também podemos usar séries de Fourier para encontrar soluções. 

Suponham que u(x,y) define uma função harmónica em 0<x<2e0<y<le com 
condições de fronteira u(x,0) = 0, u(x,l) = (e n — e~ n ) sin(7ix), u(0,y) = Se^sm^ y), 
u{2,y) = 3e~ n sin(7ty). Com os métodos que vimos até aqui, não é possível exprimir di¬ 
retamente a solução de uma tal equação. Com efeito, o que sempre fizemos foi escolher a 
variável com um domínio limitado e escolhido o tipo de série com base nas condições de 
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fronteira relativas a essa variável. Ora, nenhuma destas condições de fronteira é compatível 
com as séries que vimos. 

Felizmente, a equação de Laplace é linear, portanto podemos escrever u = + u 2 , 

onde ií 1 e u 2 também satisfazem a equação de Laplace, mas cada uma delas tem as condi¬ 
ções de fronteira relativas a uma variável. Por exemplo, para Ui impomos u 1 (v,0) = 0, 
u 1 (x, 1) = (e n — e ~ K ) sin(7tx), e lí^O.y) = u 1 (2,y) = 0, enquanto para u 2 impomos 
u 2 (x,0) = u 2 (x, 1) = 0, u 2 (0,y) = 3e" sin(7ty), e u 2 (2,y) = 3e _7r sin(7tj). 


(7.32) Exercício. Vamos determinar u a partir de u 1 e u 2 . 

(a) Resolvam a equação para u 1 , depois de mostrar que a série adequada é 


tíi íx,y) = ^ b áy) sin 


nkx 


k> 0 


(b) Resolvam a equação para u 2 , depois de mostrar que a série adequada é 

u 2 0,y) = Xi kfcM sin(Ttfcy). 


k>0 


(c) Verifiquem que Ui + u 2 satisfaz as condições de fronteira de u e determinem u. 


(7.33) Exercício. Considerem a EDP 

d 2 u d 2 u _ 
dx 2 + dy 2 

onde u(x,y) está definida para 0<x<27te0<y<7t. Em cada alínea, encontrem a 
solução satisfazendo as condições indicadas. 

(a) u(x, 0) = 3 sin(2x), iz(x, 7t) = 3e 2n sin(2x), u(0, y) = siny, u{2n,y) = e~ 2n smy, 

(b) u(x,0) = -e“ 4 " cos(4x), u(x, n) = - cos(4x), (0,y) = |^(27t,y) = 0; 

(c) u(x, 0) = u(x, ti) = 0, f^(0,y) = sin(2y), |^(27t,y) = e 4n sin(2y); 

(d) u(x, 0) = 0, u(x, 7t) = [e 2n — 1) cosx, |^(0,y) = sin(2y), |^(2 ti, y) = e 4 " sin(2y); 

(e) u(x, 0) = 0, u(x, 7t) = 3(e" — 1) cos(x/2), §^(0,y) = |^(27t,y) = 0; 

(í) u(x,0) = 0, u(x, n) = x - n, §^(0, y) = j^[2n,y) = 0; 

(g) u(x,0) = 1, u(x,7i) = 0, u(0,y) = |^(27r,y) = 0. 

(7.34) Exercício. Consideramos agora, para uma constante A, a EDP 

d 2 u d 2 u _ 
dx 2 + dy 2 

com condições de fronteira u(x,0) = u(x, 2) = u(0,y) = u(n,y) = 0, onde 0 < x < n e 
0 < y < 2. Mostrem que se u ^ 0, então A = —k 2 — yí 2 , onde k e i são inteiros positivos, 
eué um múltiplo de x,y <-* sin(fcx) sin ( ~y. 
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(7.35) Respostas de exercícios. 


Q7.4D (a) cosx + xy 2 + 1. (b) e 2x +x 2 + y. 
(c) —icos y + |e 2 *cosy. (d) 4x 2 — xy + y 2 . 
(e) (e*/y) + 3x. (f) — y cos y/x. 

m (a) x - |y + ±y 2 . (b) e 2y (x - |y). 

(c) e y2/4 (x - |y). 

(d) x-iy + ixy-|y 2 -|y 3 . 


( [7725] ) (c) b 1 (0) = 1; b 2 (0) = -3; b'(0) = 4. 

(d) b 1 (t)=i(e' /2t + e-' /2t ); 

= 4sint — 3cost. (e) 

^(e'^ 2t + e~^ 2t ) sinx + (4sint —3cost) sin(2x). 


( |7.26P (2 + 3t)sin(2x) — 3cos(\/5t)sin(3x) + 
4cos(\/l2 t) sin(4x). 


d7.8D (a) x 2 + y +e 2x y + 1. (b) xy + e 2x y. 
(c) xy - |x 3 . (d) |x 2 + |x 3 - 2xy + 1. 

(e) exp(|x 2 +y). (0 /x 2 — 2sinx+y. 

(g) x + (y — sinx) 2 . (h) 2sinx e x — 2 ye x . 

(i) sinx + cos(y/(l + x 2 )). 


( |7. 1 lD (a) t •—> (x 0 e f , — e 7 ). R. 

(b) {(x,y) e® 2 :y<0j. não. (c) x. 


( |7.12| ) (a) t •—> (x 0 cos t,x 0 sin t). R. (b) xy. 

(c) u(r cos 6, r sin 6) = 6 se r > 0 e — n < 9 < n. 
não. 


( |7.21| ) (a)e 4t sinx. (b) 3e 4t sinx + e 7 sin(2x). 
(c) e 19t ^ 4 sin(x/2). (d) ^ yy e (5_ T )7 sin y. 


n> 0 
n ímpar 

i (-l)" ( 5 -^)t „■ nx 


(e) S„ >0 -4 „ e 5 4 sin y. (f) 3e 7 cos(2x). 

n ímpar H 


(h) íte 57 + X]n>o —^e (5 4 )7 cosy. 

n ímpar nn 2 2 

® è e5t + S n >o^sin2 e (5 -T^cosf. 


2tz 2Lm>0 nn 2 

sin y 

n 2 —64 2 


®Sn>o f^sin^^-^^sin^. 

n ímpar r 

V f 96/tt • nn 20n/n \ nx 

ck, Ç£7^ sm ““^) e “ s,n -' 

COE„>o 


n impar 


(m) 2 n >0 ye (5 I6 }t sin^. 

n impar 

(n) Xj?>° ^ sin f e(5 ”^ )tcosí 7 - 

n ímpar nn z ^ 

a V" 1 f 4n/n . 20n/n \ • nn — nx 

(0) Çs;„tó + ^ )sm " e “ cos " 

CP)S;>« 4ntsmfe®-á)*cosf. 

n impar n z ^ 


H7.22D Sn >0 e5t " E a nW cos(2nx). 

7 ie St — e 1+5t_t4 cos(2x). 

( |723l ) (a) e 1_cost_2t s in (b) -f-e- 72 sin§. 


(c) e 3í / 4 sin | — t sinx. 


( |7.27P (a) 3sinx cos(3t) + sin(2x) cos(6t). 


(b) | sin | sin y . (c) 


3 sin x cos(31) + sin(2x) cos(6t) + | sin | sin y. 

(d) Z„>o -4^ sin f cos? f- 

(e) i cos(2x) sin(6t). (f) nt. (g) 
cosx cos(3t) + 7i t — Xi">o yy cos y sin ^y. 
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(h) sin f sin ^ + sin 77 cos 7- 

(i) Sn >0 yy sin y sin yi. 

n impar J7tn H H 

(j) —6cos | cosy. (k) —8cos| siny. 

( 1 ) —6 COS 7 COS 77 — 8 COS 7 sin 77. 

4 4 4 4 


( |7.28[ ) sin(3x) cos(3t) + cosx sin(2t). 


d7.29D (a) £ = mi, com n e Z. 

(b) Acos(£ct) sin(£x) = 

^ sin(£ct + £x) — | sin(£ct — £x). 

(c) £ sin(£ ct) sin(£x) = 

y; COS(£ct - £x) - y; C0S(£ct + £x). 

(d) (A„cos(n 7 ict) + — sin(n 7 i C t)) sin(ri 7 tx). 


( TÃ3ÜD (a)t=^;x=3±‘;u(^ > 5±6). 


d7.32P (a) e 7ly sin(7ix) — e - " 3 ' sin(7ix). 

(b) Se 71-71 ^ sin( 7 iy). 

(c) e 7137 sin( 7 tx) — e -713, sin( 7 ix) + Se 71-7177 sin( 7 iy). 

d7.33P (a) 3e 2y sin(2x) + e~ x siny. 

(b) —e 4y_471 cos(4x). (c) ^e 2x sin(2y). 

(d) e n+y cosx — e K ~ y cosx + |e 2x sin(2y). 

(e) 2e n ^ 2 (e y / 2 - e“ y / 2 ) cos(x/2). 

rc\ y_8_ exp(ny/2)-exp(-ny/2) nx 

n ímpar 717,2 exp(n7r/2)-exp(-nji/2) 2 

V 4 exp(n(jr—y)/4)—exp(n(y—70/4) nx 

(gJ Zj'?>0 “- 


n impar 


_ 

exp(n7r/4)—exp(— nn/4) 4 
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